| 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.1 / 13.0 | 44 | I441V (ATC→GTC) | mpl | UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase |
| Reads supporting (aligned to +/- strand): ref base A (1/24); new base G (19/0); total (20/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.42e-11 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACG > NC_000913/4456961‑4457244 | ccGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAt < 2:413844/150‑1 (MQ=255) ggcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCAT‑c > 2:367524/1‑147 (MQ=255) ggcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGaa < 1:49585/151‑1 (MQ=255) caggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGTGGGGGCCAT‑CAGaaa > 1:820044/1‑150 (MQ=255) aggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAAct < 2:865998/151‑1 (MQ=255) ggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCAT‑TAGaaa > 2:159540/1‑148 (MQ=255) gaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCAT‑CAGAAAAAGCTGGAt > 2:218879/1‑150 (MQ=255) agtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCAT‑CAGAAACTGCtgggtggg > 2:463064/1‑151 (MQ=255) agtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATggg < 2:383096/151‑1 (MQ=255) gaggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCGTTGTGGGGGTCCATACAGAAACGACTGGagg > 1:681680/3‑147 (MQ=255) aagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCgaa < 2:77576/151‑1 (MQ=255) cGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGGGGGGGCCAT‑CAGAAACTGCTGGGTGGGGTGGGGGAGAAAGCg > 1:71947/1‑151 (MQ=255) tCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGa < 1:801891/150‑1 (MQ=255) cAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGGCCCT‑CAGAAAATGCTGGATGGGGGGGCggagaaaggg > 1:382751/1‑143 (MQ=255) cAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGaa < 1:719804/150‑1 (MQ=255) gCCTGCACACTGGAGTGGCGAGGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGcc < 1:536352/151‑1 (MQ=255) ccTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCCT‑CAGAAAATGCTGGGTGGGGTGGCggagaaagggg > 1:383628/1‑140 (MQ=255) cacaCTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCa < 1:749131/151‑1 (MQ=255) aGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCAT‑CAGAAACTGCTGGATGGGGTGGCGAAGAAGGCGGGAGCCGCGCAGGAAtt > 2:834465/1‑149 (MQ=255) gCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCt < 2:395454/151‑1 (MQ=255) gCGATGTGGATACGCTGGCAGATATGGGGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGGGGGGTCCCT‑CAGAAACTGCTGGGTGGGGTGGCGAAGAAAGCGGAAGGCGCGCAAGAATTTgg > 1:154624/1‑148 (MQ=255) gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCAT‑CAGAAAATGCTGGATGGGGTGGGGAAGAAGGGGGGAACCGCGCAGTAAATCgggcct > 1:474461/1‑146 (MQ=255) gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCtt < 1:319875/150‑1 (MQ=255) aTGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAg < 2:118697/151‑1 (MQ=255) ggATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAg < 1:513164/151‑1 (MQ=255) aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGatat < 2:600711/151‑1 (MQ=255) aCGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGata < 2:569151/150‑1 (MQ=255) gCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAg < 1:351633/151‑1 (MQ=255) gCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAg < 2:84203/151‑1 (MQ=255) gCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATg < 1:118305/151‑1 (MQ=255) ggtgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGa < 2:233521/151‑1 (MQ=255) tgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCAT‑CAGAAAATGCTGGATGGGGTGGCGAAGAAGGCGGGAGCCGCGCCGTAATTCGGCCCTTGCCAGATATAGTATGCGTAGGcc > 2:208192/1‑146 (MQ=255) tgAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCAT‑CAGAAACTGCTGGATGGGGTGGCGGAGAAGGGGGGAGCCGCGCCGGAATTCGGCCTTAGCCCGATATAGTATGCGTAGGCCgggt > 1:111697/1‑148 (MQ=255) aaaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCAT‑CAGAAACTTCTGGGTGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGGAATTCGGCCTTAGCCAGGTATAGTATGCGGAGGCCGGATa > 2:718821/1‑149 (MQ=255) cAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCa < 1:679012/149‑1 (MQ=255) gCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGGCCAT‑CAGAAAATGCTGGGTGGGGTGGCGGAGAAGGGGGGAGCCGCGCAGGAATTCGGCCTTAGCCAGGTATAGTATGCGTAGGCCGGGTAAGGCGTTCAcgcc > 1:465194/1‑151 (MQ=255) ccTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCAT‑CAGAAAATGCTGGGTGGGGCGGCGGAGAAAGCGGGAGCCGCGCCGGAAATTGGCCTTAGCCAGGTATAGTATGCGTAGGCCGGGTAAGGCGGTCAcgcc > 1:574753/1‑150 (MQ=255) cGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGGGGGGTCCCT‑CAGAAACTGCTGGGTGGGGTGGGGAAGAAAGCGGGAGCCGCGCAGTAAATTGGCCTTAGCCAGATATAGTATGCGTAGGCCGGGTAAGGGGTTCACGCCGCCTc > 1:637476/1‑150 (MQ=255) ggTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCa < 2:379011/151‑1 (MQ=255) gTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCAt < 2:730722/151‑1 (MQ=255) tgatgaGCAACGGCGGGTTTGGGGGGGTCCAT‑CAGAAAATGCTGGGTGGGCCGGCGCAGAAGGCGGGAGCCGCGCAGAAATTCGGCCTTTGCCAGGTATAGTATGCGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAATTGAGc > 1:711372/1‑148 (MQ=37) gatgaGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATa < 2:286419/150‑1 (MQ=255) ggcggTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGAt < 1:392605/149‑1 (MQ=255) gTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACg > 1:75315/1‑151 (MQ=255) | CCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCAT‑CAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACG > NC_000913/4456961‑4457244 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |