Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 598,245 C→T 39.9% G196G (GGC→GGT dmsA → dimethyl sulfoxide reductase, anaerobic, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE598,2450CT39.9% 16.8 / 23.0 40G196G (GGC→GGTdmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base C (12/12);  new base T (8/8);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTA  >  minE/598190‑598301
                                                       |                                                        
gCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCAtt                                             >  1:858150/1‑69 (MQ=255)
gCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCAtt                                             >  1:128063/1‑69 (MQ=255)
gCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCAtt                                             >  1:51933/1‑69 (MQ=255)
gCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCAtt                                             >  1:2357733/1‑69 (MQ=255)
gCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCAtt                                             >  1:1872194/1‑69 (MQ=255)
                     cccTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCcgcg                         >  1:224691/1‑68 (MQ=255)
                     cccTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCcgcg                         >  1:1849776/1‑68 (MQ=255)
                     cccTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCcgc                          >  1:2168812/1‑67 (MQ=255)
                        tGGTCGCGCGGCTAATGAACTGCTGCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCaaa                     >  1:976480/1‑69 (MQ=255)
                        tGGTCGCGCGGCTAATGAACTGCTGCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCaaa                     >  1:70364/1‑69 (MQ=255)
                        tGGTCGCGCGGCTAATGAACTGCTGCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCaaa                     >  1:2833711/1‑69 (MQ=255)
                        tGGTCGCGCGGCTAATGAACTGCTGCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCaaa                     >  1:2573647/1‑69 (MQ=255)
                        tGGTCGCGCGGCTAATGAACTGCTGCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCaaa                     >  1:1743985/1‑69 (MQ=255)
                        tGGTCGCGCGGCTAATGAACTGCTGCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCaaa                     >  1:1425882/1‑69 (MQ=255)
                        tGGTCGCGCGGCTAATGAACTGCTGCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCaaa                     >  1:151911/1‑69 (MQ=255)
                        tGGTCGCGCGGCTAATGAACAGCTGCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCaaa                     >  1:2870499/1‑69 (MQ=255)
                          gTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGGACCATTACGGCGACTACTCCTCCGCGCAAAtt                   >  1:1643637/1‑69 (MQ=255)
                          gTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAAtt                   >  1:1678160/1‑69 (MQ=255)
                          gTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAAtt                   >  1:2261557/1‑69 (MQ=255)
                          gTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAAACATTACGGCGACTACTCCTCCGCGCAAAtt                   >  1:1765870/1‑69 (MQ=255)
                                        gaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTg       <  1:1927173/67‑1 (MQ=255)
                                        gaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTg       <  1:788063/67‑1 (MQ=255)
                                        gaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTg       <  1:536917/67‑1 (MQ=255)
                                        gaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTg       <  1:1475076/67‑1 (MQ=255)
                                        gaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTg       <  1:2653944/67‑1 (MQ=255)
                                        gaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTg       <  1:1725796/67‑1 (MQ=255)
                                        gaACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTg       <  1:216511/67‑1 (MQ=255)
                                         aaCTGCTGCGGCGGCTATCTGCACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAAc    <  1:1035475/69‑1 (MQ=255)
                                         aaCTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAAc    <  1:2520035/69‑1 (MQ=255)
                                         aaCTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAAc    <  1:153515/69‑1 (MQ=255)
                                         aaCTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAAc    <  1:3130460/69‑1 (MQ=255)
                                         aaCTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAAc    <  1:611823/69‑1 (MQ=255)
                                           ctgctgCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTCTGAACta  <  1:2499329/69‑1 (MQ=255)
                                           ctgctgCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTCTGAACta  <  1:2332049/69‑1 (MQ=255)
                                           ctgctgCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTCTGAACta  <  1:1457156/69‑1 (MQ=255)
                                           ctgctgCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTCTGAACta  <  1:626833/69‑1 (MQ=255)
                                           ctgctgCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTCTGAACta  <  1:1885844/69‑1 (MQ=255)
                                           ctgctgCGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTCTGAACta  <  1:909878/69‑1 (MQ=255)
                                            tgctgcGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTCTGAACta  <  1:2123543/68‑1 (MQ=255)
                                            tgctgcGGCGGTTATCTGAACCATTATGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTCTGAACta  <  1:1050960/68‑1 (MQ=255)
                                                       |                                                        
GCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTA  >  minE/598190‑598301

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