Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,590,221 C→T 100% V25V (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,2210CT87.0% 198.8 / 18.1 77V25V (GTG→GTApdxKpyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base C (10/0);  new base T (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.93e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAAACAACAACAA  >  minE/1590182‑1590286
                                       |                                                                 
aGCGGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:2211755/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:2160310/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:217099/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:314216/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:1985239/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1977022/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:612362/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:985483/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:928313/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:835520/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:762893/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:743756/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:73864/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:624866/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:620822/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:2080114/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:581657/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:504729/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:47983/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:435643/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:40653/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:361056/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:333990/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1046851/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1315904/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1440229/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1281460/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:148463/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1228568/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:1736045/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1744474/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1126662/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1102636/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1795392/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1080144/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1797740/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1818790/67‑1 (MQ=255)
aGCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1061568/67‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTTCACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1960028/66‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:2140351/66‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:520302/66‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1902348/66‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:316397/66‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:293570/66‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1640522/66‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1157490/66‑1 (MQ=255)
 gCAGGCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:2190066/66‑1 (MQ=255)
     gCACCGCAATGCTGTTGCCTACGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:132576/62‑1 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCTCCTGCAGTGCCCTAct                       >  1:175907/1‑66 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct                       >  1:911470/1‑66 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct                       >  1:968662/1‑66 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct                       >  1:748715/1‑66 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct                       >  1:1536104/1‑66 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct                       >  1:512356/1‑66 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct                       >  1:1484194/1‑66 (MQ=255)
                  gTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTAct                       >  1:1922995/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)
                       ctaCGCTGCCGTACACTACCTGCGACTGCACGGCAACGATATcc                                        <  1:1874795/42‑1 (MQ=255)
                                      ccaccTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAAAcaacaacaa  <  1:2110933/67‑1 (MQ=255)
                                      ccaccTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAAAcaacaacaa  <  1:195172/67‑1 (MQ=255)
                                      ccaccTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAAAcaacaacaa  <  1:105971/67‑1 (MQ=255)
                                      ccaccTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAAAcaacaacaa  <  1:919813/67‑1 (MQ=255)
                                      ccaccTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAAAcaacaacaa  <  1:826513/67‑1 (MQ=255)
                                      ccaccTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAAAcaacaacaa  <  1:1297955/67‑1 (MQ=255)
                                      ccaccTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAAAcaacaacaa  <  1:223767/67‑1 (MQ=255)
                                      ccaccTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAAAcaacaacaa  <  1:1177769/67‑1 (MQ=255)
                                      ccaccTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAAAcaacaacaa  <  1:553532/67‑1 (MQ=255)
                                      ccaccTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAAAcaacaacaa  <  1:407874/67‑1 (MQ=255)
                                       caccTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAAAcaacaacaa  <  1:1803649/66‑1 (MQ=255)
                                       |                                                                 
AGCAGGCACGGCAATGCTGTTGCCCACGCTGCCGTAAACCACCTGCGACTGCACGGCGACGATATCCGCCTGCAGTGCCCTACTCTTATCGTTAAACAACAACAA  >  minE/1590182‑1590286

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: