Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,385,155 | G→T | 16.6% | intergenic (‑119/+42) | proS ← / ← oprD | proline‑‑tRNA ligase/basic amino acid specific porin OprD |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,385,155 | 0 | G | T | 16.6% | 53.5 / 2.3 | 24 | intergenic (‑119/+42) | proS/oprD | proline‑‑tRNA ligase/basic amino acid specific porin OprD |
Reads supporting (aligned to +/- strand): ref base G (9/11); new base T (2/2); total (11/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01 |
AAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCA > NC_002947/1385017‑1385297 | aaaTATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTgccgggct > 2:91244‑M2/1‑142 (MQ=255) atTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTgccgggcttttt < 2:266419‑M2/150‑13 (MQ=255) aCTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAgg > 1:720533/1‑150 (MQ=255) ccGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGcc > 2:99017/1‑150 (MQ=255) ccTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGa > 1:202152/1‑150 (MQ=255) aaGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACgccgggctttttcttggctgacaagtgtcac > 2:95397‑M2/1‑101 (MQ=255) aaGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACgccgggctttttcttggctgacaagtgtcac < 1:95397‑M2/132‑32 (MQ=255) aaaTAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGt < 1:99017/150‑1 (MQ=255) gggggggAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAg > 1:273175/1‑150 (MQ=255) ggAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTtcgt > 2:143218/1‑150 (MQ=255) gAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTtcgtc < 1:266419/150‑1 (MQ=255) aCGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACgccgggctttttcttggctgacaagtgtc < 1:346861‑M2/93‑30 (MQ=255) aCGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCAATCACgccgggctttttcttggctgacaagtgtc > 2:346861‑M2/1‑64 (MQ=255) cGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGtcgatc < 1:5931/150‑1 (MQ=255) gCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATgg < 2:202152/150‑1 (MQ=255) cAGGATGCCCCAGGGCGTGACACTTGTCCGCCAGGAAAAAGCCCGGTGTGGTTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCGCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAAt < 2:278251/150‑1 (MQ=255) aTGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCgta < 2:720533/150‑1 (MQ=255) gCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCgtagt < 1:143218/150‑1 (MQ=255) ttGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCa < 1:841612/106‑1 (MQ=255) ttGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCa > 2:841612/1‑106 (MQ=255) aaaaaGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGcc < 2:273175/150‑1 (MQ=255) aGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGac > 2:647615/1‑150 (MQ=255) aGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGac < 1:647615/150‑1 (MQ=255) ccGGCGTGATTGCCGGGCTTTTTCACGTGCGCCCGACTAAAACAGAAAGTGATGGGATGCTCAGCGGGTACTGGACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACCACATCGTAAGAACGGTAGGTGGCCTGacgca > 1:574932/1‑150 (MQ=255) gATTgccgggctttttcttggctgacaagtgtcacgccctggggcatcctgtacgccttcgtttcccccccgtacaatgcggggctatttcaacgtcttaggcaatcgacacggctaaccatgc > 2:393002‑M2/1‑4 (MQ=255) gATTgccgggctttttcttggctgacaagtgtcacgccctggggcatcctgtacgccttcgtttcccccccgtacaatgcggggctatttcaacgtcttaggcaatcgacacggctaaccatgc < 1:393002‑M2/124‑121 (MQ=255) | AAATATTGACTGGTGCGCATGGTTAGCCGTGTCGATTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCA > NC_002947/1385017‑1385297 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |