Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913_3_bme_pgi | 526,800 | A→G | 29.8% | E1180E (GAA→GAG) | rhsD → | Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913_3_bme_pgi | 526,800 | 0 | A | G | 29.8% | 71.0 / 37.6 | 57 | E1180E (GAA→GAG) | rhsD | Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor |
Reads supporting (aligned to +/- strand): ref base A (20/20); new base G (6/11); total (26/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.89e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01 |
CGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATA > NC_000913_3_bme_pgi/526652‑526945 | cgcAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGCATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAAACGGCAATACGGATTGGAGCGCGGaa > 1:162774/1‑149 (MQ=255) gggCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatg < 2:270547/147‑1 (MQ=255) tGCGGGTTGAAGGTGGAGGAACGGGCAAGACAGGTGGAGCCGGAATACACACCGGCGAGAAAAGCTCATCTTTATCACTGGGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTatga < 2:83873/149‑1 (MQ=11) tGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACGGCCGCGGGCGCTTCTCAGCGAAACGGGCATAACGGGTTGGAGCCGGGAatataa > 2:255127/1‑147 (MQ=14) gCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgat < 2:63484/149‑1 (MQ=255) cGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatgatg > 1:750747/1‑149 (MQ=255) ggCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatga < 1:347667/145‑1 (MQ=255) ggCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATatga > 2:347667/1‑145 (MQ=255) tGACGGGGGAGCAACTGGCCAGACAGGTGGATCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACTGCAATACGGCGTGGAGCGCGGAATATGATGAATgg < 2:68847/149‑1 (MQ=255) gagggTGGAGCAACGCGGCAGACAGGTGGAGCCGGAACACACCCCGGCGCGAAAAGCTCATCTTTATCACTGCGACCATCGGGGACCGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATggg < 2:483918/146‑1 (MQ=17) gACGGTGGAGGAATTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCTAAAAGCTTATCTTTATCACTGCGACCACCGGGGAATGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATggg < 1:268073/149‑1 (MQ=17) gACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATggg < 1:853129/149‑1 (MQ=255) gTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTCTGATGAATGGGGCaa > 2:868156/1‑149 (MQ=255) tGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCCTCTTTATCCCTGCGACCACCGGTGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAAc < 1:521225/149‑1 (MQ=11) ggtgcTGGCCAGAGAGGTGCGGTCGGAATACACACCGGCGCGAAAAGCTCATCTTTATGACTGCGACCACCGGGGACTGTCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGc < 2:705009/145‑1 (MQ=17) cTGGCCAGACAGGTGGATCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAACACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTaa < 2:846404/149‑1 (MQ=255) cTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTACGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTaa < 2:287867/149‑1 (MQ=11) gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCATAAtaa > 1:150651/1‑147 (MQ=255) gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATga < 1:595708/149‑1 (MQ=255) gCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACAGCAATACGGCGTGGACCGGGGAGTATGATGAATGGGCAAACCAACTAAATga > 1:163917/1‑149 (MQ=14) cctgagaGGGGGAGGCGGAATACACACCGGCGTGAAGAGCTTATCTTTATCACCGCGACGACCGGGGACTGCCGCTGGGGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGTAATATGATGAATGGGGCAACCAGCTTAATgag < 2:608802/143‑1 (MQ=17) cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATgagg > 1:723199/1‑149 (MQ=255) cAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCGGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAACGTAATaaga > 2:240550/1‑148 (MQ=255) aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCCCTGCGACGACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACTGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga < 2:435164/149‑1 (MQ=11) aGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga < 2:295171/149‑1 (MQ=255) aGACAAGGGGAGCCGGAATAAACACTGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATgagga < 2:697270/149‑1 (MQ=11) gAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGGTTAATGAGGAGAACCCGca > 2:874316/1‑149 (MQ=255) cGGAATACACACATGCGCTAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcac < 1:419226/149‑2 (MQ=11) ggAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGGAAGAAGAATGGGGCAACAAGCTTAATGAGGAGAACACGcatcacc > 1:698239/1‑147 (MQ=11) aaTACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacctg < 1:163770/149‑5 (MQ=1) tACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacctgca > 1:298536/1‑143 (MQ=11) tACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGATCGCAGAATATGATGAATGGTGCAACCAGCTTAATGAGGAGAACCAGCATAATgtgtg > 1:125892/1‑148 (MQ=17) tACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTACTGAGGAGAACACGcatcatcttta > 2:771511/1‑144 (MQ=255) tACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTAGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacctgca > 1:553659/1‑143 (MQ=11) acacacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGTGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGcatcacctgcac < 1:9584/149‑8 (MQ=11) acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCa > 1:849134/1‑149 (MQ=255) acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCa > 1:721572/1‑149 (MQ=255) acacCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCGTATTAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCa > 1:849453/1‑149 (MQ=255) gcgcGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAAAACCCGCATAATGTGTATCAGCCGTa > 1:700936/1‑149 (MQ=255) aaaaCCTCATCTTTATCACTGCGACCACCGTGGACTGCCGCTGGCGCTTATCAGCGGAGACGGCACTCCGGCGTGGACCGCGTAATAGGATGAATGAGGAAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAACCGTACAGTc > 1:481336/1‑149 (MQ=18) tcctcattATCACTGCGACCACCGGGGACCGCCGCGGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACTAGCTTAATGTGGATAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAg < 2:324208/143‑1 (MQ=17) cTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGCGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgca > 1:792532/1‑149 (MQ=255) gagcACTGCGAGCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGGGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGgcagca < 2:642352/146‑1 (MQ=255) tGCGACCACCGGCGACCGCCGCTGGCGCTTATCAGCGAAGACGGACATATGGCGTGGAGCGCTGAATATGATGAATGGGGCAACCAGCTTAATGAGGCGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCatga < 2:223829/149‑1 (MQ=17) cGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATc > 1:327127/1‑149 (MQ=255) ggggACTGGCGCTGGGCCTTATCAGCGAAGAGCGGAATACGATGTGGAGCGAGAAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGTCAGCAGCATGATGAGGAATCa < 2:199820/149‑1 (MQ=17) ggggACTGCCGCTGGCGCTTATCAGCGAAGATGGCAATACTGTGTGGGGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCa < 1:30713/149‑1 (MQ=255) ttATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATAAGGGCTGTACTATAACCGt > 1:586698/1‑149 (MQ=255) aGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTCCTATAACCGTCACCg > 2:108601/1‑149 (MQ=255) gCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCgg > 1:50709/1‑149 (MQ=255) gAAGACGGAAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAAACAGCTTAATGAGGAGCACCCGCTTAATGTGTATCAGCCGTACCGTCTGCCAGGGCTGCAGCATGATGAGGAATCATGGCTGTACTATAACCGTCACCGGta < 2:176813/149‑1 (MQ=17) tACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGGACCCGCATCCTGTGGATCAGCCGTACCGTCTGCCAGGGCAGCAGCCAGATGAGGAATCGGGGCTGTACTATAACCGCAACCGGTCCTACGATCCGt > 2:453811/1‑149 (MQ=14) actgCGTGGAGCTCCGTATATGTTGAATGGGGCAACAAGCTTAATGAGGCGGGCGCTTAGGATGTGTATCAGCCGTATCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGtt < 1:205052/146‑1 (MQ=17) gCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCa > 1:59580/1‑149 (MQ=255) cGTAGAGCGCGGAATATGATGAATGGGGCAACCAGCATAAGGAAGCGTTCACGGATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAg < 1:568227/149‑1 (MQ=17) tGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCAGCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAggg > 1:707519/1‑149 (MQ=255) gcgcGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCgg < 1:574638/149‑1 (MQ=255) cgGAATATGATGAATGGGGCAACTAGTTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGtat < 1:407243/149‑1 (MQ=255) ggAATATGATTAATGGGGTAACCAGCTTAATGAGGAGAACCTGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGtata < 1:394686/149‑1 (MQ=255) | CGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGAGCGCGGAATATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCATGTGTATCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAATCAGGGCTGTACTATAACCGTCACCGGTACTACGATCCGTTGCAGGGGCGGTATA > NC_000913_3_bme_pgi/526652‑526945 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |