Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 922,314:1 | +C | 10.9% | intergenic (+21/+52) | macB → / ← cspD | macrolide ABC transporter peremase/ATPase/inhibitor of DNA replication, cold shock protein homolog |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 922,314 | 1 | . | C | 10.9% | 246.0 / 29.0 | 64 | intergenic (+21/+52) | macB/cspD | macrolide ABC transporter peremase/ATPase/inhibitor of DNA replication, cold shock protein homolog |
Reads supporting (aligned to +/- strand): ref base . (28/29); new base C (3/4); total (31/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 = 8.71e-01 |
TGGTTTTTCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCT > NC_000913/922167‑922462 | tGGTTTTTCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑c < 3:4249/149‑1 (MQ=255) tttttCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGa < 3:96058/149‑1 (MQ=255) tttCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATc < 2:1050001/149‑1 (MQ=255) aCCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGc > 2:330999/1‑149 (MQ=255) ccGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCt < 2:1680728/149‑1 (MQ=255) cGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATCggctg < 4:1328882/149‑4 (MQ=255) ctgctgctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTg < 1:330999/149‑1 (MQ=255) gCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCAtttt > 1:1690764/1‑129 (MQ=255) gCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCAtttt < 2:1690764/129‑1 (MQ=255) tCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATga > 3:245222/1‑149 (MQ=255) cGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaaga < 4:298802/147‑1 (MQ=255) cGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaaga > 3:298802/1‑147 (MQ=255) gggATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCt > 2:884859/1‑149 (MQ=255) tttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACaa > 3:474041/1‑129 (MQ=255) tttATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACaa < 4:474041/129‑1 (MQ=255) gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGc < 1:884859/149‑1 (MQ=255) gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGc < 4:245222/149‑1 (MQ=255) gTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGc < 1:642230/93‑1 (MQ=255) gTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGc > 2:642230/1‑93 (MQ=255) ttACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGc < 4:1518736/144‑1 (MQ=255) ttACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGc > 3:1518736/1‑144 (MQ=255) gCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGt < 1:265339/104‑1 (MQ=255) gCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGt > 2:265339/1‑104 (MQ=255) gCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGcc > 2:1131186/1‑136 (MQ=255) gCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGcc < 1:1131186/136‑1 (MQ=255) cACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTa > 2:1244897/1‑97 (MQ=255) cACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTa < 1:1244897/97‑1 (MQ=255) gAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCAttt > 3:820907/1‑86 (MQ=255) gAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCAttt < 4:820907/86‑1 (MQ=255) gAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaaga < 3:1753538/118‑1 (MQ=255) gAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACagaaga > 4:1753538/1‑118 (MQ=255) gAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACggg > 4:1259461/1‑149 (MQ=255) aTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcac < 3:1259461/149‑1 (MQ=255) gCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGt < 3:1436218/92‑1 (MQ=255) gCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGt > 4:1436218/1‑92 (MQ=255) tGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACaa < 3:1556180/139‑1 (MQ=255) tGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACaa > 4:1556180/1‑139 (MQ=255) cAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTc < 4:646369/116‑1 (MQ=255) cAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTc > 3:646369/1‑116 (MQ=255) ctctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg > 4:635642/1‑148 (MQ=255) ctctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg < 3:635642/148‑1 (MQ=255) tctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATga > 3:332401/1‑78 (MQ=255) tctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATga < 4:332401/78‑1 (MQ=255) gagaGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggtg > 2:980440/1‑149 (MQ=255) gTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATga < 3:758038/66‑1 (MQ=255) gTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATga > 4:758038/1‑66 (MQ=255) tAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACtt < 1:1116938/98‑1 (MQ=255) tAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACtt > 2:1116938/1‑98 (MQ=255) tAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCtttt > 4:276478/1‑133 (MQ=255) tAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCtttt < 3:276478/133‑1 (MQ=255) ttttGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCaaa < 1:980440/149‑1 (MQ=255) aaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg > 4:1015747/1‑118 (MQ=255) aaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg < 3:1015747/118‑1 (MQ=255) tGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGAt > 2:1288977/1‑107 (MQ=255) tGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGAt < 1:1288977/107‑1 (MQ=255) gCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAAc < 1:521062/135‑1 (MQ=255) gCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAAc > 2:521062/1‑135 (MQ=255) cAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAg > 4:655131/1‑149 (MQ=255) aG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACGCAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGc > 4:1598495/1‑149 (MQ=255) aG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggt < 4:1371010/121‑1 (MQ=255) aG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggt > 3:1371010/1‑121 (MQ=255) aG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCaaa > 3:1313865/1‑131 (MQ=255) aG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCaaa < 4:1313865/131‑1 (MQ=255) aG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTc > 1:1413038/1‑145 (MQ=255) aG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTc < 2:1413038/145‑1 (MQ=255) g‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACt > 4:564853/1‑132 (MQ=255) g‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACt < 3:564853/132‑1 (MQ=255) g‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCt > 2:1764487/1‑149 (MQ=255) | TGGTTTTTCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCT > NC_000913/922167‑922462 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |