| 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% | 40.1 / 10.5 | 37 | I441V (ATC→GTC) | mpl | UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase |
| Reads supporting (aligned to +/- strand): ref base A (9/12); new base G (16/0); total (25/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.09e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTG‑GCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGC‑CAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTG > NC_000913/4456956‑4457248 | atatTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGAt < 2:589175/151‑1 (MQ=255) ggcaggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGa < 1:263071/150‑1 (MQ=255) ggtggcagAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAActg < 1:557770/151‑1 (MQ=255) cagaagtggcagaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGGGGGGTCCCTCAGAAACTGCTgg > 2:448050/1‑150 (MQ=255) gcagaagCCTGCGGTCATCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGCGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGATTTTGGGGGAGTCATCCGGAAATGCTGGGAGGGCt > 2:302379/1‑149 (MQ=255) gaagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTG‑GCga < 1:45538/151‑1 (MQ=255) aagCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAAATGCAGGGTGGGGTG‑GCg > 2:513284/1‑149 (MQ=255) tGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGGTGGGGTG‑GCGAAGAAgg > 1:511956/1‑151 (MQ=255) tGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGGGCCCACAGAAACCGCTGGATGGGGTG‑GCggggaa > 2:492/1‑144 (MQ=255) gCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTG‑GCGAAGAAGGc < 2:28916/151‑1 (MQ=255) cAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGAGTCATCCGAAAATGCTGGATGGGCCG‑GGgcagaaagggg > 1:554469/1‑143 (MQ=255) gCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTG‑GCGAAGAAGGCGGAAGCCgcgc < 1:119352/151‑1 (MQ=255) acacTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTTGGGGGGTCCATCAGAAACCGCTGGGGGGGCTG‑GCGAAGAAGGCGGAAGGcgcg > 1:124096/1‑148 (MQ=255) ggAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGTGGGGGCCATCAGAAACTGCTGGATGGGCTG‑GCGAAGAAAGAGGAAGCCGCGCAGTaa > 2:549058/1‑149 (MQ=255) gATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATTCAACAGAAACTGCTGGAGGGGGTG‑GGGGAGAAGGCGGGAGCCGCGCAGGAATTCGGGCt > 2:510235/1‑149 (MQ=255) ggATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGCTTTGGGGGGGGCCATCAAAAACTGCTGGAAGGGGCT‑GCGGAGAAGGCGGGAGCCCCGGCGGAAATCGGCCCTAg > 2:607682/1‑147 (MQ=255) gCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGTGGGGTCCCTCAGAAACTGCTGGGTGGGCTG‑GCGGAGAAAGCGGGAGCCGCGCCGTAAATCGGCCTTAg > 2:170294/1‑141 (MQ=255) cTGGCAGATATGGGGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGGTTTGGTGGGATCCCTCAGAAACTGCTGGATGGGCTG‑GCGAAGAAGGCGGGAGCCGCGCAGGAATTCCGCCTTCGC‑CAGATATAg > 2:165639/1‑150 (MQ=255) tGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTG‑GCGAAGAAGGCGGAAGCCGCGCAGGAATTTGGCCTTAGC‑CAGATATAGt > 1:282258/1‑150 (MQ=255) ggCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTG‑GCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGC‑CAGATATAGTa < 2:350939/150‑1 (MQ=255) aaaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTG‑GCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGC‑CAGATATAGTATGCGTAGGCCGGATAAg < 2:92335/151‑1 (MQ=255) aaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCCTCAGAAACTGCTGGGTGGGCTG‑GCGGAGAAGGCGGGAGCCGCGCAGTAATTCGGCCTTAGC‑CAGATATAGTATTCGTCGGCCGGGTAAgg > 2:85176/1‑151 (MQ=255) aaaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTG‑GCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGC‑CAGATATAGTATGCGTAGGCCGGATAAgg < 1:396897/151‑1 (MQ=255) aaCCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGCTG‑GCGAAGAAGGCGGGAGCCGCGCAGTAATTCGGCCTTAGC‑CAGATATAGTAAGCGTAGGCCGGATAAgg > 1:259274/1‑150 (MQ=255) cGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAAATGCTGGATGGGCTG‑GCGAAGAAGGCGGAAGCCGCGCAGGAATTCGGCCTTAGC‑CAGATATAGTATGCGTAGGCCGGATAAGGCggt > 1:95350/1‑149 (MQ=255) aGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTG‑GCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGC‑CAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCAcgc > 1:156809/1‑151 (MQ=255) cTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGGGGGGTCCATCAGAAACTGCTGGATGGGGTG‑GCGGAGAAGGCGGGAGCCGCGCAGTAATTCGGCCTTAGC‑CAGATATAGTATGCGTAGGCCGGGTAAGGCGGTCAcgcc > 1:58996/1‑149 (MQ=255) gACCATATTCTGGTGATGAGCAACGGCGGGTTTTGTGGGATACATCCGAAAATGCTGGATTGGTAGTGGAAAAAAGGCGG‑AGACGCGGCTTAATTTGGGCCTAGCTCAGTAAAAGAAT‑CGTCGGCACGATAAGGGGTTCACGCCGCAt > 1:249939/1‑148 (MQ=255) aCCATATTCTGGTGATGAGCAACGGCGGTTTTTGTGGGGTCCCACAGAAACTGCTGGATGGGCTG‑GCGAAAAAGGCGGAAGCCGCGCAGTAATTTGGCCTTAGC‑CAGATATAGTATGCGTAGGCCGGATAAGGCGTTTACGCCGCCTCCgg > 2:559784/1‑151 (MQ=255) cTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGGTG‑GCGGAGAAGGCGGAAGCCGCGCAGGAATTCGGCCCTAGC‑CAGATATAGTATGCGTAGGCCGGGTAAGGCGTTCACGCCGCATCCCGCAttt > 1:59967/1‑148 (MQ=255) cTGGTGATGAGCAACGGCGGTTTTGGTGGGGTCCATCAGAAACTGCTGGATGGGCTG‑GGGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGC‑CAGATATAGTATGCGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCAttt > 2:59074/1‑148 (MQ=255) ggcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTG‑GCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGC‑CAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATg > 2:97056/1‑150 (MQ=255) gcggTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTG‑GCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGC‑CAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCg < 1:513448/151‑1 (MQ=255) ttttGGTGGGATCCATCAGAAACTGCTGGATGGGCTG‑GCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGC‑CAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGc > 1:119547/1‑151 (MQ=255) tggtggGATCCATCAGAAACTGCTGGATGGGCTG‑GCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGC‑CAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTg > 1:371095/1‑151 (MQ=255) ggtggGATCCATCAGAAACTGCTGGATGGGCTG‑GCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGC‑CAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTg < 2:470392/150‑1 (MQ=255) ggtggGATCCATCAGAAACTGCTGGATGGGCTG‑GCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGC‑CAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTg < 1:575520/150‑1 (MQ=255) | ATATTCCGTGGCAGGTGGCAGAAGTGGCAGAAGCCTGCGTTCAGCCTGCACACTGGAGTGGCGATGTGGATACGCTGGCAGATATGGTGGTGAAAACCGCTCAGCCTGGCGACCATATTCTGGTGATGAGCAACGGCGGTTTTGGTGGGATCCATCAGAAACTGCTGGATGGGCTG‑GCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGC‑CAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTG > NC_000913/4456956‑4457248 |
| 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: ‑ |