| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,153,573 | 0 | T | G | 48.0% | 42.9 / 43.0 | 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 (1/25); new base G (24/0); total (25/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.11e-13 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAA > NC_000913/1153431‑1153719 | gCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGcaca < 1:388115/151‑1 (MQ=255) cAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACa < 2:1044755/151‑1 (MQ=255) cAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGGACAATGaa > 2:205923/1‑149 (MQ=255) aGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGcaaa > 2:1011545/1‑141 (MQ=255) aTGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGAGAGCAAATTggaaaacgg > 1:947707/1‑143 (MQ=255) gTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGc < 1:836792/151‑1 (MQ=255) aGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGTGTCCAGCACATTGAACAGCGGACCGGCAAc < 1:849066/151‑1 (MQ=255) cGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGCACAATTAAAAACGGACCGGCaa > 2:281938/1‑148 (MQ=255) cGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACCTTgcaaacggggcc > 1:441481/1‑132 (MQ=255) ggTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAAcgcg < 2:825642/151‑1 (MQ=255) ttGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAAAGCACATTTAACAGCGGAGCGGGAACGCggg > 2:176662/1‑150 (MQ=255) tgatgaTTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGAAGCACACTGAAAAGCgg > 2:566837/1‑136 (MQ=255) atTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGGAGCACACTGTACAGCGGAACGGGAACCCggggtttc > 1:215040/1‑144 (MQ=255) tCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGGAGCACATTTGAAAACGGGCCGGCGACGCggggtttc > 1:732205/1‑141 (MQ=255) tgtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGCACACTTTAAAGCGGGGCGGGGACGCggggtgt > 2:94899/1‑138 (MQ=255) gtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGAGAGCACATTAGAGAGCGGACCGGCCCCGCggggggtcccccttgc > 1:498869/1‑137 (MQ=255) gACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACg < 1:288054/150‑1 (MQ=255) gACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACg < 1:476449/150‑1 (MQ=255) cTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCACAATGAACAACGGGCCCGCAACGCGGGGTCTCC‑CCC‑TTTCGG‑c > 2:995699/1‑150 (MQ=255) tCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑Ctt < 2:955416/150‑1 (MQ=255) cTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAACAAAATTAACAACCGGACGGCAACGCCGGATCTCTCCCG‑TTACg > 1:1078132/1‑146 (MQ=255) ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATc < 2:594929/150‑1 (MQ=255) gCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTAt < 1:545972/148‑1 (MQ=255) ccTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTg < 1:582284/151‑1 (MQ=255) tCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGGAGAAAAATTAAAAACGGACCGGCCACGCGGGGTCTCT‑CCC‑tttaggggtattc > 1:964399/1‑136 (MQ=255) gCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTGCACAt < 2:555483/151‑1 (MQ=255) gAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTGCACAtttc < 1:820271/150‑3 (MQ=255) agagATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTGCACATTACGACCGt < 1:418350/151‑1 (MQ=255) agaTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTGCACATTACGACCGttt < 1:1086084/151‑1 (MQ=255) gaTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGGGAGCACAATGGAAAACGGAGCGGGGAAGCGGGTTTTCT‑CCC‑TTTCGG‑CTTTTTCTGccccattcgacccgtt > 1:376153/1‑135 (MQ=255) aaaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGACAGCACACTTAAAAGCGGGGCGGCCACGCGGAGTCTCT‑CCC‑TTTCGGCCTTTTCCT‑CACACTACgcac > 1:334371/1‑138 (MQ=255) aaaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTGCACATTACGACCGTTTGCGTaa < 2:995714/151‑1 (MQ=255) aaaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTGCACATTACGACCGTTTGCGTaa < 1:303771/151‑1 (MQ=255) tGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGGACAATTAACAGCGGAGCGGCGAAGCGGGTTCTCT‑CCG‑CTTCCG‑CCTATTCTTCAAATTACGAACGTTTTTGTGACGAggg > 1:1007944/1‑151 (MQ=255) ggCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTGCACATTACGACCGTTTGCGTAACGAggg < 2:580617/150‑1 (MQ=255) cagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCACAATGAAAAACGGACCCGCAACGCGGGTTCTCT‑CCC‑GTTCGG‑CTTATCCTGCCCCATACGACCCTTTTCGGAACGAGGGGGtt > 2:664670/1‑149 (MQ=255) agAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACg < 1:259434/151‑1 (MQ=255) gCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGGCCAATTTAAAGCGGGCCGGCAACGCGGGGTCTCT‑TCC‑TTTCGG‑GTTATTCTTCCCATTATGCACGCTTTCGTTAAGAggggggtaactt > 2:707651/1‑139 (MQ=255) gCTTAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTg < 2:974073/151‑1 (MQ=255) cTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAAAGCACAATTAAAAACGGACCGGCCACGCGGGGTTTCT‑CCG‑CTTCGG‑CTTTTTCTGCACATTACGACCGGTTTCGTAAAGAGGGGGGTACG‑TTTGGGCCAAGCCCgg > 2:613485/1‑148 (MQ=255) aaGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCt < 2:1031622/151‑1 (MQ=255) cGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGCACAATGAACAGGGGGCCGGCGACGCGGGGTGTCT‑CCG‑GTTCGG‑CTTATCCTGCCCACTTCGACCCGTTGTGTGACGAGGGGGTTACG‑TTTGGGCCAAGCGCCGGGGGGGTgggg > 1:111646/1‑151 (MQ=255) cGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCTggg < 2:1075908/150‑1 (MQ=255) gTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGCACAATGAAAAGCGGGCCGGCAAAGCGGGGTTTCT‑TCGCTTTCGG‑CTT‑TTCTGCCCCATACTAACGGTTTCGTTAAGAGGGGGGTAAG‑TTGGGGGTAAGACCCGGGCGGCTgg > 2:693471/1‑148 (MQ=255) gTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGGACATTGAAAAGAGGGCCGGGCACGCGGGGTTTCT‑CCG‑GTACGG‑GTTATACTGCACATTACGACCGGTTTCGTAAAGAGGGGGGTACG‑GTGGGGCCAAGGCCCGGGCGGCGggggg > 2:397196/1‑150 (MQ=255) cAGTCGCGGTAGTGGCGGGCGAGGGGGCAGCACCTTGAAAAGCGGGGCGGCAAAGCGGGGTCTCT‑CCC‑TTTCGG‑GCTTTCCTGCACATTAACACCGGTTTTGGAAAGAGGGGGTTTGGTTTGGGGTTAAGCCCgggggggggggggg > 2:821524/1‑135 (MQ=255) cgcgGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTc < 1:472184/151‑1 (MQ=255) gcgGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCa < 1:770587/151‑1 (MQ=255) gTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCt < 2:590595/151‑1 (MQ=255) ggTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCCGGGGCGCAATCCTCATCTTGCAGGTATTaaa > 1:926067/1‑149 (MQ=255) | GCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCT‑CCG‑TTACGG‑CTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACG‑TTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAA > NC_000913/1153431‑1153719 |
| 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: ‑ |