Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,269,872:1 | +C | 16.7% | intergenic (‑528/‑33) | PP_3740 ← / → mrdA‑I | MFS transporter/transpeptidase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,269,872 | 1 | . | C | 16.7% | 66.8 / 12.3 | 24 | intergenic (‑528/‑33) | PP_3740/mrdA‑I | MFS transporter/transpeptidase |
Reads supporting (aligned to +/- strand): ref base . (12/8); new base C (2/2); total (14/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-01 |
ATGACTCCGCAATGAACAAATCTTCATATTCCCCCTCAGCGCCCTGATCATGGGGCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGCCTGTGTGCCGTCTTGGTGGG > NC_002947/4269732‑4270014 | atGACTCCGCAATGAACAAATCTTCATATTCCCCCTCAGCGCCCTGATCATGGGGCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATTGc > 1:226423/1‑150 (MQ=255) ccGAAATGAACAAATCTGCATATTCCCCCCCTCCGCCCTGATCATGGGGCACACAGCAGATGCCAGGCAGCGGGCAGCCGCGGCTCGTCTTGACGGTGCGGTACGGTTTGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATGGCCAcccc < 2:129416/150‑1 (MQ=255) aaTCTTCATATTCCCCCTCAGCGCCCTGATCATGGGGCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATAGCCACCCCTCgggtgcgact < 2:596642/150‑8 (MQ=255) tatTCCCCCTCAGCGCCCTGATCATGGGGCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTCACAGATg > 1:209795/1‑150 (MQ=255) atTCCCCCTCAGCGCCCTGATCATGGGGCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTcaca > 1:378344/1‑145 (MQ=255) atTCCCCCTCAGCGCCCTGATCATGGGGCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGAGCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTcaca < 2:378344/145‑1 (MQ=255) gCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACg > 1:140576/1‑150 (MQ=255) cagccagcGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCAcc < 2:356406/150‑1 (MQ=255) agccagcGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAAc < 1:226269/113‑1 (MQ=255) agccagcGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAAc > 2:226269/1‑113 (MQ=255) agccagcGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCt > 2:367965/1‑150 (MQ=255) gccagcGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTg < 2:209795/150‑1 (MQ=255) cGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCCCGAAAAGGAAATGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATg < 2:226423/150‑1 (MQ=255) atcctgaggggtggcaataCGCG‑CGCGATCGGGCCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTaaa > 1:65847‑M1/20‑90 (MQ=255) atcctgaggggtggcaataCGCG‑CGCGATCGGGCCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTaaa < 2:65847‑M1/71‑1 (MQ=255) gggTTGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGAACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTc > 2:69031/1‑150 (MQ=255) ttGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTCTTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGc > 2:227907/1‑150 (MQ=255) ttGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGc > 2:405089/1‑150 (MQ=255) gTAGCATCCGCGGCCCGATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGCCtgt > 1:212214/1‑150 (MQ=255) aataCGCG‑CGCGATCGGGCCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACgaaaa > 1:83640‑M1/5‑87 (MQ=255) aataCGCG‑CGCGATCGGGCCGCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACgaaaa < 2:83640‑M1/83‑1 (MQ=255) aTCCGCGGCCCGATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGCCTGTGTGcc > 2:279154/1‑150 (MQ=255) GGCCCGATCGCGC‑GCGTATTGCCACCCCTGAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGCCTGTGTGCCGTCTTg > 2:91736‑M1/1‑150 (MQ=255) gATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGCCTGTGTGCCGTCTTGGTggg < 1:405089/150‑1 (MQ=255) | ATGACTCCGCAATGAACAAATCTTCATATTCCCCCTCAGCGCCCTGATCATGGGGCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGCCTGTGTGCCGTCTTGGTGGG > NC_002947/4269732‑4270014 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |