Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 600,670 T→C 53.5% N186N (AAT→AAC dmsB → dimethyl sulfoxide reductase, anaerobic, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE600,6700TC53.5% ‑3.5 / 33.0 43N186N (AAT→AACdmsBdimethyl sulfoxide reductase, anaerobic, subunit B
Reads supporting (aligned to +/- strand):  ref base T (2/18);  new base C (12/11);  total (14/29)
Fisher's exact test for biased strand distribution p-value = 4.01e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.78e-02

GGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAG  >  minE/600610‑600731
                                                            |                                                             
ggcCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGcc                                                            <  1:946247/64‑1 (MQ=255)
ggcCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGcc                                                            <  1:3090806/64‑1 (MQ=255)
ggcCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGcc                                                            <  1:2946783/64‑1 (MQ=255)
ggcCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGcc                                                            <  1:2565478/64‑1 (MQ=255)
ggcCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGcc                                                            <  1:579287/64‑1 (MQ=255)
ggcCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGcc                                                            <  1:2068549/64‑1 (MQ=255)
           cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:1132918/69‑1 (MQ=14)
           cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:1612339/69‑1 (MQ=14)
           cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:1694529/69‑1 (MQ=14)
           cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:1986411/69‑1 (MQ=14)
           cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:2754564/69‑1 (MQ=14)
           cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:1998667/69‑1 (MQ=14)
           cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:611161/69‑1 (MQ=14)
           cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:2294519/69‑1 (MQ=14)
            gttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:1869732/68‑1 (MQ=14)
            gttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:2023196/68‑1 (MQ=14)
                          ccTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGc                              >  1:1415494/2‑68 (MQ=17)
                          aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGc                              >  1:926649/1‑68 (MQ=17)
                          aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGc                              >  1:2614088/1‑68 (MQ=17)
                          aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                             >  1:405885/1‑69 (MQ=17)
                          aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                             >  1:628226/1‑69 (MQ=17)
                          aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                             >  1:1235497/1‑69 (MQ=17)
                          aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                             >  1:125046/1‑69 (MQ=17)
                          aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                             >  1:2634767/1‑69 (MQ=17)
                          aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                             >  1:2741310/1‑69 (MQ=17)
                          aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                             >  1:1998413/1‑69 (MQ=17)
                          aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                             >  1:1441293/1‑69 (MQ=17)
                          aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCAGCCCAACCGGGGATACCACTGGCt                             >  1:2951432/1‑69 (MQ=17)
                                     ccGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:3091324/43‑1 (MQ=17)
                                     ccGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAAc                  >  1:2320700/1‑69 (MQ=255)
                                     ccGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAAc                  >  1:1363304/1‑69 (MQ=255)
                                           atTGTTATCAAACCCAATGCCAATAGCCGCCCGACCggg                                          <  1:361086/39‑1 (MQ=255)
                                           atTGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                          <  1:1015799/39‑1 (MQ=255)
                                           atTGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                          <  1:2140658/39‑1 (MQ=255)
                                           atTGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                          <  1:70222/39‑1 (MQ=255)
                                            ttGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                          <  1:541133/38‑1 (MQ=255)
                                            ttGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                          <  1:854502/38‑1 (MQ=255)
                                            ttGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                          <  1:1249028/38‑1 (MQ=255)
                                            ttGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                          <  1:933354/38‑1 (MQ=255)
                                                 aTCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGgt       <  1:2159385/68‑1 (MQ=255)
                                                  tCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGgt       <  1:3029519/67‑1 (MQ=255)
                                                    aaaCCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTga   <  1:1057555/69‑1 (MQ=255)
                                                      aCCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTgag  <  1:2222015/68‑1 (MQ=255)
                                                            |                                                             
GGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAG  >  minE/600610‑600731

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