Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 600,685 G→A 54.8% P191P (CCG→CCA dmsB → dimethyl sulfoxide reductase, anaerobic, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE600,6850GA54.8% ‑3.9 / 32.5 42P191P (CCG→CCAdmsBdimethyl sulfoxide reductase, anaerobic, subunit B
Reads supporting (aligned to +/- strand):  ref base G (6/13);  new base A (12/11);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATG  >  minE/600621‑600750
                                                                |                                                                 
cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                                               <  1:1694529/69‑1 (MQ=14)
cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                                               <  1:611161/69‑1 (MQ=14)
cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                                               <  1:2294519/69‑1 (MQ=14)
cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                                               <  1:1998667/69‑1 (MQ=14)
cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                                               <  1:1986411/69‑1 (MQ=14)
cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                                               <  1:2754564/69‑1 (MQ=14)
cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                                               <  1:1612339/69‑1 (MQ=14)
cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                                               <  1:1132918/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:1441293/1‑69 (MQ=17)
               aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                                                >  1:125046/1‑69 (MQ=17)
               aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                                                >  1:1998413/1‑69 (MQ=17)
               aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                                                >  1:1235497/1‑69 (MQ=17)
               aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                                                >  1:2741310/1‑69 (MQ=17)
               aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                                                >  1:2634767/1‑69 (MQ=17)
               aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                                                >  1:628226/1‑69 (MQ=17)
               aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                                                >  1:405885/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:70222/39‑1 (MQ=255)
                                atTGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                                             <  1:2140658/39‑1 (MQ=255)
                                atTGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                                             <  1:1015799/39‑1 (MQ=255)
                                 ttGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                                             <  1:854502/38‑1 (MQ=255)
                                 ttGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                                             <  1:541133/38‑1 (MQ=255)
                                 ttGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                                             <  1:933354/38‑1 (MQ=255)
                                 ttGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                                             <  1:1249028/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)
                                                           cgccCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGc     >  1:3085017/1‑68 (MQ=255)
                                                           cgccCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGc     >  1:2511230/1‑68 (MQ=255)
                                                           cgccCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGc     >  1:680719/1‑68 (MQ=255)
                                                           cgccCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGc     >  1:1646155/1‑68 (MQ=255)
                                                             cccGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATg  <  1:3027727/69‑1 (MQ=255)
                                                                gACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGgt                          <  1:2983069/42‑1 (MQ=255)
                                                                |                                                                 
CGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATG  >  minE/600621‑600750

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: