Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,790,195 T→A 100% L14* (TTG→TAG)  gabD → succinate‑semialdehyde dehydrogenase I, NADP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,790,1950TA84.1% 149.4 / 23.3 63L14* (TTG→TAG) gabDsuccinate‑semialdehyde dehydrogenase I, NADP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/10);  new base A (53/0);  total (53/10)
Fisher's exact test for biased strand distribution p-value = 7.82e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.44e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCA  >  minE/1790154‑1790245
                                         |                                                  
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGAtt                                                >  1:1878891/1‑46 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:475061/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:1964218/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:1066270/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:2117908/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:2124694/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:2150859/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:250171/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:306846/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:322920/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:330661/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:368697/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:37462/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:389039/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:400753/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:20562/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:48804/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:526359/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:535866/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:621825/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:629087/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:638890/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:664818/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:695279/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:753641/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:814882/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:899652/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:977310/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:986340/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:1071648/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:1254675/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:1273603/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:1192668/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:1292837/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:1180083/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:1383812/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:193418/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:1418646/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:148199/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:1497324/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:1512397/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:1572192/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:1572882/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:1610205/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:183015/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:1893663/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGcc                           >  1:1919062/1‑67 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGc                            >  1:1449509/1‑66 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGc                            >  1:827147/1‑66 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGc                            >  1:1341174/1‑66 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGc                            >  1:1301527/1‑66 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGc                            >  1:1072128/1‑66 (MQ=255)
gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTAGATTAACGGGGAATGGCTGGACGc                            >  1:751906/1‑66 (MQ=255)
                          tCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCa  <  1:1155062/66‑1 (MQ=255)
                          tCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCa  <  1:924463/66‑1 (MQ=255)
                          tCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCa  <  1:1418228/66‑1 (MQ=255)
                          tCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCa  <  1:647824/66‑1 (MQ=255)
                          tCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCa  <  1:1369466/66‑1 (MQ=255)
                          tCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCa  <  1:483222/66‑1 (MQ=255)
                          tCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCa  <  1:1880676/66‑1 (MQ=255)
                          tCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCa  <  1:1926453/66‑1 (MQ=255)
                           ccGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCa  <  1:180875/65‑1 (MQ=255)
                           ccGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCa  <  1:214064/65‑1 (MQ=255)
                                         |                                                  
GATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCA  >  minE/1790154‑1790245

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: