Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 6,104,110 | T→A | 23.6% | intergenic (+32/‑82) | PP_5353 → / → PP_5354 | hypothetical protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 6,104,110 | 0 | T | A | 23.6% | 60.2 / 19.1 | 34 | intergenic (+32/‑82) | PP_5353/PP_5354 | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (11/15); new base A (4/4); total (15/19) | |||||||||||
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.50e-01 |
GAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGAGCT > NC_002947/6103967‑6104253 | gaagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg > 1:101792/1‑150 (MQ=255) gaagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCgg > 1:594623/1‑150 (MQ=255) aagaaTTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAAAACCCCGGGCAAGCCCggg < 2:427420/150‑1 (MQ=255) agaaTTCGCGCTTCGACAGAGGTGTAGTATTTAATAAAATCTTGCATATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCAAGCCCgggg < 2:223697/150‑1 (MQ=255) gcTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgctt > 1:222310/1‑150 (MQ=255) ttCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTttttgctttt < 2:75154/150‑1 (MQ=255) ccAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGTCGAGAAGCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCaa > 2:87516/1‑150 (MQ=255) gggCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGcc < 1:132673/145‑1 (MQ=255) gggCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGcc > 2:132673/1‑145 (MQ=255) gAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGc < 1:551522/150‑1 (MQ=255) tGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGgcgc > 1:30852/1‑150 (MQ=255) ggCACATGGTTCGCGCCCAAGGCGAGACACAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCa < 2:101792/150‑1 (MQ=255) caTGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGTTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGt > 1:497628/1‑150 (MQ=255) aTGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGtt < 2:594623/150‑1 (MQ=255) ggCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCaga < 1:89407/150‑1 (MQ=255) cAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAg < 1:341368/135‑1 (MQ=255) cAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAg > 2:341368/1‑135 (MQ=255) aGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAAc < 1:609730/150‑1 (MQ=255) ttGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTg > 1:448672/1‑150 (MQ=255) ttGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGCGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTg < 2:135493/150‑1 (MQ=255) ccGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGAt < 2:448672/150‑1 (MQ=255) tCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTc > 1:250632/1‑150 (MQ=255) aCAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGa < 1:74313/150‑1 (MQ=255) gCTAATTGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGc < 2:497628/150‑1 (MQ=255) aaCACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGgc < 2:22485/66‑1 (MQ=255) aaCACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGgc > 1:22485/1‑66 (MQ=255) aaaaaCCCCGGGCAAGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGccc > 2:419516/4‑150 (MQ=255) ccccGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGa < 2:250632/150‑1 (MQ=255) ccccGGGCAAGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGaa > 2:221352/1‑126 (MQ=255) ccccGGGCAAGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGaa < 1:221352/126‑1 (MQ=255) cccGGGCAAGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGaa > 2:1464/1‑108 (MQ=255) cccGGGCAAGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGaa < 1:1464/108‑1 (MQ=255) cccGGGCAAGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGAg > 2:570049/1‑150 (MQ=255) cGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGAGCt < 1:87516/150‑1 (MQ=255) | GAAGAATTCGCGCTTCGACCAGGGCGAAGTGTTCAGCGAAAGCTGGCACATGGTTCGCGCCCAAGGCGAGAACCAGCCTTGGCTGGTAGCCGGTATCCGTCAAAACAGCTAATCGCTGCGTGCTAAACACAGAACCCCGGGCTTGCCCGGGGTTTTTGCTTTTGCCAAAGTGGGCTACTAGGGTATAAGGCGCAGCGTTGTCATCTAGGTCAGAGGAAGTGAACCGTGGAAGAAGTGATCGAACAACTCCGTGAAGCGAATGAGCCGGTACCGGTGCCCCTGGAGCT > NC_002947/6103967‑6104253 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |