Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 600,553 G→A 28.4% L147L (CTG→CTA dmsB → dimethyl sulfoxide reductase, anaerobic, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE600,5530GA28.4% 58.6 / 20.3 56L147L (CTG→CTAdmsBdimethyl sulfoxide reductase, anaerobic, subunit B
Reads supporting (aligned to +/- strand):  ref base G (23/17);  new base A (10/6);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACAT  >  minE/600495‑600598
                                                          |                                             
gttgttATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTgg                                     >  1:2644615/1‑69 (MQ=255)
gttgttATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTgg                                     >  1:2721789/1‑69 (MQ=255)
gttgttATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTgg                                     >  1:782895/1‑69 (MQ=255)
       tGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTc                                >  1:566454/1‑67 (MQ=255)
       tGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCgg                              >  1:1386871/1‑69 (MQ=255)
       tGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCgg                              >  1:1511350/1‑69 (MQ=255)
        gACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCgg                              <  1:1778081/68‑1 (MQ=255)
        gACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCgg                              <  1:1613414/68‑1 (MQ=255)
         aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGc                             >  1:234379/1‑68 (MQ=255)
         aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGc                             >  1:1443125/1‑68 (MQ=255)
         aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGc                             >  1:811681/1‑68 (MQ=255)
              gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc                       <  1:1454588/69‑1 (MQ=255)
              gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc                       <  1:3057055/69‑1 (MQ=255)
              gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc                       <  1:1825213/69‑1 (MQ=255)
              gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc                       <  1:3182408/69‑1 (MQ=255)
              gtTGCGGAAGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc                       <  1:2990846/69‑1 (MQ=255)
                          aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg            >  1:1918571/1‑68 (MQ=255)
                          aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg            >  1:372127/1‑68 (MQ=255)
                          aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg            >  1:937652/1‑68 (MQ=255)
                          aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg            >  1:817901/1‑68 (MQ=255)
                          aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg            >  1:3142960/1‑68 (MQ=255)
                          aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg            >  1:148853/1‑68 (MQ=255)
                          aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg            >  1:1787895/1‑68 (MQ=255)
                          aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg            >  1:3034058/1‑68 (MQ=255)
                              aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa       >  1:2944792/1‑69 (MQ=255)
                              aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa       >  1:2563538/1‑69 (MQ=255)
                              aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa       >  1:638323/1‑69 (MQ=255)
                              aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa       >  1:2137095/1‑69 (MQ=255)
                              aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa       >  1:1677454/1‑69 (MQ=255)
                              aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa       >  1:1654183/1‑69 (MQ=255)
                                ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa     <  1:2293559/69‑1 (MQ=255)
                                ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa     <  1:3245806/69‑1 (MQ=255)
                                ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa     <  1:2517649/69‑1 (MQ=255)
                                ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa     <  1:2056622/69‑1 (MQ=255)
                                ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa     <  1:2052734/69‑1 (MQ=255)
                                ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa     <  1:660962/69‑1 (MQ=255)
                                ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa     <  1:1201173/69‑1 (MQ=255)
                                   aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGCTTTCGGACCTATCGACGAGCTGCGTAAAAAACa   >  1:2233302/1‑68 (MQ=255)
                                   aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGAt                                   >  1:2156433/1‑36 (MQ=19)
                                   aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac  >  1:976327/1‑68 (MQ=255)
                                   aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac  >  1:3232381/1‑68 (MQ=255)
                                   aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac  >  1:3133524/1‑68 (MQ=255)
                                   aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac  >  1:2829044/1‑68 (MQ=255)
                                   aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac  >  1:2440226/1‑68 (MQ=255)
                                   aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac  >  1:2142610/1‑68 (MQ=255)
                                   aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac  >  1:1741115/1‑68 (MQ=255)
                                   aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac  >  1:1206470/1‑68 (MQ=255)
                                     cTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa       <  1:2754409/62‑1 (MQ=255)
                                     cTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa       <  1:1151059/62‑1 (MQ=255)
                                          ttGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa     <  1:3000979/59‑1 (MQ=255)
                                                 ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa     <  1:2793834/52‑1 (MQ=255)
                                                 ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa     <  1:2467035/52‑1 (MQ=255)
                                                 ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa     <  1:2096578/52‑1 (MQ=255)
                                                 ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa     <  1:1734795/52‑1 (MQ=255)
                                                 ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa     <  1:1324400/52‑1 (MQ=255)
                                                 ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa     <  1:1202253/52‑1 (MQ=255)
                                                          |                                             
GTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACAT  >  minE/600495‑600598

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