Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,421,437 | 0 | T | C | 43.2% | 47.0 / 12.2 | 37 | I38T (ATC→ACC) | folX | D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase |
Reads supporting (aligned to +/- strand): ref base T (18/3); new base C (0/16); total (18/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.71e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATG‑TGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGC > NC_000913/2421287‑2421584 | tattatTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATG‑TGACGAt > 2:252448/1‑151 (MQ=255) aGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATG‑TGACGATCCACTACCCCGCCGa < 1:802165/151‑1 (MQ=255) aacaacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATG‑TGACGATCCACTACCCCGCCGATaa > 2:738644/1‑149 (MQ=255) caacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATG‑TGACGATCCACTACCCCGCCGATAAAgcg > 2:375847/1‑151 (MQ=255) cGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATG‑TGACGATCCACTACCCCGCCGATAAAgcgcgc > 2:580682/1‑151 (MQ=255) gagaACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATG‑TGACGATCCACTACCCCGCCGATAAAGCGCGCAc > 1:707243/1‑151 (MQ=255) ccTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATG‑TGACGATCCACTACCCCTCCGATAAAGCGCGCACTAg > 1:15957/1‑149 (MQ=255) cTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAAGATATTGTTATCAATG‑TGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCg > 2:190029/1‑150 (MQ=255) acacCCTGCCGTTATTATTCGAAAAAAGACCCTTCGTTTGCGTACGTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGAAAAAATATTGTTATAAAGG‑AGACACCCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATc < 2:21853/148‑1 (MQ=255) cccTGCCGTTATTATTAGTAAAAAAAACCTTCGTTTGCGTACTTTTATCGAAATAAAGAAAGAAGAAATAAACCCCCGTCAAGATTTTTTAATCAATG‑TGCCACTCCACGCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAAc < 2:435875/147‑1 (MQ=255) cccgCCGTTATTATTCGTAAAAAGACCCTTCGTTTGCGTACGTTTATCGAAATTAAGAAAAAAAAAATAAACACCCGTCAGGATTTTGTTAAAAAGG‑TGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATg < 1:616118/148‑1 (MQ=255) cgtttattatTCGTAAAAAGACCCTTCTTTTGCGTACGTTTATGGGAATAAAGGAAGAAAAAATAACCACCCGTAAGGATATGGTTATCAAGG‑TGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCt < 1:394901/147‑1 (MQ=255) ttattCTAAAAAAGAACCTTCTTTTGCGTACGTTTATGGAAATAAAGAAAAAAAAAATAAACACCCGTCAGAATTTTGTTATAAAGG‑TGACCACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGaa < 1:677593/148‑1 (MQ=255) ttCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATG‑TGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg > 1:743445/1‑151 (MQ=255) cgaaaaaaGAACCTTCGTTGGCCTCCTTTTACCGAAAGAAGGAAAGAAAAAATTAACACCCGTCAGGATTTTGTTATAAATG‑TGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATc < 1:293032/143‑1 (MQ=255) tataAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATG‑TGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAcc > 1:217008/1‑151 (MQ=255) aaaaaaaaCCTTTTTTTGCGTATGTTTATGGAAATTAAGAAAAAAAAAATAAACACCCGTCAGGATATTTTTATAAATG‑TGACATCCCACACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATc < 2:430051/143‑1 (MQ=255) cTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATG‑TGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa > 2:758702/1‑150 (MQ=255) ttCTTTTGCGTACTTTTATGGAAATAAAGAAAAAAAAAATAAACAACCGTCAGGATTTTTTTATAAATG‑TGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaa < 1:320314/148‑1 (MQ=255) ttCTTTGGCGTTCGTTTATGGAAATAAAGAAAAAAAAAATAACCACCCGTCAGGATTTTTTAATAAAGG‑TGCGGCCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaa < 2:833469/148‑1 (MQ=255) ttCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATG‑TGACGATCCACTACCCCGCCGATAACGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa > 2:194905/1‑150 (MQ=255) ttCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATG‑TGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa > 2:4976/1‑150 (MQ=255) cGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATG‑TGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc > 2:70823/1‑150 (MQ=255) gTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATG‑TGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc > 2:90368/1‑150 (MQ=255) aTCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATG‑TGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGGTATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt > 2:188151/1‑149 (MQ=255) aaTAAAGAAAGAAGAAATAAACACCCGTCAGGATTTTGTTAAAAATG‑TCACGACCCCCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt < 1:744107/144‑1 (MQ=255) aaTAAAGAAAAAAAAAATAACCACCCGGAAGGATATTTTTATAAATG‑TGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGaa < 2:500396/150‑1 (MQ=255) taaaGGAAGAAAAAAAAAACACCCGTAAGGATTTTTTTAT‑AATGTTGACAACCCCCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATa < 2:213186/148‑1 (MQ=255) aaaaaaTAAACACCCGTCAGGATTTTTTTACAAATG‑GAACGCCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtctc < 1:216823/148‑1 (MQ=255) aacaccCGTCAGAATTTTTTTATAAATG‑TGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg < 2:707046/151‑1 (MQ=255) aacCGTCAGGATATTGTTATCAATG‑TGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaa > 1:814031/1‑150 (MQ=255) ccctccAGGATTTTTTTAACAATG‑TGACCCCCCCCACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGc < 1:582347/139‑1 (MQ=255) ttttATCAATG‑TGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGa < 1:580832/148‑1 (MQ=255) aTCAATG‑TGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGc > 2:179800/1‑151 (MQ=255) tCAATG‑TGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCt > 1:625407/1‑151 (MQ=255) aaTG‑TGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTc < 1:190064/150‑1 (MQ=255) tgACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTc < 1:278521/146‑1 (MQ=255) gATCCACTACCCCGCCGATAACGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGTATGAGCTCGATACCGc > 2:303102/1‑150 (MQ=255) | TATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATG‑TGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGC > NC_000913/2421287‑2421584 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |