Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,066,146 | G→C | 9.7% | Q165E (CAA→GAA) | PP_4463 ← | carbon‑nitrogen hydrolase family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,066,146 | 0 | G | C | 9.7% | 124.5 / 4.7 | 42 | Q165E (CAA→GAA) | PP_4463 | carbon‑nitrogen hydrolase family protein |
Reads supporting (aligned to +/- strand): ref base G (18/19); new base C (2/2); total (20/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01 |
GTTTGACCTGGAGCAAGGAAATAGCACTGCGTTTCAATTGCCCGAGCGCGGCAGAGTACATCCCAATGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTGGCGCTACTGCACTGGATTCGCCATAGCGCATACCGTCAGGAGTGACGATGTCGAA > NC_002947/5066007‑5066287 | gTTTGACCTGGAGCAAGGAAATAGCACTGCGTTTCAATTGCCCGAGCGCGGCAGAGTACATCCCAATGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTcc > 1:750872/1‑150 (MQ=255) cTGGAGCAAGGAAATAGCACTGCGTTTCAATTGCCCGAGCGCGGCAGAGTACATCCCAATGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGc < 2:380854/150‑1 (MQ=255) cTGGAGCAAGGAAATAGCACTGCGTTTCAATTGCCCGAGCGCGGCAGAGTACATCCCAATGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGc > 1:380854/1‑150 (MQ=255) cTGGAGCAAGGAAATAGCACTGCGTTTCAATTGCCCGAGCGCGGCAGAGTACATCCCAATGATCTTTACCAGTTTGAAGGGTAAATGCAGCAGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGc > 2:150927/1‑150 (MQ=255) gAGCAAGGAAATAGCACTGCGTTTCAATTGCCCGAGCGCGGCAGAGTACATCCCAATGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGaa > 1:572580/1‑150 (MQ=255) gAAATAGCACTGCGTTTCAATTGCCCGAGCGCGGCAGAGTACATCCCAATGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTa < 2:1073431/150‑1 (MQ=255) gTTTCAATTGCCCGAGCGCGGCAGAGTACATCCCAATGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAg < 2:750872/150‑1 (MQ=255) gTTTCAATTGCCCAAGCGCGGCAGAGTACATCCCAATGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACAGCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAg < 1:56692/150‑1 (MQ=255) ttCAATTGCCCGAGCGCGGCAGAGTACATCCCAATGATCTTTACCAGTTTGAAGGGTAAAGGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGcc < 2:331890/150‑1 (MQ=255) gCCCGAGCGCGGCAGAGTACATCCCAATGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTa < 2:502939/128‑1 (MQ=255) gCCCGAGCGCGGCAGAGTACATCCCAATGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTa > 1:502939/1‑128 (MQ=255) cccGAGCGCGGCAGAGTACATCCCAATGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTAGTAGCAGATCGCAAAGCCATACTTaa < 1:176144/150‑1 (MQ=255) gAGCGCGGCATAGTACATCCCAATGATCTTTACCAGTTTTAAGGGTAAATGCCGCGGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGc < 1:303922/150‑1 (MQ=255) agagTACATCCCAATGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGc < 1:1026474/124‑1 (MQ=255) agagTACATCCCAATGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGc > 2:1026474/1‑124 (MQ=255) tCCCAATGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGaccacc > 1:740440/1‑150 (MQ=255) cccAATGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCg > 1:463612/1‑150 (MQ=255) cAATGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGAt > 1:361219/1‑150 (MQ=255) aaggATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACTCAACCAGCTTTTTGAATAGCTCCTGGAAGCGAATGTCGTAGCAGCTCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATa < 2:89913/147‑1 (MQ=255) aaTGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATa > 1:185525/1‑150 (MQ=255) ccAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAATCGAATGTCGTAGCAGATCGCACAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCt < 2:572580/150‑1 (MQ=255) aaGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTg > 1:936547/2‑150 (MQ=255) tAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACAGTTATCTCTGTACCTGGCGCTACTGCAc < 1:150927/150‑1 (MQ=255) cAATTACATCAGCACCAAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTGGCGCTACTGCACTGGATTCGCCATAGCGCa < 2:626654/150‑1 (MQ=255) cATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTGGCGCTACTGCACTGGATTCGCCATAGCGCATACCGt < 2:936547/150‑1 (MQ=255) aTCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTg < 1:1086971/113‑1 (MQ=255) aTCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTg > 2:1086971/1‑113 (MQ=255) aTCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTGGCGCTACTGCACTGGATTCGCCATAGCGCATACCGTc > 1:287866/1‑150 (MQ=255) aTCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTGGCGCTACTGCACTGGATTCGCCATAGCGCATACCGTc > 1:822526/1‑150 (MQ=255) tCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTGGCGCTACTGCACTGGATTc < 2:873425/131‑1 (MQ=255) tCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTGGCGCTACTGCACTGGATTc > 1:873425/1‑131 (MQ=255) ccAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTGGCGCTACTGCACTGGATTCGCCATAGCGCATACCGTCAGGAGTGa < 2:463612/150‑1 (MQ=255) aaCGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAAt > 2:444500/1‑38 (MQ=255) aaCGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAAt < 1:444500/38‑1 (MQ=255) aaCGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTGGCGCTACTGCACTGGATTCGCCATAGCGCATACCGTCAGGAGTGACg > 2:772121/1‑150 (MQ=255) ccggaGCTTTTCCAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTGGCGCTACTGCACTGGATTCGCCATAGCGCATACCGTCAGGAGTGAc < 1:104953/140‑1 (MQ=255) ccggaGCTTTTCCAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTGGCGCTACTGCACTGGATTCGCCATAGCGCATACCGTCAGGAGTGAc > 2:104953/5‑144 (MQ=255) aaCCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTGGCGCTACTGCACTGGATTCGCCATAGCGCATACCGTCAGGAGTGACGATGTc < 2:361219/150‑1 (MQ=255) aaCCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTGGCGCTACTGAACTGGATTCGCCATAGCGCATACCATCAGGCGTGACGATGTc > 2:219337/1‑150 (MQ=255) cggaGCTTTTCCAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGaccac > 1:996502/4‑80 (MQ=255) cggaGCTTTTCCAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGaccac < 2:996502/77‑1 (MQ=255) cAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTGGCGCTACTGCACTGGATTCGCCATAGCGCATACCGTCAGGAGTGACGATGTCGaa < 2:940462/150‑1 (MQ=255) | GTTTGACCTGGAGCAAGGAAATAGCACTGCGTTTCAATTGCCCGAGCGCGGCAGAGTACATCCCAATGATCTTTACCAGTTTGAAGGGTAAATGCCGCCGGCAGTACAATTACATCAGCACCCAACGCAACCAGCTTTTGGAAAAGCTCCGGGAAGCGAATGTCGTAGCAGATCGCAAAGCCATACTTAAGGCCTTCAATATCGACCACCGATACCTCTGTACCTGGCGCTACTGCACTGGATTCGCCATAGCGCATACCGTCAGGAGTGACGATGTCGAA > NC_002947/5066007‑5066287 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |