Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 600,547 T→C 22.1% C145C (TGT→TGC dmsB → dimethyl sulfoxide reductase, anaerobic, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE600,5470TC22.1% 93.8 / 20.0 72C145C (TGT→TGCdmsBdimethyl sulfoxide reductase, anaerobic, subunit B
Reads supporting (aligned to +/- strand):  ref base T (36/20);  new base C (10/6);  total (46/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

TGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACAT  >  minE/600480‑600598
                                                                   |                                                   
tGACCAAATGGGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc                                                    <  1:2575951/69‑1 (MQ=255)
tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc                                                    >  1:104665/1‑69 (MQ=255)
tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc                                                    <  1:984901/69‑1 (MQ=255)
tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc                                                    >  1:732957/1‑69 (MQ=255)
tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc                                                    >  1:725863/1‑69 (MQ=255)
tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc                                                    >  1:375619/1‑69 (MQ=255)
tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc                                                    >  1:3005265/1‑69 (MQ=255)
tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc                                                    >  1:2778517/1‑69 (MQ=255)
tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc                                                    >  1:2748287/1‑69 (MQ=255)
tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc                                                    <  1:2648700/69‑1 (MQ=255)
tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc                                                    >  1:2540920/1‑69 (MQ=255)
tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc                                                    >  1:1913020/1‑69 (MQ=255)
tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc                                                    >  1:1652910/1‑69 (MQ=255)
tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc                                                    >  1:1009705/1‑69 (MQ=255)
tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc                                                    >  1:1288943/1‑69 (MQ=255)
tGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTc                                                    >  1:1484668/1‑69 (MQ=255)
               gttgttATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTgg                                     >  1:2721789/1‑69 (MQ=255)
               gttgttATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTgg                                     >  1:782895/1‑69 (MQ=255)
               gttgttATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTgg                                     >  1:2644615/1‑69 (MQ=255)
                      tGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTc                                >  1:566454/1‑67 (MQ=255)
                      tGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCgg                              >  1:1511350/1‑69 (MQ=255)
                      tGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCgg                              >  1:1386871/1‑69 (MQ=255)
                       gACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCgg                              <  1:1778081/68‑1 (MQ=255)
                       gACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCgg                              <  1:1613414/68‑1 (MQ=255)
                        aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGt                                                     >  1:2549692/1‑44 (MQ=255)
                        aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGc                             >  1:811681/1‑68 (MQ=255)
                        aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGc                             >  1:1443125/1‑68 (MQ=255)
                        aCCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGc                             >  1:234379/1‑68 (MQ=255)
                             gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc                       <  1:3182408/69‑1 (MQ=255)
                             gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc                       <  1:1825213/69‑1 (MQ=255)
                             gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc                       <  1:3057055/69‑1 (MQ=255)
                             gtTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc                       <  1:1454588/69‑1 (MQ=255)
                             gtTGCGGAAGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATc                       <  1:2990846/69‑1 (MQ=255)
                                         aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg            >  1:148853/1‑68 (MQ=255)
                                         aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg            >  1:817901/1‑68 (MQ=255)
                                         aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg            >  1:937652/1‑68 (MQ=255)
                                         aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg            >  1:1787895/1‑68 (MQ=255)
                                         aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg            >  1:372127/1‑68 (MQ=255)
                                         aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg            >  1:3034058/1‑68 (MQ=255)
                                         aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg            >  1:1918571/1‑68 (MQ=255)
                                         aaaaaGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCg            >  1:3142960/1‑68 (MQ=255)
                                             aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa       >  1:2944792/1‑69 (MQ=255)
                                             aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa       >  1:638323/1‑69 (MQ=255)
                                             aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa       >  1:1654183/1‑69 (MQ=255)
                                             aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa       >  1:1677454/1‑69 (MQ=255)
                                             aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa       >  1:2563538/1‑69 (MQ=255)
                                             aGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa       >  1:2137095/1‑69 (MQ=255)
                                               ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa     <  1:1201173/69‑1 (MQ=255)
                                               ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa     <  1:660962/69‑1 (MQ=255)
                                               ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa     <  1:2517649/69‑1 (MQ=255)
                                               ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa     <  1:2293559/69‑1 (MQ=255)
                                               ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa     <  1:3245806/69‑1 (MQ=255)
                                               ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa     <  1:2056622/69‑1 (MQ=255)
                                               ccGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa     <  1:2052734/69‑1 (MQ=255)
                                                  aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGCTTTCGGACCTATCGACGAGCTGCGTAAAAAACa   >  1:2233302/1‑68 (MQ=255)
                                                  aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGAt                                   >  1:2156433/1‑36 (MQ=19)
                                                  aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac  >  1:1741115/1‑68 (MQ=255)
                                                  aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac  >  1:976327/1‑68 (MQ=255)
                                                  aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac  >  1:3232381/1‑68 (MQ=255)
                                                  aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac  >  1:3133524/1‑68 (MQ=255)
                                                  aTCTGCGTCGAATCTTGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTAAAAAACac  >  1:1206470/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)
                                                    cTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa       <  1:2754409/62‑1 (MQ=255)
                                                    cTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaa       <  1:1151059/62‑1 (MQ=255)
                                                         ttGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTaaaaaa     <  1:3000979/59‑1 (MQ=255)
                                                                ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa     <  1:1202253/52‑1 (MQ=255)
                                                                ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa     <  1:2793834/52‑1 (MQ=255)
                                                                ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa     <  1:1324400/52‑1 (MQ=255)
                                                                ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa     <  1:1734795/52‑1 (MQ=255)
                                                                ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa     <  1:2467035/52‑1 (MQ=255)
                                                                ttGCCCGCTACGGGCGCTGGATTTCGGACCTATCGACGAGCTGCGTaaaaaa     <  1:2096578/52‑1 (MQ=255)
                                                                   |                                                   
TGACCAAATGCGATGGTTGTTATGACCGTGTTGCGGAGGGTAAAAAGCCGATCTGTGTTGAATCTTGTCCGCTGCGGGCGCTGGATTTCGGCCCTATCGACGAGCTGCGTAAAAAACAT  >  minE/600480‑600598

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