Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,912,914 | G→A | 16.2% | E43K (GAG→AAG) | spuI → | glutamylpolyamine synthetase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,912,914 | 0 | G | A | 16.2% | 76.8 / 8.1 | 37 | E43K (GAG→AAG) | spuI | glutamylpolyamine synthetase |
Reads supporting (aligned to +/- strand): ref base G (16/15); new base A (4/2); total (20/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.67e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01 |
AAAGGTGTTGCATCATGTCGGTACCCCCGCGTGCCGTTCAGCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATCGGCGATGCTGACCGCATCTGCTATCCAATCCCC > NC_002947/5912774‑5913060 | aaaGGTGTTGCATCATGTCGGTACCCCCGCGTGCCGTTCAGCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCa > 2:212860/1‑150 (MQ=255) aGGTGTTGCATCATGTCGGTACCCCCGCGTGCCGTTCAGCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGc > 2:750201/1‑150 (MQ=255) ggTGTTGCATCATGTCGGTACCCCCGCGTGCCGTTCAGCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGcc < 1:212860/150‑1 (MQ=255) tCGGTACCCCCGCGTGCCGTTCAGCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTAc > 2:699426/1‑150 (MQ=255) tCGGTACCCCCGCGTGCCGTTCAGCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAgg < 2:583154/145‑1 (MQ=255) tCGGTACCCCCGCGTGCCGTTCAGCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAgg > 1:583154/1‑145 (MQ=255) ccGCGTGCCGTTCAGCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGc > 1:688287/1‑150 (MQ=255) cTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCAcc < 1:252431/108‑1 (MQ=255) cTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCAcc > 2:252431/1‑108 (MQ=255) tcaggatgTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTg > 1:90497/8‑150 (MQ=255) cGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGAGTTGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTg < 2:116835/150‑1 (MQ=255) gcgTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTg < 1:699426/150‑1 (MQ=255) gTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCt > 2:665353/1‑147 (MQ=255) gTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCt < 1:665353/147‑1 (MQ=255) gTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGa > 1:252445/1‑150 (MQ=255) aaCATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGtt < 1:750201/150‑1 (MQ=255) gAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTc > 2:88664/1‑150 (MQ=255) aGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCg > 2:474034/1‑150 (MQ=255) gTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGaa > 1:288066/1‑150 (MQ=255) tctACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGTGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCa > 2:292779/1‑150 (MQ=255) gTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCggg > 2:683822/1‑150 (MQ=255) ttGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGc < 1:474034/150‑1 (MQ=255) aCCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTgg > 2:522977/1‑150 (MQ=255) ttCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGcc > 1:729456/1‑150 (MQ=255) ttCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGcc > 2:569334/1‑150 (MQ=255) aTTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCAACAGGACCCACAAGGTTTACGTGGAAGGTATGGATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGAc > 2:313640/1‑150 (MQ=255) tGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACAt < 1:569334/150‑1 (MQ=255) gCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATc < 1:326830/150‑1 (MQ=255) gATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATCGGc < 2:729456/150‑1 (MQ=255) atatGAATGGTGTGGTGCGTGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATCGGCg > 1:327401/1‑150 (MQ=255) tggtgtggtgCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCCGGGCCTGGACATCGGCGATGCTGa < 2:90497/150‑1 (MQ=255) tggtgCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGCCCTTTACTGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATCGGCGATGCTGACCGCa < 2:755867/150‑1 (MQ=255) tGGCAAGCGCATAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATCGGCGATGCTGACCGCATCTGCTa < 1:683822/150‑1 (MQ=255) cgcATAGAGCGCACCATCCTCCACACGGATTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTTCACCGTCGAAAGCACCGGGCTGGGCCTGGAGATCGGCGATGCTGACCGCATCTGCTATCCAATc < 1:771285/150‑1 (MQ=255) cATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATCGGCGATGCTGACCGCATCTGCTATCCAATccc < 2:288066/150‑1 (MQ=255) aTAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATCGGCGATGCTGACCGCATCTGCTATCCAATcccc < 1:88664/150‑1 (MQ=255) aTAAAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACCTCAACGGTTCCACCGTCGAAAGCACCGGACTGGGCCTGGACATCGGCGATGCTGACCGCGTCTGCTGTCCACTcccc > 2:92397/1‑150 (MQ=255) | AAAGGTGTTGCATCATGTCGGTACCCCCGCGTGCCGTTCAGCTTAACGAAGCGAACGCGTTCCTTAAGGAACATCCTGAGGTTCTCTACGTTGACCTTCTGATTGCAGATATGAATGGTGTGGTGCGTGGCAAGCGCATAGAGCGCACCAGCCTCCACAAGGTTTACGAGAAAGGCATCAATCTGCCGGCCTCCCTCTTTGCCCTGGACATCAACGGTTCCACCGTCGAAAGCACCGGGCTGGGCCTGGACATCGGCGATGCTGACCGCATCTGCTATCCAATCCCC > NC_002947/5912774‑5913060 |
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: ‑ |