Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,590,191 G→C 100% A35A (GCC→GCG pdxK ← pyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,590,1910GC92.8% 220.6 / 7.2 69A35A (GCC→GCGpdxKpyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base C (0/64);  total (2/67)
Fisher's exact test for biased strand distribution p-value = 4.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCC  >  minE/1590133‑1590248
                                                          |                                                         
tCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCt                                                   >  1:1963175/1‑67 (MQ=255)
tCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCt                                                   >  1:629370/1‑67 (MQ=255)
         ccGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTg                                               <  1:141968/62‑1 (MQ=255)
         ccGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTg                                               <  1:1363648/62‑1 (MQ=255)
                cacggcaAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCg                                   <  1:4605/67‑1 (MQ=255)
                                                 aGCGGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:2211755/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:2210821/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:333990/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:314216/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:308042/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:296833/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:225517/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:435643/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:2187410/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:2178575/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:217099/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:2160310/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:2080114/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1034366/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:361056/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:40653/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:2042338/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:47983/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:504729/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:581657/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:612362/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:620822/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:624866/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:73864/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:743756/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:762893/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:835520/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:928313/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:985483/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1685572/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1046851/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1061568/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1080144/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1102636/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1126662/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1228568/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1281460/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1315904/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1440229/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:148463/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1556276/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1606390/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:166879/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:2042272/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1716855/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1736045/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:2041878/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1995600/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1985239/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1977022/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1818790/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1797740/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1795392/67‑1 (MQ=255)
                                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1744474/67‑1 (MQ=255)
                                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTTCACTACCTGCGACTGCACGGCAACGATATcc  <  1:1960028/66‑1 (MQ=255)
                                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:520302/66‑1 (MQ=255)
                                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:293570/66‑1 (MQ=255)
                                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1902348/66‑1 (MQ=255)
                                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:2190066/66‑1 (MQ=255)
                                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1157490/66‑1 (MQ=255)
                                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:1640522/66‑1 (MQ=255)
                                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:2140351/66‑1 (MQ=255)
                                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:316397/66‑1 (MQ=255)
                                                      gCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc  <  1:132576/62‑1 (MQ=255)
                                                          |                                                         
TCAGCAATACCGTCGGCACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCC  >  minE/1590133‑1590248

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: