Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,590,218 A→C 100% V26V (GTT→GTG pdxK ← pyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,590,2180AC87.0% 235.6 / 16.6 77V26V (GTT→GTGpdxKpyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base A (10/0);  new base C (0/67);  total (10/67)
Fisher's exact test for biased strand distribution p-value = 9.12e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACT  >  minE/1590182‑1590265
                                    |                                               
aGCGGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:2211755/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:217099/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:2210821/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:2187410/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:2178575/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:333990/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:2160310/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:2080114/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:2042338/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:2042272/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:1046851/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:361056/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:40653/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:435643/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:1797740/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1606390/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:166879/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1685572/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1556276/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1716855/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1736045/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1744474/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:148463/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1440229/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1795392/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1818790/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1315904/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1281460/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1228568/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1126662/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1102636/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1080144/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1061568/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1034366/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1977022/67‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTTCACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1960028/66‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1640522/66‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:293570/66‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:316397/66‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:2190066/66‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1902348/66‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:520302/66‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:2140351/66‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1157490/66‑1 (MQ=255)
     gCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:132576/62‑1 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCTCCTGCAGTGCCCTAct  >  1:175907/1‑66 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct  >  1:1922995/1‑66 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct  >  1:968662/1‑66 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct  >  1:911470/1‑66 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct  >  1:748715/1‑66 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct  >  1:1822590/1‑66 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct  >  1:512356/1‑66 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct  >  1:1484194/1‑66 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct  >  1:1536104/1‑66 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGAGCactact  >  1:213137/1‑66 (MQ=255)
                     gCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:2125895/46‑1 (MQ=255)
                     gCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:2054985/46‑1 (MQ=255)
                       ctaCGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1874795/42‑1 (MQ=255)
                                    |                                               
AGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACT  >  minE/1590182‑1590265

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