Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,744,420 Δ1 bp 100% coding (1315/1623 nt) nadB → quinolinate synthase, L‑aspartate oxidase (B protein) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,744,4200A.96.0% 101.4 / ‑2.1 25coding (1315/1623 nt)nadBquinolinate synthase, L‑aspartate oxidase (B protein) subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (0/24);  total (0/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAA  >  minE/1744374‑1744478
                                              |                                                          
tGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTg                                        <  1:1363033/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:995892/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:1223313/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:948395/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:834417/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:767587/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:649057/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:613595/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:609250/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:585930/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:572765/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:487253/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:275356/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:2126640/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:1887239/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:1787372/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:176986/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:1738208/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:1688921/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:1679903/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:1675973/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:1400153/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATAACGTTGGCATTGTGCGCACAACGaa  <  1:2068498/67‑1 (MQ=255)
                                     aTTCAGCATACTGGCACGAGCTACGTCTGTTTATCTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:817480/67‑1 (MQ=255)
                                      ttCAGCATACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGaa  <  1:1536268/66‑1 (MQ=255)
                                              |                                                          
TGAAAGCCGCGTTGAGAACCCTGACGAACGGGTAGTAATTCAGCATAACTGGCACGAGCTACGTCTGTTTATGTGGGATTACGTTGGCATTGTGCGCACAACGAA  >  minE/1744374‑1744478

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