Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,003,504 | A→G | 9.7% | D414G (GAT→GGT) | PP_4409 → | site‑specific recombinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,003,504 | 0 | A | G | 9.7% | 108.4 / 3.1 | 41 | D414G (GAT→GGT) | PP_4409 | site‑specific recombinase |
Reads supporting (aligned to +/- strand): ref base A (21/16); new base G (2/2); total (23/18) | |||||||||||
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.14e-01 |
GCAACATCTCGCGTCGAAGCTTGGACTACATAGCCTTCGCCATAGCAACATTTTCTTTCTAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGCCTGG > NC_002947/5003370‑5003644 | gCAACATCTCGCGTCGAAGCTTGGACTACATAGCCTTCGCCATAGCAACATTTTCTTTCTAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGa < 1:1263063/135‑1 (MQ=255) gCAACATCTCGCGTCGAAGCTTGGACTACATAGCCTTCGCCATAGCAACATTTTCTTTCTAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGa > 2:1263063/1‑135 (MQ=255) gCAACATCTCGCGTCGAAGCTTGGACTACATAGCCTTCGCCATAGCAACATTTTCTTTCTAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCa < 1:502327/150‑1 (MQ=255) cAACATCTCGCGTCGAAGCTTGGACTACATAGCCTTCGCCATAGCAACATTTTCTTTCTAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCtt > 2:236611/1‑147 (MQ=255) cAACATCTCGCGTCGAAGCTTGGACTACATAGCCTTCGCCATAGCAACATTTTCTTTCTAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCtt < 1:236611/147‑1 (MQ=255) aTCTCGCGTCGAAGCTTGGACTACATAGCCTTCGCCATAGCAACATTTTCTTTCTAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAAt > 1:481098/1‑150 (MQ=255) ctcGCGTCGAAGCTTGGACTACATAGCCTTCGCCATAGCAACATTTTCTTTCTAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATAc < 2:1401753/150‑1 (MQ=255) tcGCGTCGAAGCTTGGACTACATAGCCTTCGCCATAGCAACATTTTCTTTCTAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGaaa > 1:407979/1‑146 (MQ=255) tcGCGTCGAAGCTTGGACTACATAGCCTTCGCCATAGCAACATTTTCTTTCTAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGaaa < 2:407979/146‑1 (MQ=255) ttGGACTACATAGCCGCCGGCATAGCAACATTTTTTTTCTAAAGAACTACCTAGAACACAGTCAAGGGCAACGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAACTACGCGAAGGTCGAcc > 2:476011/1‑150 (MQ=255) aGCCTTCGCCATAGCAACATTTTCTTTCTAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGGTATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTAttg > 1:78699/1‑150 (MQ=255) aGCCTTCGCCATAGCAACATTTTCTTTCTAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGGTATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTAttg > 1:78741/1‑150 (MQ=255) ccTTCGCCATAGCAACATTTTCTTTCTAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTAttgtt < 1:1632114/150‑1 (MQ=255) ttCGCCATAGCAACATTTTCTTTCTAAAGAACTACCTAGAGCACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGATGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTgg > 1:1572320/1‑150 (MQ=255) aaCATTTTCTTTCTAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGcc > 2:1115341/1‑150 (MQ=255) aCATTTTCTTTCTAAAGAACTACCTAGAACACAGACAAGGACAAGGTTTGACAGATATCGAATTGCTATTGATTACATGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCt > 2:546738/1‑150 (MQ=255) ttttCTTTCTAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACg < 2:481098/150‑1 (MQ=255) tAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACgg < 1:1495702/142‑1 (MQ=255) tAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACgg > 2:1495702/1‑142 (MQ=255) tAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTccc > 1:1157375/1‑150 (MQ=255) tAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTccc > 1:912576/1‑150 (MQ=255) aaaGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCt < 2:1157375/150‑1 (MQ=255) aGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTaa < 1:728754/150‑1 (MQ=255) cTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAAt < 2:220110/150‑1 (MQ=255) tAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATg > 1:164699/1‑150 (MQ=255) aGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTAttgtt < 2:134536/104‑1 (MQ=255) aGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTAttgtt > 1:134536/1‑104 (MQ=255) aaGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGATATATTAAAt > 1:1304514/1‑150 (MQ=255) aaGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCACTAAGAAAAAATGGAAATATTAAAt > 2:92028/1‑150 (MQ=255) ggCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGaaaaaa < 1:819783/133‑1 (MQ=255) ggCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGaaaaaa > 2:819783/1‑133 (MQ=255) ggCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCaa > 1:433081/1‑150 (MQ=255) cAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGc > 2:38778/1‑150 (MQ=255) tttGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGa < 2:433081/150‑1 (MQ=255) tttGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGa > 1:170573/1‑150 (MQ=255) ttGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGGTATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGaa < 2:78741/150‑1 (MQ=255) ttGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGGTATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGaa < 2:78699/150‑1 (MQ=255) gCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTGAAATTCAGtgctt < 1:92028/150‑5 (MQ=255) ttGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCa > 2:690234/1‑140 (MQ=255) ttGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTgaa < 1:1115341/150‑1 (MQ=255) ttGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTAGTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCa < 1:690234/140‑1 (MQ=255) cTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGcc > 2:987671/1‑150 (MQ=255) gCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGCCTgg > 2:685126/1‑150 (MQ=255) | GCAACATCTCGCGTCGAAGCTTGGACTACATAGCCTTCGCCATAGCAACATTTTCTTTCTAAAGAACTACCTAGAACACAGTCAAGGGCAAGGTTTGACTGATAGCGAATTGCTATTGATTACTGGCCACGCTGATATCAGAAGTCTTCAGAAATACGCGAAGGTCGACCGAGAATTATTGTTGGAAAAGATCGCCTACGGAAACTCCCTAAGAAAAAATGGAAATATTAAATCAAGCACTGAATTTCAAATTCAGTACCTTGAAGAACGCCTGG > NC_002947/5003370‑5003644 |
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: ‑ |