Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,457,105 | 0 | A | G | 43.2% | 58.3 / 12.7 | 45 | I441V (ATC→GTC) | mpl | UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase |
Reads supporting (aligned to +/- strand): ref base A (3/22); new base G (19/0); total (23/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.19e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTG > NC_000913/4456961‑4457248 | ccGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATc < 2:86529/151‑1 (MQ=255) gcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGGTT‑TGGGGGGGGCCCTCAGaa > 1:345055/1‑150 (MQ=255) aggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGGTT‑TTGTGGGGTGCATAAGaaa > 1:562878/1‑149 (MQ=255) agaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGa < 2:1608/150‑1 (MQ=255) aagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATgg < 2:720374/151‑1 (MQ=255) gtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGc < 1:492426/151‑1 (MQ=255) ggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCt < 1:708838/150‑1 (MQ=255) gcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCTgg < 2:574983/151‑1 (MQ=255) gaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCga < 1:874103/151‑1 (MQ=255) gCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCgaaga < 2:576007/151‑1 (MQ=255) tGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTT‑TGGGGGGGTCC‑TTAGGAACTGCTGGGTGGGGTGGCGAAGAAgg > 2:563779/1‑151 (MQ=255) gCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGaagaag < 2:674365/149‑1 (MQ=255) cGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TTGTGGGGCCCACCAGAAAATGTCGGGGGGGGTGGCGCAGAAgggg > 1:808727/1‑149 (MQ=255) cAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGGTT‑TGGGGGGGGGCATCCGAAAAAGCTGGGTGGGGTGGCGAAGAAAGCggg > 1:132757/1‑148 (MQ=255) gCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGGTCCCACCGAAAATGCTGGGTGGGGTGGGGGAGAAAGCGGAAGcc > 1:892202/1‑151 (MQ=255) tGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TTGGGGGGTCCATCAGAAAATGCTGGATGGGCTGGGGAAGAAGGCGGGAGCCGCGCAGTAATTCgg > 2:748481/1‑150 (MQ=255) cGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCtt < 2:130296/151‑1 (MQ=255) gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGGGCCCTCAGAAACTGCTGGATGGGCGGGGGAAGAAGGCGGGAGGCGCGCAGGAATTCGGccct > 1:371291/1‑148 (MQ=255) gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTa < 1:699319/151‑1 (MQ=255) aTGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TTGGGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAAGCGGGAGCCGCGCAGTAATTCGGCCtttg > 2:533625/1‑149 (MQ=255) ggATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TTGGGGGGGCCCACCGAAACTGCTTGAGGGGGTGGCGGAGAAGGGGGGAGCCGCGCAGGAATTCCGCCTTAGcc > 1:294585/1‑149 (MQ=255) aTACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGat < 1:662928/151‑1 (MQ=255) tACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TTGTGGGGTCCATCAGAAACTGCTGGGTGGGGTGGGGAAGAAGGGGGGAGGCGCGCCGTAATTTGGCCTTAGCCAGat > 2:571585/1‑150 (MQ=255) tACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGGTCCATCAGAAAATGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGGAATTCGGCCTTAGCCAGat > 2:83252/1‑150 (MQ=255) gCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAg < 2:758775/151‑1 (MQ=255) gCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAg < 2:614818/151‑1 (MQ=255) aGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCg < 2:52052/151‑1 (MQ=255) aaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGGGCCCTCAGAAACTGCTGGGTGGGCTGGCGGAGAAGGCGGGAGCCGCGCCGTAAATTGGCCTTAGCCCGATATAGTATGCGTAGGCCgggtta > 1:194481/1‑144 (MQ=255) ccGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGGTCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGGGGAAGCCGCGCAGTAATTCGGCCCTTGCCCGATATAGTATGCGTAGGCCGGGTAAgg > 1:762953/1‑148 (MQ=255) cGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGGTCCATCAGAAAATGCTGGGTGGGGTGGCGGAGAAAGCGGGAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAATATGCGTAGGCCGGATAAgg > 2:338442/1‑147 (MQ=255) gCTCAGCCTGGCGACCATATTCTGGTGATGTGCAACGGC‑GGTTT‑TGGTGGGGTCCATCAGAAACTGCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAgg > 2:400699/1‑146 (MQ=255) gCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGGTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCTGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTTAGGCggt > 2:697479/1‑148 (MQ=255) tCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGGTCCCTCAGAAAATGCTGGGTGGGGTGGGGGAGAAGGGGGGAGCCGCGCAGGAATTCGGCCTTTGCCAGATATAGTATGCGTAGGCCGGGTAAGGGGGTCa > 1:723227/1‑150 (MQ=255) cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCAcg < 1:857466/151‑1 (MQ=255) gCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTTATGCTGGG‑TGCCCCAAACCCCCCTGGCTGGCCCGGCGCCGAAAGCGGCAGCCGCCCAATAATTCGGCCTTCGCCAGAGATTGGGTGCGTGGGACGGGGAACGCGTTCAc > 1:284219/1‑148 (MQ=255) tGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TTGGGGGGTCCCTTAGAAAATGCTGGGTGGGGTGGGGAAGAAGGGGGGAGCCGCGCAGGAATTCGGCCTTAGGCCGATATAGTATGCGTAGGGCGGGTAAGGGGGTCAcgcc > 1:198567/1‑148 (MQ=255) gCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TTGGGGGAACCCGCAAAAACTGCTGGATGGGCGGGCGAAGAAGGGGGGAGCCGCGCAGGAATTCGGCCTTAGCCAGGTATAAGATGCGTAGGCCGGATAAGGGGGTTAcgccgactc > 1:219392/1‑147 (MQ=255) gACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg < 2:816369/151‑1 (MQ=255) ttCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAttt < 1:580519/150‑1 (MQ=255) ttCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAAATGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAttt > 2:195562/1‑150 (MQ=255) cTGGTGATGAGCAACGGC‑GGGTT‑TGGGGGGGGTCCTCAGAAACTGCTGGGAGGGGCGGCGGAGAAGGGGGAAGCCGCGCCGTAATTCTGCCTTAGCCAGGTATAGTATGCGTAGGCCGGGTAAGGCGGTCACGCCGCATCCGGCAttt > 1:9898/1‑148 (MQ=255) cAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTg < 1:218817/151‑1 (MQ=255) ggTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGAc < 2:563026/151‑1 (MQ=255) tt‑tGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGc < 2:367902/150‑1 (MQ=255) ggtggGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTg < 2:191783/150‑1 (MQ=255) | CCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGC‑GGTTT‑TGGTGGGATCCATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTG > NC_000913/4456961‑4457248 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |