Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,328,408 (T)7→8 coding (1390/1650 nt) yfaQ ← tandem DUF2300 domain‑containing protein YfaQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,328,4011.T93.3% 46.3 / ‑3.1 15coding (1397/1650 nt)yfaQtandem DUF2300 domain‑containing protein YfaQ
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (11/3);  total (11/4)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGAGACTAAACGGCAAACCGCAAACACGTCTGGTTCGTTGTACCCCGTCTCAGCTTGTAATATACGCCGCCACTGCGGCATTTTTTTCGCCAGCCAGGCTTTCGCTTTGGGTAATAACTGGCAGGTTGAGCGTGGCGCTCCCCAGCGGCTAAGAC  >  NC_000913/2328321‑2328476
                                                                                 |                                                                            
cTGAGACTAAACGGCAAACCGCAAACACGTCTGGTTCGTTGTACCCCGTCTCAGCTTGTAATATACGCCGCCACTGCGGCa                                                                               <  2:366646/81‑1 (MQ=255)
   agaCTAAACGGCAAACCGCAAACACGTCTGGTTCGTTGTACCCCGTCTCAGCTTGTAATATACGCCGCCACTGCGGCAttc                                                                           <  2:366647/81‑2 (MQ=255)
        aaaCGGCAAACCGCAAACACGTCTGGTTCGTTGTACCCCGTCTCAGCTTGTAATATACGCCGCCACTGCGGCAtttttttt                                                                     >  2:366661/1‑80 (MQ=255)
         aaCGGCAAACCGCAAACACGTCTGGTTCGTTGTACCCCGTCTCAGCTTGTAATATACGCCGCCACTGCGGCAttttttttc                                                                    <  1:366648/81‑3 (MQ=255)
         aaCGGCAAACCGCAAACACGTCTGGTTCGTTGTACCCCGTCTCAGCTTGTAATATACGCCGCCACTGCGGCAttttttttc                                                                    >  1:366675/1‑79 (MQ=255)
            ggCAAACCGCAAACACGTCTGGTTCGTTGTACCCCGTCTCAGCTTGTAATATACGCCGCCACTGCGGCAttttttttcgc                                                                  >  1:366676/1‑76 (MQ=255)
                  ccGCAAACACGTCTGGTTCGTTGTACCCCGTCTCAGCTTGTAATATACGCCGCCACTGCGGCATTTTTTTTCGccagccag                                                            >  1:366672/1‑81 (MQ=255)
                     cAAACACGTCTGGTTCGTTGTACCCCGTCTCAGCTTGTAATATACGCCGCCACTGCGGCATTGTTATTCGCCAGCCAGGCt                                                         <  2:366649/81‑1 (MQ=255)
                     cAAACACGTCTGGTTCGTTGTACCCCGTCTCAGCTTGTAATATACGCCCCACTGCGGCATTTTTTTTCGCCAGCCAGGCtt                                                        >  1:366679/1‑81 (MQ=255)
                       aaCACGTCTGGTTCGTTGTACCCCGTCTCAGCTTGTAATATACGCCGCCACTGCGGCATTTTTTTTCGCCAGCCAGGCttt                                                       >  2:366683/1‑81 (MQ=255)
                               tGGTTCGTTGTACCCCGTCTCAGCTTGTAATATACGCCGCCACTGCGGCATTTTTTTTCGCCAGCCAGGCTTTCGCTTTgg                                               >  2:366684/1‑81 (MQ=255)
                               tGGTTCGTTGTACCCCGTCTCAGCTTGTAATATACGCCGCCACTGCGGCATTTTTTTTCGCCAGCCAGGCTTTCGCTTTgg                                               >  1:366680/1‑81 (MQ=255)
                                               gTCTCAGCTTGTAATATACGCCGCCACTGCGGCATTTTTTTTCGCCAGCCAGGCTTTCGCTTTGGGTAATAACTGGCAGGt                               >  1:366673/1‑81 (MQ=255)
                                                          tAATATACGCCGCCACTGCGGCATTTTTTTTCGCCAGCCAGGCTTTCGCTTTGGGTAATAACTGGCAGGTTGAGCGTGgc                     <  2:366650/80‑1 (MQ=255)
                                                           aaTATACGCCGCCACTGCGGCATTTTTTTTCGCCAGCCAGGCTTTCGCTTTGGGTAATAACTGGCAGGTTGAGCGTGgcgc                   >  1:366689/1‑81 (MQ=255)
                                                              ataCGCCGCCACTGCGGCATTTTTTTTCGCCAGCCAGGCTTTCGCTTTGGGTAATAACTGGCAGGTTGAGCGTGGCGCTcc                >  1:366688/1‑81 (MQ=255)
                                                                 cgacgccACTGCGGCATTTTTTTTCGCCAGCCAGGCTTTCGCTTTGGGTAATAACTGGCAGGTTGAGCGTGGCGCTcccc              <  2:366651/77‑1 (MQ=255)
                                                                         cTGCGGCATTTTTTTTCGCCAGCCAGGCTTTCGCTTTGGGTAATAACTGGCAGGTTGAGCGTGGCGCTCCCCAGCGGCTaa     <  2:366652/81‑1 (MQ=255)
                                                                           gCGGCATTTTTTTTCGCCAGCCAGGCTTTCGCTTTGGGTAATAACTGGCAGGTTGAGCGTGGCGCTCCCCAGCGGCTAAGa   >  1:366697/1‑81 (MQ=255)
                                                                           gCGGCATTTTTTTTCGCCAGCCAGGCTTTCGCTTTGGGTAATAACTGGCAGGTTGAGCGTGGCGCTCCCCAGCGGCTAAGa   >  2:366681/1‑81 (MQ=255)
                                                                            cGGCATTTTTTTTCGCCAGCCAGGCTTTCGCTTTGGGTAATAACTGGCAGGTTGAGCGTGGCGCTCCCCAGCGGCTAAGAc  >  1:366686/1‑81 (MQ=255)
                                                                                 |                                                                            
CTGAGACTAAACGGCAAACCGCAAACACGTCTGGTTCGTTGTACCCCGTCTCAGCTTGTAATATACGCCGCCACTGCGGCATTTTTTTCGCCAGCCAGGCTTTCGCTTTGGGTAATAACTGGCAGGTTGAGCGTGGCGCTCCCCAGCGGCTAAGAC  >  NC_000913/2328321‑2328476

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: