Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 600,640 T→C 40.7% F176F (TTT→TTC dmsB → dimethyl sulfoxide reductase, anaerobic, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE600,6400TC40.7% 28.7 / 37.0 54F176F (TTT→TTCdmsBdimethyl sulfoxide reductase, anaerobic, subunit B
Reads supporting (aligned to +/- strand):  ref base T (11/21);  new base C (12/10);  total (23/31)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-02

CTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCT  >  minE/600573‑600704
                                                                   |                                                                
cTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGTGCTCACTTTa                                                                 <  1:2569101/69‑1 (MQ=255)
cTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTa                                                                 <  1:751555/69‑1 (MQ=255)
cTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTa                                                                 <  1:2316222/69‑1 (MQ=255)
cTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTa                                                                 <  1:2313159/69‑1 (MQ=255)
cTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTa                                                                 <  1:1805159/69‑1 (MQ=255)
cTATCGACGAGCTGCGCAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTa                                                                 <  1:3235764/69‑1 (MQ=255)
   tCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCa                                                              >  1:3254663/1‑69 (MQ=255)
   tCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCa                                                              >  1:2789216/1‑69 (MQ=255)
   tCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCa                                                              >  1:2545236/1‑69 (MQ=255)
   tCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCa                                                              >  1:2388409/1‑69 (MQ=255)
   tCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCa                                                              >  1:2174998/1‑69 (MQ=255)
   tCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCa                                                              >  1:3165755/1‑69 (MQ=255)
   tCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCa                                                              >  1:329078/1‑69 (MQ=255)
   tCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCa                                                              >  1:1789396/1‑69 (MQ=255)
   tCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCa                                                              >  1:3338724/1‑69 (MQ=255)
   tCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCa                                                              >  1:838285/1‑69 (MQ=255)
   tCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCa                                                              >  1:539101/1‑69 (MQ=255)
          gCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGaa                                                       <  1:1642728/69‑1 (MQ=255)
          gCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGaa                                                       <  1:2944504/69‑1 (MQ=255)
                          ggCGATCTGGTGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAAccc                                       <  1:1022455/69‑1 (MQ=255)
                          ggCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAAccc                                       <  1:2261339/69‑1 (MQ=255)
                          ggCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAAccc                                       <  1:518891/69‑1 (MQ=255)
                          ggCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAAccc                                       <  1:1615599/69‑1 (MQ=255)
                          ggCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAAccc                                       <  1:1361765/69‑1 (MQ=255)
                          ggCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAAccc                                       <  1:2581829/69‑1 (MQ=255)
                          ggCGATATGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAAccc                                       <  1:483429/69‑1 (MQ=255)
                                     ggcCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGcc                                 <  1:2946783/64‑1 (MQ=255)
                                     ggcCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGcc                                 <  1:579287/64‑1 (MQ=255)
                                     ggcCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGcc                                 <  1:3090806/64‑1 (MQ=255)
                                     ggcCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGcc                                 <  1:946247/64‑1 (MQ=255)
                                     ggcCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGcc                                 <  1:2565478/64‑1 (MQ=255)
                                     ggcCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGcc                                 <  1:2068549/64‑1 (MQ=255)
                                                cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                 <  1:1132918/69‑1 (MQ=14)
                                                cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                 <  1:2754564/69‑1 (MQ=14)
                                                cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                 <  1:611161/69‑1 (MQ=14)
                                                cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                 <  1:2294519/69‑1 (MQ=14)
                                                cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                 <  1:1612339/69‑1 (MQ=14)
                                                cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                 <  1:1694529/69‑1 (MQ=14)
                                                cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                 <  1:1998667/69‑1 (MQ=14)
                                                cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                 <  1:1986411/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:2614088/1‑68 (MQ=17)
                                                               aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGc   >  1:926649/1‑68 (MQ=17)
                                                               aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt  >  1:1998413/1‑69 (MQ=17)
                                                               aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt  >  1:2634767/1‑69 (MQ=17)
                                                               aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt  >  1:405885/1‑69 (MQ=17)
                                                               aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt  >  1:1441293/1‑69 (MQ=17)
                                                               aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt  >  1:628226/1‑69 (MQ=17)
                                                               aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt  >  1:125046/1‑69 (MQ=17)
                                                               aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt  >  1:1235497/1‑69 (MQ=17)
                                                               aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt  >  1:2741310/1‑69 (MQ=17)
                                                               aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCAGCCCAACCGGGGATACCACTGGCt  >  1:2951432/1‑69 (MQ=17)
                                                                   |                                                                
CTATCGACGAGCTGCGTAAAAAACATGGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCT  >  minE/600573‑600704

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: