Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,555,432 | Δ1 bp | 10.2% | intergenic (+28/+46) | sbcB → / ← PP_1366 | exodeoxyribonuclease I/transcriptional regulator MvaT |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,555,432 | 0 | G | . | 10.2% | 171.4 / 15.3 | 59 | intergenic (+28/+46) | sbcB/PP_1366 | exodeoxyribonuclease I/transcriptional regulator MvaT |
Reads supporting (aligned to +/- strand): ref base G (23/29); new base . (3/3); total (26/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01 |
GCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTT > NC_002947/1555286‑1555569 | gCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGTCTGATGTCATATAAACAAAAATGCCGGCA‑ATGTGc > 2:682270/1‑150 (MQ=255) gCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGc < 2:1099372/150‑1 (MQ=255) gagttCGAAAAGGCGCGCCAGCAAGCCTGGGAGGCCGCCGACGAGGGGGGCAGGGGCGTGTGGGAGGCCTGGCAACTGCATGCGTGACAACTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAGCGCCGGCA‑ATGTGCCg < 2:447144/147‑1 (MQ=255) gACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCg > 2:1087164/1‑150 (MQ=255) tCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACAT‑TGCCGGCGt < 1:544943/150‑1 (MQ=255) cGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGtt < 1:864137/150‑1 (MQ=255) gAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGACGCAGGCAGGCGCGTGTTGGAGGCCTGGCACCTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGttt > 1:312706/1‑150 (MQ=255) gcgcgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTttgttg < 1:1087164/150‑1 (MQ=255) gcgcCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTtgtg > 1:131484/1‑150 (MQ=255) cAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGaa < 2:131484/150‑1 (MQ=255) cAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGCAGGCCTGGCAACTGCATGCTAGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGaa > 2:635749/1‑150 (MQ=255) cAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGtttt < 2:1322231/136‑1 (MQ=255) cAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGtttt > 1:1322231/1‑136 (MQ=255) cTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAgg > 1:1212069/1‑150 (MQ=255) ggCGCAGACGAGGGAGACCGCCGGGTGCGTGAGGGCTGGCAAGTGCATGCTGGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGAGAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATg < 2:874200/150‑1 (MQ=255) gcgcCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGc < 2:53775/150‑1 (MQ=255) cgaGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAAt < 2:1791727/150‑1 (MQ=255) aggcaggcaggcGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGACGACG‑ACGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCa < 1:1572696/150‑1 (MQ=255) gcgcGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACTAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGCATGCGGATAATCAGTCCAGCAgg < 2:950855/150‑1 (MQ=255) gcgcGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGCATGCGGATAAGTAGTACAGCAgg > 1:390534/1‑150 (MQ=255) gtgtTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTg > 1:1904968/1‑150 (MQ=255) gtTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGcc < 2:855037/150‑1 (MQ=255) ggAGGCCCGGCAACTGCACGCTCAACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAg < 2:136792/150‑1 (MQ=255) ggCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACAT‑TGCCGGCGTTTttgt > 2:144586/1‑95 (MQ=38) ggCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCACAT‑TGCCGGCGTTTttgt < 1:144586/95‑1 (MQ=38) gCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCttt < 1:1080423/150‑1 (MQ=255) ccTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTc < 1:871942/150‑1 (MQ=255) cTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCa > 1:1722388/1‑150 (MQ=255) cTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCa > 2:91292/1‑150 (MQ=255) cTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCAAGCTGTCa > 1:83699/1‑150 (MQ=255) gCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTcg < 2:1212069/150‑1 (MQ=255) gCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTcgccgc < 2:867858/150‑1 (MQ=255) tCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGa > 1:1159582/1‑104 (MQ=255) tCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGa < 2:1159582/104‑1 (MQ=255) tCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGccc > 1:1330494/1‑150 (MQ=255) aCAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAACCGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCaa < 1:682270/150‑1 (MQ=255) ttGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTATTTCAATGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTgg < 1:1861131/150‑1 (MQ=255) aaCGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCa > 1:1620353/1‑150 (MQ=255) aaCGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCa > 2:1490320/1‑150 (MQ=255) cAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAg < 1:466557/150‑1 (MQ=255) aaTCGGCTGATGGCAGATAAACAAAAACGCCGCCA‑ATCTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGGGAATGCTGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAgg < 1:423004/150‑1 (MQ=255) tCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAgg < 1:1241765/94‑1 (MQ=255) tCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAgg > 2:1241765/1‑94 (MQ=255) gCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGc > 1:1390309/1‑72 (MQ=255) gCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGc < 2:1390309/72‑1 (MQ=255) gCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCgg > 2:1039701/1‑67 (MQ=255) gCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCgg < 1:1039701/67‑1 (MQ=255) gCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTg > 1:710276/1‑150 (MQ=255) cAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgc < 1:199515/150‑1 (MQ=255) aTAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTgccgc < 2:1904968/150‑1 (MQ=255) acaaCAAAAACGCCGGCACAT‑TGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTg > 2:436257/3‑115 (MQ=38) acaaCAAAAACGCCGGCACAT‑TGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTg < 1:436257/113‑1 (MQ=38) aCAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCttt < 1:1160579/150‑1 (MQ=255) aCAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGGTGCACAGCTATCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCttt < 2:312706/150‑1 (MQ=255) aaaaCGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGt > 2:1148225/1‑150 (MQ=255) aaaCGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGtt > 1:1273876/1‑150 (MQ=255) aaCGCCGGCACAT‑TGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGttt > 1:1802136/1‑150 (MQ=255) gCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTc > 2:1733472/1‑110 (MQ=255) gCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTc < 1:1733472/110‑1 (MQ=255) | GCGACTTCGAAAAGGCGCGCCAGCAAGCCTGGGAGGGCGCCGACGAGGCAGGCAGGCGCGTGTTGGAGGCCTGGCAACTGCATGCTCGACAATTGCAGGAACGCTATTCAATCGGCTGATGGCAGATAAACAAAAACGCCGGCA‑ATGTGCCGGCGTTTTTGTTGTGCGAAGCAAGAGGGCGTATGCGGATAATCAGTCCAGCAGGGTTGCCCAGCTTTCAACCACGTCGCCGCCCCAAGTGGCTTTCCACTCTTTCAGGGTCTTGTGGTTGCCGCCTTTGGTTT > NC_002947/1555286‑1555569 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |