Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,913,204 | G→A | 10.8% | intergenic (‑186/‑96) | PP_4317 ← / → PP_4318 | hypothetical protein/transposase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,913,204 | 0 | G | A | 10.8% | 95.0 / 3.7 | 37 | intergenic (‑186/‑96) | PP_4317/PP_4318 | hypothetical protein/transposase |
Reads supporting (aligned to +/- strand): ref base G (18/15); new base A (2/2); total (20/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01 |
GAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTG > NC_002947/4913079‑4913312 | gAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTat > 1:64606/1‑135 (MQ=21) gAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAg > 1:147554/1‑136 (MQ=35) aaGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTat > 1:15779/2‑135 (MQ=18) ggggACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAAt > 1:277099/1‑136 (MQ=40) ggggACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAAt > 1:514752/1‑136 (MQ=40) gACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCa > 1:222673/1‑136 (MQ=38) cGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAg > 2:506176/1‑136 (MQ=37) cGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCAAGTCAg > 2:55137/1‑136 (MQ=37) gAGCGCAGCGCTGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGa < 2:358098/136‑1 (MQ=37) gcgcagcgcaGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAg > 1:457104/1‑136 (MQ=37) gcagcgcagGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGAc > 1:379505/1‑136 (MQ=37) agcgcagGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGAccc > 1:374651/1‑136 (MQ=34) gcagGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAg > 2:491821/1‑136 (MQ=34) ggggCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGa < 2:346570/136‑1 (MQ=34) ggCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGAtg > 2:361681/1‑136 (MQ=34) ccGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGaa > 1:536658/1‑108 (MQ=35) ccGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGaa < 2:536658/108‑1 (MQ=35) ccGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGAtgtg < 1:498962/136‑1 (MQ=34) ccGGATGATGGGAGCGCAGCGTTTTTTGGTTACTCTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGAtgtg > 1:383134/1‑136 (MQ=21) ggATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTggg < 2:141389/136‑1 (MQ=34) gatgatGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGATCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTggg < 1:442724/135‑1 (MQ=21) gatgatGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGATCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTggg < 1:440208/135‑1 (MQ=21) gatgatGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGATCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTggg > 2:442724/1‑135 (MQ=21) gatgatGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGATCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTggg > 2:440208/1‑135 (MQ=21) gCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGcgcg > 1:40795/1‑136 (MQ=32) gCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGcgcg > 1:386975/1‑136 (MQ=32) cGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGa < 2:188878/136‑1 (MQ=32) cGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGa > 1:401757/1‑136 (MQ=32) tttttGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTc < 1:361681/136‑1 (MQ=32) aCTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGaaaaaa < 2:200203/136‑1 (MQ=32) aCTTTTTGTCGCGGTTGACAAAAAGTGACTGGCCTTAAGGGCCATACCTTTCTTTGTCTTCTTTAGATTTCGCTGAACTCCAGTCAGGTTGACCGGAGCTATGTGGGCTGATGTACTAGCGCGATTCGATaaaaaa < 2:249109/136‑1 (MQ=11) aaGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCa < 1:140980/136‑1 (MQ=32) tgTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCaa < 2:194838/135‑1 (MQ=17) aGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCaa < 2:109063/136‑1 (MQ=32) gACTCGCCGTAATGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACt < 2:207269/136‑1 (MQ=17) cTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTg < 2:350545/136‑1 (MQ=31) cTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTg < 2:420440/136‑1 (MQ=31) | GAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTG > NC_002947/4913079‑4913312 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |