Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 910,914 | 0 | A | G | 43.9% | 77.8 / 11.1 | 57 | S46P (TCC→CCC) | poxB | pyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding |
Reads supporting (aligned to +/- strand): ref base A (6/26); new base G (25/0); total (31/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.52e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTC > NC_000913/910768‑911062 | cGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATc < 1:327821/151‑1 (MQ=255) aCGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATc > 1:805251/1‑147 (MQ=255) tGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATgg < 1:765667/150‑1 (MQ=255) gCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTg < 2:337341/151‑1 (MQ=255) cAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATCCAACCGGTGGGGcccccgc > 2:931853/1‑147 (MQ=255) aaTCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑g < 1:885325/151‑1 (MQ=255) aaTCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGGCTCCACTcccgtctgcccccg > 2:978157/1‑136 (MQ=255) aTCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑Ga < 2:416060/151‑1 (MQ=255) aTCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCCTCCACCCGATGGTGgcccttc > 2:1110663/1‑145 (MQ=255) aaCAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GAt < 2:579209/148‑1 (MQ=255) cAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGGGATCCACCCGCAGGTGGCCacacg > 1:967330/1‑140 (MQ=255) cGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATCCACCCGATGGTGCCCATGC‑GATTTAGACTGt > 1:335103/1‑150 (MQ=255) cGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGACATCCACTCCATGGGGCCCATGC‑GATTAAGaacg > 2:647407/1‑146 (MQ=255) gTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGGCATTCCCTCGATGGGGCCCCTGC‑GATTAAAACTGt > 1:634839/1‑149 (MQ=255) ttGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCa < 2:619859/150‑1 (MQ=255) tGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGGAATCCCCTCGATGGTGCCCCTTT‑GACTAAGACTGGCACt > 1:756253/1‑151 (MQ=255) ttAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATCCACCCGATGGGGGCCCCGC‑GGTTAAGAATGGCACTaa > 2:744410/1‑150 (MQ=255) tAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAg < 1:397038/150‑1 (MQ=255) gCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGt < 2:155289/150‑1 (MQ=255) cAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGtt < 2:574186/150‑1 (MQ=255) cAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGTTc < 1:636668/151‑1 (MQ=255) cAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATCCACCCCATGGGGCCCCTGC‑GAGTAAGAATGTCACTAAGAACGTTc > 2:398047/1‑151 (MQ=255) ggTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGTTCag < 1:759349/151‑1 (MQ=255) gCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGTTCagaga < 2:433168/150‑1 (MQ=255) ccGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATCCACTCGATGGGGCCCCTGC‑GATTAAGACTGTCACTAAGACCGTTCagaga > 2:433020/1‑149 (MQ=255) cGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATCCACCCCCTGGGGCCCATGC‑GGTTAAGACTGTCACTAAGACCGGTCagaga > 2:93045/1‑148 (MQ=255) cGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGACATCCACTCCCTGGTGGCCCAGG‑GGTTAAGAATGTTCCTAAGACCGGTc > 2:433943/1‑143 (MQ=255) ggggCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATTCACTCCATGGTGCCCCTGC‑GATTAAAAACGTCACCAAGaaccgtcaca > 1:491025/1‑136 (MQ=255) cACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCCCTCGATGGGGGCCCTGC‑GATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGcc > 1:158358/1‑146 (MQ=255) gATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGGCATCCCCCCGATGGTGCCCATGC‑GATTAAGACTGGCACTAAGACCGTTCAGAGAGGCGCCCGTGACTc > 1:640222/1‑151 (MQ=255) tCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATCCCCTCGATGGGGCCCCCGC‑GCTTTAAACTGTCACTCAGACCCGTCAGAGAGTCGCCTGTGTCTccc > 2:731777/1‑151 (MQ=255) ccGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCCTTCTCCCGGATGTGCCCATGCGGGTTAAAACGGTC‑CTAAAACCGGTCCGAGAGGCGCCTCTGGCTcccc > 2:832719/1‑151 (MQ=255) ggCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATCCCCTCGATGGGGCCCATGC‑GGTTTAAAATGTTACTAAAAACGTTTAGAGAGGCGCCCGTGACTcccc > 1:441049/1‑149 (MQ=255) ggCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGCGGGCATCCACTCGCTGGTGGCCCCGC‑GATTTAGACTGTCACTAAGACCGTTCAGAGAGTCGCCCGTGGATcccc > 2:466065/1‑149 (MQ=255) gcgcAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGGAACCCCAAGAGGGTGCCCCAGC‑GGCTAAGCACGATACCAACACCGTTCAGAGAGGTGACCGTGACcccc > 2:159032/1‑147 (MQ=255) cGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTc < 2:129293/151‑1 (MQ=255) ttGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGGATCCACTCGATGGGGCCCATGC‑GATTAAAACTGTCACTAAAAACGTTTAGAGAGTCGCCTGTGACTCCCCAGATGCGGTTTACCCCCGCCGGtttggg > 1:268101/1‑146 (MQ=255) gtgCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAgtgt < 2:335041/151‑1 (MQ=255) tgCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGGATCCACTCGATGGGGGCCATGC‑GATTAAAAATGTCACTAAAAACGTTCAGAGAGTCGCCCGTGACTCCCCCGATGCGGTTTACCCCCGCCGATTCgcggggt > 1:388964/1‑145 (MQ=255) gCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGtt < 1:471789/150‑1 (MQ=255) cTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGGCATTCACCCGATGGGGCCCCTGC‑GATTAAAACTGTCACTAAAACCGTTCAGAGAGTCGCCCGTGGCTCCCCCGAAGCGGTTCACCCCCGCCGATTCGAGTGtttt > 1:131076/1‑151 (MQ=255) ttCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGAAATCCACTCGGAGGTGCCCCTGC‑GAATAAAAATGTCACTAAGACCGTTCTGAGAATCGCCTGTGACTCCCCAGATGCGGTTTACCCCTGCCGGGTCGAGGGgttt > 1:282070/1‑150 (MQ=255) aGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAATCCCGTTCAGAGAGTCTCCTGTGACTCCCCAGATGCGTTTCCCCCCTGCCGATTCGAGTGtttt < 2:739351/147‑1 (MQ=255) aGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGAGGAATACCACTCGCTCGTGCCCAAGC‑GGTTAAACAAGTTCCTAAGACCCTTTCGAGAGGCGCCCGTGGCCCCCCAACGGCGTTTCACCCCCCGCCAGTCGAGTGTCTTgggg > 2:935176/1‑149 (MQ=255) cgcCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACAACCACTCTATTTTGCCCCTGC‑TATTTAGACTAGTACTAACACCCTTCATAGAGAAGCCCATGACTCCCCACAATCTGTTCTCCCCTTTCCATTCTAGTGtttt > 1:486756/1‑145 (MQ=255) cgcCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGGGGGCATCCACTCGATGGTGGCCCCGC‑GCTTTAAACTGTCACTAAGACCGTTCAGAGAGTGGCCCGTGGCTCCCCCGATGCGTTTTACCCCTCCCGCTTCGCGGGtttttggg > 2:54378/1‑145 (MQ=255) gcCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGata < 1:1017280/151‑1 (MQ=255) cAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTTGGCATCCCCCCCATTGTGCCCACTC‑GGTTTAAAATGTCCCTAAGACCGGTCCGAGAGGCGCCTGTTACTCCCCCGAGGCGTTTTACCCCCTCCGGTTCGAGTGgttt > 2:313462/1‑142 (MQ=255) gCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCTCTAAGACCGTTCAGAGAGTCTCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGc < 2:413724/149‑1 (MQ=255) gCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTg < 1:368171/151‑1 (MQ=255) aaaGGCCGCCACTTCTTCGTGGCGGGTGGGCATCCCCCCCGTGGTGCCCCACC‑GGTTAAGACTGTCACTAAGAACGTTCAGAGAGGCGCCCGTGACTCCCCACATGCGGTTTACCCCCGCCCGTTCGAGGGGTTTTGCGgtat > 2:615802/1‑143 (MQ=255) cGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATgg < 2:615585/150‑1 (MQ=255) cGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGt < 2:377076/151‑1 (MQ=255) ggCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTtc < 1:54384/150‑1 (MQ=255) gggTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCa < 2:367937/151‑1 (MQ=255) ggACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCAtctc < 1:932092/151‑1 (MQ=255) ggACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCAtct < 1:1080125/150‑1 (MQ=255) | CGGAACGTGATTGCGGTGGCAATCGAACAGGCCGTTGATTAAGTGCAGGTTGCCGGGGCCGCACGATCCGGCGCAGACCGCCAGTTCTCCGCTAAGTTGTGCTTCAGCGCCAGCGGCAAAGGCCGCCACTTCTTCGTGGCGGGTGGACATCCACTCGATGGTGCCCATGC‑GATTAAGACTGTCACTAAGACCGTTCAGAGAGTCGCCTGTGACTCCCCAGATGCGTTTCACCCCTGCCGATTCGAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTC > NC_000913/910768‑911062 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |