| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,153,573 | 0 | T | G | 60.0% | 17.0 / 72.4 | 50 | Y92D (TAC→GAC) | pabC | 4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex |
| Reads supporting (aligned to +/- strand): ref base T (3/17); new base G (30/0); total (33/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.16e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-10 | |||||||||||
GACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACA > NC_000913/1153433‑1153721 | gACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGGAC‑Att > 1:84324/1‑151 (MQ=255) ggATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGc < 2:801136/151‑1 (MQ=255) tCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGCAC‑ATTTAACa > 1:703891/1‑140 (MQ=255) cAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGcac‑aattaaaaacgg > 1:386812/1‑132 (MQ=255) aGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGGAGCAC‑ATTTaaaaaggggcc > 1:222030/1‑136 (MQ=255) ttGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTTAAAAGCGGACCGGCCACACgg > 1:863990/1‑150 (MQ=255) gatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑AATGAAAAGCGGACCGGCAACGCGGAGtct > 1:408145/1‑150 (MQ=255) gatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGCCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATtct < 1:477904/150‑1 (MQ=255) ccTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGt < 2:559856/151‑1 (MQ=255) tgtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGCAC‑AATGAACAGAGGACCGGCCACGCGGGGTTTCTcc > 1:349202/1‑147 (MQ=255) tgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAGAGCAC‑AATGAACAGCGGGCACGGAACGCGGGGtctctc > 2:215250/1‑144 (MQ=255) gACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACg < 2:163092/150‑1 (MQ=255) gACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCAC‑ATTTAACAACGGGCCGGGGAAGCGGGTTCTCTCCCTTTCgg > 1:678304/1‑151 (MQ=255) aCTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCAC‑ATTTAAAAACGGACCCGGGACGCGGGGTCTCTCCCTTTCg > 2:519477/1‑149 (MQ=255) aCTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑AATGAACAGCGGACCGGCAACGCGGGGTCTCTCCGTTTCgg > 2:565608/1‑150 (MQ=255) cTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAAAGCAA‑AATTAAAAGCGGGCCGGGCACGCGGGGTTTCTCCGTTTCgggctattc > 2:711102/1‑141 (MQ=255) cAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGCAGCCC‑ATTTAACAGCGGAGCGGGAACACGGGGTTTCTCCCCTTCTGGTTATCCTg > 1:594447/1‑148 (MQ=255) tGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGACAGCCC‑AATGAAAAGGGGACCGGGGACCCGGAGTCTCTCCGTTACGGCTTATCCTGCACAtt > 2:674334/1‑149 (MQ=255) aaCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGGCAGGAC‑AATGGAAAGCGGACCGGGGACGCGGGGTCTCTTCCTTTCCGCTTAttatccacaa > 2:216926/1‑137 (MQ=255) aaCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGa < 2:691932/151‑1 (MQ=255) agagATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCAC‑AATGTACAACGGACCGGGGAAGCGGGGTCTCTCCCTTTTGGGTTTTTCTGCACATTACg > 2:607029/1‑146 (MQ=255) aTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTTAACAGCGGACCCGCAACGCGGGTTTTCTCCCTTTCGGCCTTTTCTTCCCATTACGAcccgtt > 1:836484/1‑145 (MQ=255) aCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACg < 2:252538/151‑1 (MQ=255) cGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGGACAATTGAA‑AACGGGCCGGCAACGCGGGGTTTCTCCGTTTCGGCTTATTCTGCACATTAAGACCGTTTGCg > 2:181455/1‑145 (MQ=255) gCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAg < 2:261736/151‑1 (MQ=255) agcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGCAGGAA‑ACTGAAAAGCGGACCGGCCACGCGGGGTCTCTCCCTTTCGGCCTTTCCTGCACACTACGCCCGCTTGCGTAAAGAgggggg > 2:190262/1‑148 (MQ=255) agAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGGAGCAC‑ATTGAAAAACGGACCGGCAACCCGGGGTCTCTCCGCTATGGCTTATTCTCCACAATACGACCGCTTGCGTGAAGAGGGGGTTACg > 2:739910/1‑151 (MQ=255) gAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGt < 1:508103/151‑1 (MQ=255) aaCAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGAAC‑ATTGAAAAGCGGACCGGCAACGCGGGGTTTCTCCCTTACGGCTTATCCTGCACATTAAGACCGGTTGCGTAAAGAGGGGGTTACGtt > 1:499824/1‑151 (MQ=255) aGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGgc < 1:216949/151‑1 (MQ=255) tGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGCAC‑AATTAAAAACGGGCCGGCAACGCGGGTTCTCTCCGTTTCGGCTTATTCTGCACATTACGACCGCTTGCGTGAAGAGGGGGGTACGTTTGCGC‑Caa > 2:582587/1‑149 (MQ=255) ggTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGc < 2:354261/150‑1 (MQ=255) tgtgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGACAGCAC‑AATTTACAGCGGGCCGGCCACCCGGGTTCTCTCCCCTTCCGGTTTTCCCTCCCAATACCACCCTTTTCGTAACGGGGGGGTTAAGtttggggcaaacc > 1:756868/1‑138 (MQ=255) gtgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTTAACAGCGGGGCGGCAACGCGGGTTCTCTCCGTTACGGCTTATTCTGCACATTACGACCGGTTGCGTAACGAGGGGGGTTCGTTTGCGC‑TAAGCCCgg > 1:330117/1‑151 (MQ=255) gtgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCg < 1:120122/150‑1 (MQ=255) aaGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTGAAAAGCGGACCGGCAACGCGGGTTTTCTCCCTTTCCGCTTTTTCTGCACATTACGACCGTTTTCGGAACGAGGGGGTTACGTTGGCGC‑TAAACCCggggggg > 1:66803/1‑144 (MQ=255) tCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGGAC‑ATTTAAAAACGGACCCGGAACGCGGGGTCTCTCCGTTTCGGCTTATCCTGCACCTTAAGACCGGTTGCGTGACGAGGGGGGTACGTTGGCGC‑TGAGACCCGGGCGGCgggg > 2:110138/1‑151 (MQ=255) gTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGgc < 1:611026/151‑1 (MQ=255) gTGATCAGTCGCGGTAGTGGCGGGCGAGGGGATAGGAC‑ATTGTACAGCGGGCCGGCAACACGGGGTCTCTCCCCTAATGGTTTTCCTGGACCATACGACCGTTTGCGTTACGAGGGGGGTACGTTTGGGCGTAGGCCCgggcgggttgggg > 1:672844/1‑140 (MQ=255) aTCAGTCGCGGTAGTGGCGGGGGGGGGGAGAGCAC‑ATTTAAAAACGGACCGGCAAAACGGGGTCTCTCCCCTTCGGCTTTTTCTGCACACTACGACCGCTTTCGTAACGAGGGGGGTATGGTTGCGC‑GCAACCCCGGGCGGCGgggg > 2:566777/1‑147 (MQ=255) aTCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCa < 2:597076/151‑1 (MQ=255) cAGTCGCGGTAGTGGCGGGCGGGGGGACAGCAC‑ATTGAACAACGGGACGGCAACGCGGGTTCTCTCCGGTACGGCCTATCCTGCACATTACGACCGTTTGTGTAACGAGGGGGTTACGTTGGCGC‑TAAGCCCCGGGGGGCTgggggg > 1:83204/1‑145 (MQ=255) aGTCGCGGTAGTGGCGGGCGGGGGGACAACAC‑ATTGAAAAGCGGGCCGGCCACGCGGGGTTTCTCCCCTTTGGCCTATTCTTCCCAATACGACCGGTTTCGGAAAGAGGGGGGTTCGTTGGCGC‑CAAGCCCCGGGCGGGGGGGgggcga > 1:659573/1‑147 (MQ=255) cgcgGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCAATCCTc < 2:607135/151‑1 (MQ=255) tGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGc > 2:277104/1‑150 (MQ=255) ggcgggcgAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGc > 2:357622/1‑149 (MQ=255) ggcgggcgAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAg < 1:441169/151‑1 (MQ=255) gAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTa > 1:598403/1‑150 (MQ=255) aGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTaaa < 1:399579/151‑1 (MQ=255) ggTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACa < 2:289274/151‑1 (MQ=255) | GACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCAC‑ATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGC‑TAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACA > NC_000913/1153433‑1153721 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |