| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,153,573 | 0 | T | G | 45.6% | 63.8 / 19.4 | 46 | Y92D (TAC→GAC) | pabC | 4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex |
| Reads supporting (aligned to +/- strand): ref base T (5/20); new base G (21/0); total (26/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.04e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAA > NC_000913/1153423‑1153718 | cATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGt < 2:704458/151‑1 (MQ=255) gATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCg < 2:696406/151‑1 (MQ=255) gCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCACA‑ATGAACAGCg > 1:461103/1‑148 (MQ=255) tGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGCAAA‑ATGAAAAGCgg > 1:647515/1‑146 (MQ=255) cAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCa < 1:342698/150‑1 (MQ=255) gTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCACA‑ATGAACAGCGGAGCGGCCACGCgg > 2:993356/1‑151 (MQ=255) atgatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAGAGCACA‑TTGAACAGCGGGCCGGGAACGCgg > 1:223486/1‑147 (MQ=255) atgatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACA‑TTTAAAAGCGGACCGGCaa > 1:1001038/1‑142 (MQ=255) tgatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCACA‑TTGAAAAACGGGACGGCGAAGCGGGttgt > 1:523663/1‑149 (MQ=255) ttCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGCACA‑TTGAAAAACGGAGCGGCAACGCGGGATTT‑CTcc > 1:450708/1‑151 (MQ=255) gACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑g < 1:979343/150‑1 (MQ=255) aCTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAACACA‑CTGAACAGCGGACCGGGAACGCCGGTTTT‑CTcccctt > 1:529038/1‑143 (MQ=255) cTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGGACA‑ATTTAAAACGGACCGGGAAAGGGGGGTCTCCTCCGttc > 2:574467/1‑146 (MQ=255) tCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCt < 1:353822/149‑1 (MQ=255) cAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGca < 1:504387/150‑1 (MQ=255) gCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGAAGCACA‑ATGAAAAAAGGGGCGGCAAAGCGGGGTTT‑CTCCCTTTC‑GGGCTATTCTGcaca > 1:876673/1‑150 (MQ=255) gCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACAt < 2:843149/151‑1 (MQ=255) cTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGAAGGACA‑TTTTACAGGGGGCCGGCGACGCGGAGTTT‑CTCCCTTTC‑GGCCTATTCtcccc > 1:262210/1‑144 (MQ=255) gAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACC‑ATTAACAACGGACCGGCAACGCGGGGTCT‑CTCCCCTTC‑GGCTTATcc > 1:185215/1‑140 (MQ=255) aaaaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTa < 2:355055/151‑1 (MQ=255) aaaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTaa < 1:60324/151‑1 (MQ=255) aaaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGGCAGCACA‑ATGGAAAGCGGGCCGGCCACGCGGGTTCT‑CTCCCTTAC‑GGCCTATTCTTCACCTTAccaac > 2:36642/1‑137 (MQ=255) cagcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACC‑TTTAAAAACGGGGCGGGCACGCGGGGTTT‑CTCCCTTTC‑GGCCTATTCTTCCCATTACCAACGTTTGCGGTAAGAggggg > 1:772874/1‑149 (MQ=255) agcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGAtt < 2:1038158/151‑1 (MQ=255) agcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGAtt < 2:274111/151‑1 (MQ=255) cagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAAAACACC‑ATGAAAAACGGGCCCGGGAAGCGGGGTCT‑CTCCGTTCC‑GGGTTATCCTCCCCAATTCGAACCGTTGGGGAACCAggggggta > 1:742476/1‑146 (MQ=255) agAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGGAACACA‑ATTAAAAGCGGGCCGGGAACGCGGGGTTT‑CTCCCTTTC‑GGCCTATTTTCCCCATTACGACCGCTTTGGGAACGAggggggt > 2:158382/1‑145 (MQ=255) aaCAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGtt < 1:137550/151‑1 (MQ=255) gCAAAATGGTGTGCTGAAAGTCGTGATCGGTCGCGGTAGTGGCGGGCGGGGGGACAGCACA‑TTTAAAAGCGGGACCGCCACCCGGGTTTT‑CTCCCTTTC‑GGCTTATTCTGCACCTTAAGACCGTTTTCGTGACGAGGGGGGTACGTTTgcg > 1:255287/1‑151 (MQ=255) gCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAACACACTTTAAAACGGGCCCGGCAA‑GCGGGGTTT‑CTCCCTTTC‑GGCTTTTTCTGCCCATTAAGACCGTTTTCGTAACGAGGGGGGTACGTTTgcg > 1:345753/1‑151 (MQ=255) cAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAAAACACA‑ATTAACAACGGGCCCGCGACACGGATGCT‑CTTCGTTACTGGCTT‑TCCCTCCCATTACAACCGGTTGCGTAAGGAggg > 1:880071/1‑137 (MQ=255) aaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCT‑aa < 2:468512/151‑1 (MQ=255) gtgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCT‑AAGCCCgg < 1:633978/151‑1 (MQ=255) tGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACA‑TTGTAAAGCGGACCGGCAACGCGGGGTTT‑CTCCGCTAC‑GGCCTATTCCGCACAATACGACCCTTTTGGGAAAGAGGGGGGTAAGTTGGCGCT‑AAGccccggggg > 1:10701/1‑145 (MQ=255) aaaGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAAGACA‑ATTAAAAAAGGGCCGGCCACGCGGAGTTT‑CTCCCTTAC‑GGCTTTTCCCTCCCATTACCACCCTTTGCGTCACGAGGGGGGTACGTTGGCGCTCAAACCCgggggg > 1:852371/1‑146 (MQ=255) aGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTg < 1:564593/151‑1 (MQ=255) cGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTggg < 1:372177/150‑1 (MQ=255) cgcgGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTGGGGCGCAATCCTc < 1:785419/151‑1 (MQ=255) gcgGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTGGGGCGCAATCCTCa < 1:432051/151‑1 (MQ=255) gcgggcgAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAgg < 1:267835/151‑1 (MQ=255) gcgggcgAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAgg < 1:214975/151‑1 (MQ=255) gcgAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTAtt > 2:486349/1‑151 (MQ=255) cgAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTAtt > 2:735357/1‑150 (MQ=255) gAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTa > 2:758898/1‑150 (MQ=255) gAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTa < 2:8513/150‑1 (MQ=255) aGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTaa > 1:60495/1‑150 (MQ=255) aGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTaa > 1:16466/1‑150 (MQ=255) | CATATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACA‑TTGAACAGCGGACCGGCAACGCGGATTCT‑CTCCGTTAC‑GGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCT‑AAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAA > NC_000913/1153423‑1153718 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |