Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 266,349 | 0 | C | A | 56.0% | 30.4 / 84.4 | 50 | V69F (GTT→TTT) | yafY | lipoprotein, inner membrane; degP regulator; CP4‑6 prophage |
Reads supporting (aligned to +/- strand): ref base C (4/18); new base A (27/1); total (31/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.79e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.82e-05 |
CACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAT > NC_000913/266209‑266499 | cACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAgggg > 2:871899/1‑149 (MQ=255) cctcctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTggg < 2:587143/151‑1 (MQ=255) cctcctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTgg < 2:118430/150‑1 (MQ=255) ctcctGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGTAGTAAAAACTGCCAgggggg > 1:738153/1‑150 (MQ=255) ggTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGa < 1:898716/151‑1 (MQ=255) ggTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAAGTAAAAACTGCCAGGGGGGGCCgg > 1:282930/1‑150 (MQ=255) aTTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAAACCGCAAGGGGGGGCCgggaa > 2:661445/1‑147 (MQ=255) tCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAAACCGCAAGGGGGGGCCGGGTaa > 1:60002/1‑149 (MQ=255) cAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCAAGGGGGGGCCgggaaa > 1:642272/1‑143 (MQ=255) tGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAGGGGGGGCCGGATAAGTAGTCAA‑cc > 1:765942/1‑149 (MQ=255) gTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCtgt < 2:901701/151‑1 (MQ=255) aTCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAACCCGTAAGGGGGGGCCCCAGAAAGAGTCAA‑CCagt > 1:640080/1‑148 (MQ=255) cAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTAAGTAAAAACTGCAAGGTGGGACCGGATAAATAGTCAA‑ACCGTGACCAtt > 1:709216/1‑150 (MQ=255) aaTACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCAAGGGGGGGGCGGGGAAATAGGCAA‑ACTTTGACCCtt > 1:277404/1‑149 (MQ=255) aTACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGGAAAAACCGCAAAGGGGGGCCGGGAAAAAAGTCAA‑ACTGTGAACcttctt > 1:497391/1‑151 (MQ=255) acacAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAAGGGGGGCCGGATAAAAGGGCAA‑CCCGTGAGCcttctt > 2:952383/1‑149 (MQ=255) acaGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACCTAAAAACTGCAAAGGGGGACCGGGTAAGGAGTCAA‑CCTGTGACCAttctt > 1:941243/1‑147 (MQ=255) caGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTg < 2:427958/151‑1 (MQ=255) aGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGa < 2:228675/151‑1 (MQ=255) aGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCACCAAGTGGGGCCGGGAAAGGAGTCAA‑ACTCTGACCCttcttctccg > 1:503464/1‑146 (MQ=255) gcATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGAcat < 2:228902/151‑1 (MQ=255) catcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAAGGGGGGCCGGGTAAAGAAGCAA‑CCTGTGACCATTCTTTCCCGACATCATcc > 1:895246/1‑151 (MQ=255) catcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACTTCAAAGTGGGACCCGCTAAGGAGTGAA‑ACTCGGACCCCTCTTCGCTgccggca > 1:800777/1‑142 (MQ=255) tGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACCGCAAGGGGGGGCCGGATAAGTAGGCAA‑CCCGTGGCCAATCTTCGCTGACATCATCCttgtgg > 2:214949/1‑151 (MQ=255) tGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCCAGGGGGGGCCGGGGAAATAGTCAA‑CCCGTGACCCTTCTTCGCCGACATCATCCtggtgg > 2:82366/1‑151 (MQ=255) tCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAACCTGCCAGGGGGGACCGGGGAAATAGGCAA‑ACTGTGACCCTTCTTCGCTGACATCCTCCTCGTGGCTGGCCACAt > 1:869434/1‑145 (MQ=255) tGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTATAAACTACAAAGGGGGGCCGGGAAAGAAGTCAA‑CCTGTGAACATTCTTTGCTGACATCATCCTGGTGGGTGGCCACAACAATGaaa > 2:668636/1‑150 (MQ=255) gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCGCAAAGGGGGAGCGGGGAAAGAGTCAA‑CCCGTGACCATTCTTCGCTGGCAAAATCCCGGTGGGGGACCACATCACCGaaa > 2:103998/1‑149 (MQ=255) atatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAAcgcg < 2:642082/151‑1 (MQ=255) tatCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACCGCAAGGTGGGGCCGGGGAAGTAGTCAACCCTGT‑ACCCTTTTTCGCCTACACCCTCCTGGTGGCCGACCACCACACCTAAACGCgg > 1:80526/1‑151 (MQ=255) atCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCgg < 1:952180/150‑1 (MQ=255) atCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGt < 2:632281/151‑1 (MQ=255) ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAACCTGCAAAGGGGGACCCGATAAGTAGTCAA‑ACCGTGACCATTCTTCGCCGACATCATCCCGGTGGCTGACCACATCACTGAAACGCGgggg > 1:186327/1‑151 (MQ=255) cACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTggg < 1:845811/151‑1 (MQ=255) aCGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAGGGGGGACCGGGTAAGTAGTCAA‑CCTGTGAGCATTCTTTGCTGACATCATCCTGGTGGGTGGCCACCTCACTGAAACGCGgggggg > 1:620145/1‑151 (MQ=255) ggTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGa < 2:582264/150‑1 (MQ=255) agtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGTAAAAACTGCCAGGGGGGGCCGGATAAGTAGGCAA‑CCTGTGACCATTCTTTGCTGACATCATCCTGGGGGGTGACCACATCACTGAAAACCGggggggggggg > 1:65818/1‑144 (MQ=255) tagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCCGCCAGGGGGGGCCGGATAAGAAGGCAA‑CCCGTGACCCTTCTTCGCTGACATCATCCTGGGGGCTGACCACATCACTGAAACGCGggggggggggggcc > 1:115992/1‑142 (MQ=255) aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGCAAAAACCGCAAGGAGGGGCCGGAAAAGGAGGCAA‑ACCGTGGCCCTTCTTTGCCGACCTCAACCCCGTGGCTGGGCCCATCACTGAAAAGCGggggggggggggccct > 2:373045/1‑138 (MQ=255) gCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAAtt < 1:93895/149‑1 (MQ=255) aTCCTCGCAGCCATCGGGCGACTGGCGGGTATGTAAAAACTGCGAGGGGGGGGCGGATAAATAGGCAA‑ACCGTGGCCCTTCTTCTCCGTCATCATCCCGGGGGCTTCCCACAACACCGGAACCGGGGGGGGGGGGGGCCTTtaa > 1:14270/1‑144 (MQ=255) cccgcgGCCATCGGGCGCCTGGCCGGTACGTAAAGACTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACGATTCTTCGCGGACATCATCCTGATGTCTGCCCACAGCACTGAAGCGGGGTGGGGTGTCGTCCTTATAtttc < 2:937917/143‑2 (MQ=255) cTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAAGTGGGGCCGGATAAATAAGCAA‑ACTGTGACCATTCTTCGCTGACATCATCCTGGGGGCTGGCCACACCACTGAAAAGCGGTGGGGGGAGGGCCTTAAATTTTCTg > 2:452678/1‑148 (MQ=18) cGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGc < 2:59989/151‑1 (MQ=255) cAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAGGGGGGGCCGGAGAAATAGTCAA‑CCTGTGACCCTTCTTCGCTGACATCATCCTGGTGGCGGAACCCATCACTGAAACGCGGGGGGGGGGGGGCCTTAAAATTACTGAttggc > 2:343919/1‑146 (MQ=18) ccATCGGGCGACTGGCGGGTACGTAAAAACTGCAAGGGGGGGCCGGGTAAGTAGTCAA‑CCTGTGACCATTTTTCGCTGACATCATCCTGGGGGGTGACCACATCACTGAAACGCGGGGGGGGGGGGGCCTTTAATTTACTGATGGCCtt > 1:325227/1‑149 (MQ=255) aTCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAgg < 1:587088/151‑1 (MQ=255) gCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtcatc < 2:577647/151‑1 (MQ=255) ggTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCa < 2:123408/151‑1 (MQ=255) cGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGc < 2:277332/151‑1 (MQ=255) ccTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAt > 1:541324/1‑151 (MQ=255) | CACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAA‑CCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAT > NC_000913/266209‑266499 |
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: ‑ |