Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,590,206 C→T 100% V30V (GTG→GTA pdxK ← pyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,590,2060CT85.7% 191.7 / 20.7 77V30V (GTG→GTApdxKpyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base C (10/1);  new base T (0/66);  total (10/67)
Fisher's exact test for biased strand distribution p-value = 1.00e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACT  >  minE/1590149‑1590265
                                                         |                                                           
cacggcaAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCg                                                    <  1:4605/67‑1 (MQ=255)
                                 aGCGGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:2211755/67‑1 (MQ=255)
                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:2178575/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:361056/67‑1 (MQ=255)
                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:217099/67‑1 (MQ=255)
                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:2160310/67‑1 (MQ=255)
                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1034366/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:333990/67‑1 (MQ=255)
                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:2080114/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:1977022/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:1685572/67‑1 (MQ=255)
                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1716855/67‑1 (MQ=255)
                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1985239/67‑1 (MQ=255)
                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1744474/67‑1 (MQ=255)
                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1795392/67‑1 (MQ=255)
                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1797740/67‑1 (MQ=255)
                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1736045/67‑1 (MQ=255)
                                 aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1818790/67‑1 (MQ=255)
                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTTCACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1960028/66‑1 (MQ=255)
                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:520302/66‑1 (MQ=255)
                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1902348/66‑1 (MQ=255)
                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1157490/66‑1 (MQ=255)
                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:2140351/66‑1 (MQ=255)
                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:316397/66‑1 (MQ=255)
                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:293570/66‑1 (MQ=255)
                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:2190066/66‑1 (MQ=255)
                                  gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:1640522/66‑1 (MQ=255)
                                      gCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                   <  1:132576/62‑1 (MQ=255)
                                                   gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCTCCTGCAGTGCCCTAct  >  1:175907/1‑66 (MQ=255)
                                                   gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct  >  1:1536104/1‑66 (MQ=255)
                                                   gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct  >  1:1484194/1‑66 (MQ=255)
                                                   gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct  >  1:512356/1‑66 (MQ=255)
                                                   gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct  >  1:748715/1‑66 (MQ=255)
                                                   gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct  >  1:1922995/1‑66 (MQ=255)
                                                   gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct  >  1:911470/1‑66 (MQ=255)
                                                   gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct  >  1:968662/1‑66 (MQ=255)
                                                   gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct  >  1:1822590/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)
                                                         |                                                           
CACGGCAAAGACATTCAGGCCGTTCTGTTTGATAGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACT  >  minE/1590149‑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: