Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 600,667 C→T 53.5% P185P (CCC→CCT dmsB → dimethyl sulfoxide reductase, anaerobic, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE600,6670CT53.5% ‑1.7 / 34.7 43P185P (CCC→CCTdmsBdimethyl sulfoxide reductase, anaerobic, subunit B
Reads supporting (aligned to +/- strand):  ref base C (2/18);  new base T (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 = 9.80e-01

GGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAG  >  minE/600599‑600731
                                                                    |                                                                
ggCGATCTGGTGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAAccc                                                                  <  1:1022455/69‑1 (MQ=255)
ggCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAAccc                                                                  <  1:518891/69‑1 (MQ=255)
ggCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAAccc                                                                  <  1:1361765/69‑1 (MQ=255)
ggCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAAccc                                                                  <  1:2581829/69‑1 (MQ=255)
ggCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAAccc                                                                  <  1:2261339/69‑1 (MQ=255)
ggCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAAccc                                                                  <  1:1615599/69‑1 (MQ=255)
ggCGATATGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAAccc                                                                  <  1:483429/69‑1 (MQ=255)
           ggcCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGcc                                                            <  1:946247/64‑1 (MQ=255)
           ggcCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGcc                                                            <  1:579287/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:2068549/64‑1 (MQ=255)
                      cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:1694529/69‑1 (MQ=14)
                      cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:1986411/69‑1 (MQ=14)
                      cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:611161/69‑1 (MQ=14)
                      cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:1998667/69‑1 (MQ=14)
                      cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:1132918/69‑1 (MQ=14)
                      cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:1612339/69‑1 (MQ=14)
                      cgttgcCGCGTGCGCACTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCg                                            <  1:2754564/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:1441293/1‑69 (MQ=17)
                                     aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                             >  1:628226/1‑69 (MQ=17)
                                     aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                             >  1:1998413/1‑69 (MQ=17)
                                     aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                             >  1:2634767/1‑69 (MQ=17)
                                     aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                             >  1:2741310/1‑69 (MQ=17)
                                     aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                             >  1:125046/1‑69 (MQ=17)
                                     aCTTCACCAAACCGAATATTGTGATCAAACCTAACGCCAATAGCCGCCCAACCGGGGATACCACTGGCt                             >  1:1235497/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:1363304/1‑69 (MQ=255)
                                                ccGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAAc                  >  1:2320700/1‑69 (MQ=255)
                                                      atTGTTATCAAACCCAATGCCAATAGCCGCCCGACCggg                                          <  1:361086/39‑1 (MQ=255)
                                                      atTGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                          <  1:2140658/39‑1 (MQ=255)
                                                      atTGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                          <  1:70222/39‑1 (MQ=255)
                                                      atTGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                          <  1:1015799/39‑1 (MQ=255)
                                                       ttGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                          <  1:1249028/38‑1 (MQ=255)
                                                       ttGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                          <  1:541133/38‑1 (MQ=255)
                                                       ttGTGATCAAACCCAATGCCAATAGCCGCCCGACCggg                                          <  1:854502/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)
                                                                    |                                                                
GGCGATCTGGCGGCCGTTGCGCCGTTGCCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAG  >  minE/600599‑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: