Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,900,322 C→T 42.9% intergenic (‑49/‑270) ygcW ← / → yqcE putative SDR family oxidoreductase/putative MFS transporter, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,900,3220CT42.9% 1.2 / 13.6 14intergenic (‑49/‑270)ygcW/yqcEputative SDR family oxidoreductase/putative MFS transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base T (2/4);  total (6/8)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

CCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTAT  >  NC_000913/2900188‑2900428
                                                                                                                                      |                                                                                                          
ccGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAg                                                                                               >  4:484451/1‑148 (MQ=255)
 cGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATTGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCTCTCGCGGGGGAAGt                                                                                              <  2:1083152/148‑1 (MQ=255)
                     aTTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTg                                                                          <  3:484451/148‑1 (MQ=255)
                             tttACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAgg                                                                                                   <  3:267033/115‑1 (MQ=255)
                             tttACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAgg                                                                                                   >  4:267033/1‑115 (MQ=255)
                                    ttCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCTCTCGCGGGGGAAGtt                                                                                             <  1:1018100/114‑1 (MQ=39)
                                    ttCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCTCTCGCGGGGGAAGtt                                                                                             >  2:1018100/1‑114 (MQ=39)
                                            gaAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGa                                                                                                 <  1:688035/102‑1 (MQ=255)
                                            gaAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGa                                                                                                 >  2:688035/1‑102 (MQ=255)
                                                                  ttGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCTCTCGCGGGGGa                                                                                                 >  4:606261/1‑80 (MQ=37)
                                                                  ttGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCTCTCGCGGGGGa                                                                                                 <  3:606261/80‑1 (MQ=37)
                                                                      gagaTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCgaga                                                                                                     <  4:550811/72‑1 (MQ=255)
                                                                      gagaTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCgaga                                                                                                     >  3:550811/1‑72 (MQ=255)
                                                                                                                                 ttCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGt    <  1:167069/110‑1 (MQ=39)
                                                                                                                                      ccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCt                                   >  2:864923/1‑74 (MQ=255)
                                                                                                                                      ccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCt                                   <  1:864923/74‑1 (MQ=255)
                                                                                                                                      ccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTAt  <  4:355904/107‑1 (MQ=255)
                                                                                                                                      ccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTAt  >  3:355904/1‑107 (MQ=255)
                                                                                                                                      |                                                                                                          
CCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTAT  >  NC_000913/2900188‑2900428

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: