Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,228,383 | G→A | 13.3% | intergenic (‑127/+409) | PP_5624 ← / ← PP_3707 | hypothetical protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,228,383 | 0 | G | A | 13.3% | 67.8 / 3.3 | 30 | intergenic (‑127/+409) | PP_5624/PP_3707 | hypothetical protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (14/12); new base A (2/2); total (16/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-02 |
CCCACATGTTGAGCTCAGATGAATAGGAGCCAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCACTACGACATTTCCGT > NC_002947/4228250‑4228517 | cccacATGTTGAGCTCAGATGAATAGGAGCCAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATg > 3:272191/1‑149 (MQ=255) ccacATGTTGAGCTCAGATGAATAGGAGCCAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGc > 3:173426/1‑149 (MQ=255) aTGTTGAGCTCAGATGAATAGGAGCCAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGAc > 3:296417/1‑148 (MQ=255) tGTTGAGCTCAGATGAATAGGAGCCAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACaa < 5:161477/149‑1 (MQ=255) tCAGATGAATAGGAGCCAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTcga > 5:100647/1‑149 (MQ=255) aGATGAATAGGAGCCAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTcgac > 2:5963/1‑148 (MQ=255) gAGCCAGGGATACTTTCCTCCCCTGTCCTTGGAGTGGCGAAGCTAGAAAGGTACGGTGATTTGCCCAACGATTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATCCCATTAg < 7:121068/149‑1 (MQ=255) ccAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGc > 4:215664/1‑148 (MQ=255) gATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCAc < 4:173426/149‑1 (MQ=255) ccTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGa < 6:100647/149‑1 (MQ=255) tCCTTGGAGTTGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGGCATAACCTTGTTAAATGGAGTAt < 2:268031/149‑1 (MQ=255) tCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTAt < 4:296417/149‑1 (MQ=255) gCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGAtt > 3:176126/1‑149 (MQ=255) gCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGAtt > 1:202477/1‑149 (MQ=255) gCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGAt > 1:160592/1‑148 (MQ=255) aGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCAGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGcc > 6:136835/1‑149 (MQ=255) aGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCAGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGcc > 6:141291/1‑149 (MQ=255) ggTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCtt < 4:205423/106‑1 (MQ=255) ggTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTc < 4:272191/149‑1 (MQ=255) ggTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATAACCtt > 3:205423/1‑106 (MQ=255) cccAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGt < 3:215664/148‑1 (MQ=255) ttttCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGa > 7:42601/1‑148 (MQ=255) gTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGc < 1:5963/146‑1 (MQ=255) gTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCAc > 2:278678/1‑149 (MQ=255) cTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCACTa < 1:278678/148‑1 (MQ=255) ccAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCACTACGa > 1:258656/1‑148 (MQ=255) ccAAGCTGGCGTCTCCTATGCCAGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCACTACGa < 5:136835/148‑1 (MQ=255) ccAAGCTGGCGTCTCCTATGCCAGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCACTACGa < 5:141291/148‑1 (MQ=255) gCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCACTACGACAttt > 5:297400/1‑149 (MQ=255) gCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCACTACGACATTTCCGt < 6:297400/149‑1 (MQ=255) | CCCACATGTTGAGCTCAGATGAATAGGAGCCAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCACTACGACATTTCCGT > NC_002947/4228250‑4228517 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |