Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,073,654 | T→C | 13.8% | intergenic (‑148/+476) | mdtA ← / ← PP_3586 | multidrug transporter membrane fusion protein/transposase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,073,654 | 0 | T | C | 13.8% | 67.2 / 3.8 | 29 | intergenic (‑148/+476) | mdtA/PP_3586 | multidrug transporter membrane fusion protein/transposase |
Reads supporting (aligned to +/- strand): ref base T (12/13); new base C (2/2); total (14/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01 |
TTGAACGAACCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGC > NC_002947/4073510‑4073793 | ttGAACGAACCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTcc < 2:351044/150‑1 (MQ=255) aacgaacCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGAACTGCGATCGGTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGTCGGAATTTTAAATGAAAACGGCCCGGAGCGTGTTCCggg > 2:42359/1‑150 (MQ=255) gaacCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCg > 1:419971/1‑150 (MQ=255) ttGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACa > 2:93508/1‑150 (MQ=255) gTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACACGTTCCGGGCCGttttcat < 1:345698/150‑6 (MQ=255) gcGGAATGACGACGTGAGGAGGATGAACTATATGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGgt < 2:192210/149‑1 (MQ=255) ggAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGgtgt > 2:504003/1‑150 (MQ=255) ttGAGAAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCg > 2:274233/1‑150 (MQ=255) aggaTGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCg < 2:419971/150‑1 (MQ=255) ggaTGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCgg < 2:162156/150‑1 (MQ=255) aaCGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTg > 1:220127/1‑150 (MQ=255) tAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACACGTTCCGGGc > 2:335552/1‑109 (MQ=255) tAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACACGTTCCGGGc < 1:335552/109‑1 (MQ=255) cATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCt < 2:366042/150‑1 (MQ=255) cAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCt > 2:415576/1‑150 (MQ=255) ggCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTcgc < 1:274233/150‑1 (MQ=255) cAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAAcc > 2:322597/1‑150 (MQ=255) aCCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTAt < 2:498049/150‑1 (MQ=255) ccTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATcaaacaa > 1:233282/1‑150 (MQ=255) ccTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATcaaacaa > 2:47372/1‑150 (MQ=255) tGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACaaaa < 2:220127/150‑1 (MQ=255) aaTTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCTGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGGTTGGTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGttctttc < 1:415576/150‑1 (MQ=38) aTTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGttct > 1:177816/1‑146 (MQ=255) aTTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGttct < 2:177816/146‑1 (MQ=255) ttttAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTa < 2:233282/150‑1 (MQ=255) ttttAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTa < 2:16572/150‑1 (MQ=255) tGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGTCTGGTATGGGCCTTTCGGGGGTGAACCCGCTTCTATCAAACAAAAGATAAACCACTGTTCTTTCTAGCGaaaa > 2:259093/1‑150 (MQ=38) aCGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGACCCCCCACTGTTCTTTCTAGCGAAAAAAACa < 1:93508/150‑1 (MQ=255) gCCCGGAACACGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGgc > 2:337554/1‑150 (MQ=255) | TTGAACGAACCTTGTTCGGAATGACGACTTGAGGAGGATGAACGATAAGCACTGCCATCGTTTTAGCAAAGGGCCTTTACCAGAGGTTTACCTTTGCCTGACGTTGCCAAGAGGCGGAATTTTAAATGAAAACGGCCCGGAACGTGTTCCGGGCCGTTACAGGGTGTTGCTTGGAGCGGTGGCGGGCCTGAAGGTTCGGTTGCCTGGTCTGGCCCTTTCGCGGGTGAACCCGCTTCTATCAAACAAAAGATAACCCACTGTTCTTTCTAGCGAAAAAAACAGGC > NC_002947/4073510‑4073793 |
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: ‑ |