Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 601,323 G→A 100% V198M (GTG→ATG)  dmsC → dimethyl sulfoxide reductase, anaerobic, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE601,3230GA82.0% 46.7 / 3.7 33V198M (GTG→ATG) dmsCdimethyl sulfoxide reductase, anaerobic, subunit C
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base A (0/27);  total (1/32)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCT  >  minE/601260‑601387
                                                               |                                                                
gATGTCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGAGTCGTGATGTCa                                                             <  1:2700124/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:1834069/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:2695973/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:2663611/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:2767115/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:2437765/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:2424441/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:2296463/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:2248369/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:2242669/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:2894875/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:2174528/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:2110302/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:2915357/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:2954829/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:1669738/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:3300349/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:1533111/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:465020/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:1285133/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:1200161/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:1171902/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:702713/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTCATGTCa                                                             <  1:441740/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCCATTTCTGTTCTGGCACTGGTAGTGAGGGGCGTGATGTCa                                                             <  1:2915567/69‑1 (MQ=255)
gATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCa                                                             >  1:1298157/1‑69 (MQ=255)
                   tgctgcCAGCCATTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:3005822/50‑1 (MQ=255)
                     ctgcCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCt                                         <  1:2195345/68‑1 (MQ=255)
                     ctgcCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCt                                         <  1:576384/68‑1 (MQ=255)
                       gccagccaTTTCTGTTCTGGCACTGGTAGTGAGTGGCGTGATGTCa                                                             <  1:1118055/46‑1 (MQ=255)
                       gcCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCt                                         <  1:2528714/66‑1 (MQ=255)
                       gcCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCt                                         <  1:1888017/66‑1 (MQ=255)
                       gcCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCt                                         <  1:297145/66‑1 (MQ=255)
                                                            gtggtgTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCt  <  1:2195230/68‑1 (MQ=255)
                                                            gtggtgTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCt  <  1:1792384/68‑1 (MQ=255)
                                                            gtggtgTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCt  <  1:2945402/68‑1 (MQ=255)
                                                            gtggtgTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCt  <  1:1749194/68‑1 (MQ=255)
                                                            gtggtgTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCt  <  1:1736818/68‑1 (MQ=255)
                                                            gtggtgTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCt  <  1:1622306/68‑1 (MQ=255)
                                                            gtggtgTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCt  <  1:1490974/68‑1 (MQ=255)
                                                            gtggtgTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCt  <  1:1441404/68‑1 (MQ=255)
                                                            gtggtgTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCt  <  1:1316421/68‑1 (MQ=255)
                                                            gtggtgTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCt  <  1:707509/68‑1 (MQ=255)
                                                            gtggtgTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCt  <  1:737421/68‑1 (MQ=255)
                                                            gtggtgTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCt  <  1:832513/68‑1 (MQ=255)
                                                             tggtgTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCt  <  1:1114290/67‑1 (MQ=255)
                                                               |                                                                
GATGGCTGGGCGATGCGTCTGCTGCCAGCAATTTCTGTACTGGCACTGGTAGTGAGTGGCGTGGTGTCAGTGATGCAGGGCGCAGAGCTGGCGACTATTCATAGTTCTGTGCAGCAGGCCGCAGCGCT  >  minE/601260‑601387

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