| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 4,269,861 | Δ1 bp | 16.7% | intergenic (‑517/‑44) | 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,861 | 0 | G | . | 16.7% | 61.5 / 9.9 | 24 | intergenic (‑517/‑44) | PP_3740/mrdA‑I | MFS transporter/transpeptidase |
| Reads supporting (aligned to +/- strand): ref base G (12/8); new base . (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 = 4.42e-01 | |||||||||||
ACCCCTGACTCGTCATCATGACTCCGCAATGAACAAATCTTCATATTCCCCCTCAGCGCCCTGATCATGGGGCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGCCTGTGTGCCGTCTTG > NC_002947/4269715‑4270009 | accCCTGACTCGTCATCATGACTCCGCAATGAACAAATCTTCATATTCCCCCTCAGCGCCCTGATCATGGGGCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGccc < 1:70536/150‑1 (MQ=255) 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) | ACCCCTGACTCGTCATCATGACTCCGCAATGAACAAATCTTCATATTCCCCCTCAGCGCCCTGATCATGGGGCACACAGCCGATGCCAGCCAGCGGGCAGCCGGGGCTCGTCTTGAGGGTGCGGTACGGGTTGCGTAGCATCCGCGGCCCGATCGCGC‑GCGTATTGCCACCCCTCAGGATCGATTCACAGATGCCGCAACCCATCCCGCTAAAGGACCACGAAAAGGAAAAGCACCTGGTCAACCGTCGGCTTCTGGCCTGCGCTGCCATGGTGTTCAGCCTGTGTGCCGTCTTG > NC_002947/4269715‑4270009 |
| 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: ‑ |