Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,231,788 | C→A | 10.0% | intergenic (‑4/+80) | PP_3709 ← / ← PP_5625 | hypothetical protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,231,788 | 0 | C | A | 10.0% | 115.0 / 3.6 | 40 | intergenic (‑4/+80) | PP_3709/PP_5625 | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base C (17/19); new base A (2/2); total (19/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01 |
ACTACCAAGGTAGGCAGTAATGGGCTTGTCGAAGCCCTGATGTTTCGCTTTGGCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGTAGAGGGTCTGCAGTCTACTCTGCACGGATCATCGCCGATAGGAAGATAGGCCTTTGGTTGGT > NC_002947/4231665‑4231915 | aCTACCAAGGTAGGCAGTAATGGGCTTGTCGAAGCCCTGATGTTTCGCTTTGGCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCt > 2:126159/1‑150 (MQ=255) tACCAAGGTAGGCAGTAATGGGCTTGTCGAAGCCCTGATGTTTCGCTTTGGCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTa < 2:18041/125‑1 (MQ=255) tACCAAGGTAGGCAGTAATGGGCTTGTCGAAGCCCTGATGTTTCGCTTTGGCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTa > 1:18041/1‑125 (MQ=255) ccAAGGTAGGCAGTAATGGGCTTGTCGAAGCCCTGATGTTTCGCTTTGGCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGc > 1:285935/1‑150 (MQ=255) cAAGGTAGGCAGTAATGGGCTTGTCGAGGCCCTGATGTTTCGCTTTGGCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCt > 1:240608/1‑150 (MQ=255) cAGTAATGGGCTTGTCGAAGCCCTGATGTTTCGCTTTGGCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACAcc < 2:285935/150‑1 (MQ=255) aTGGGCTTGTCGAAGCCCTGATGTTTCGCTTTGGCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAAt > 1:369021/1‑150 (MQ=255) gTCGAAGCCCTGATGTTTCGCTTTGGCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACct < 1:217867/150‑1 (MQ=255) aGCCCTGATGTTTCGCTTTGGCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTgag < 1:186295/150‑1 (MQ=255) ccgATGTTTCGCTTTGGCAAGCCGCTTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGAcgc < 2:332955/148‑1 (MQ=255) ttCGCTTTGGCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAAt < 1:524783/95‑1 (MQ=255) ttCGCTTTGGCAATCCGCCTCAGATCTAGAACGATTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAAt > 2:524783/1‑95 (MQ=255) cTTTGGCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTcc > 1:483628/1‑105 (MQ=255) cTTTGGCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTcc < 2:483628/105‑1 (MQ=255) cTTTGGCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGTAGAGGGTc < 1:126159/150‑1 (MQ=255) tGGCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAGACTCAGTTCAACCGGTTCCTGCATTTTACTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGTAGAGGGTCTGc > 2:807531/1‑150 (MQ=255) tGGCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTACTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGTAGAGGGTCTGc > 2:807353/1‑150 (MQ=255) ggCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGAACAAGACCGTCGCGTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGTAGAGGGTCTGCa < 2:543043/150‑1 (MQ=255) gCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGTAGAGGGTCTGCAGTCTActc < 2:369021/150‑1 (MQ=255) tCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGTAGAGGGTCTGCAGTCTACTCTGc > 1:89626/1‑150 (MQ=255) tCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGGAGAGGGTCTGCAGTCTaagatcg > 1:137081/1‑144 (MQ=255) gAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATa < 2:785233/99‑1 (MQ=255) gAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATa > 1:785233/1‑99 (MQ=255) ttAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGTAGAGGGTCTGCAGTCTACTCTGCACGGATCATCGCCGa > 1:827574/1‑150 (MQ=255) aGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCATACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGTAGAGGGTCTGCAGTCTACTCTGCACGGATCATCGCCGATAg < 2:240608/150‑1 (MQ=255) gCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGc < 2:160366/67‑1 (MQ=255) gCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGc > 1:160366/1‑67 (MQ=255) ctctTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGTAGAGGGTCTGCAGTCTACTCTGCACGGATCATCGCCGATAGGa < 2:827574/150‑1 (MQ=255) ttGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGTAGAGGGTCTGCAGTCTACTCTGCACGGATCATCGCCGATAGGAAGa < 2:89626/150‑1 (MQ=255) ccGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGTAGAGGGTCTGCAGTCTACTCTGCACGGATCATCGCCGATAGGa < 1:110119/144‑1 (MQ=255) ccGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGTAGAGGGTCTGCAGTCTACTCTGCACGGATCATCGCCGATAGGa > 2:110119/1‑144 (MQ=255) aGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGac < 2:91342/73‑1 (MQ=255) aGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGac > 1:91342/1‑73 (MQ=255) aGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACctct < 2:541943/79‑1 (MQ=255) aGTTCAACCGGTTCCTGCATTTTCATACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACctct > 1:541943/1‑79 (MQ=255) aGTTCAACCGGTTCCTGCATTTTACTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGTAGAGGGTCTGCAGTCTACTCTGCACGGATCATCGCCGATAGGAAGATAGGCCTttggttgg < 1:807531/150‑1 (MQ=255) aGTTCAACCGGTTCCTGCATTTTACTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGTAGAGGGTCTGCAGTCTACTCTGCACGGATCATCGCCGATAGGAAGATAGGCCTttggttgg < 1:807353/150‑1 (MQ=255) gTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGTAGAGGGTCTGCAGTCTACTCTGAAGGGATCATCGCCGATAGGAAGATAGGCCTTtggttggt > 1:225038/1‑150 (MQ=255) gTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTagag < 2:121873/49‑1 (MQ=255) gTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTagag > 1:121873/1‑49 (MQ=255) | ACTACCAAGGTAGGCAGTAATGGGCTTGTCGAAGCCCTGATGTTTCGCTTTGGCAATCCGCCTCAGATCTAGAACGCTTAAGCTCTTGCCGTGAAAACTCAGTTCAACCGGTTCCTGCATTTTCCTACTCCATGTAATCGCTAGACGCCTCCGCTAGAGACACCGAAATAGACACCTCTGAGACGCGGTAGAGGGTCTGCAGTCTACTCTGCACGGATCATCGCCGATAGGAAGATAGGCCTTTGGTTGGT > NC_002947/4231665‑4231915 |
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: ‑ |