Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,189,042 | C→T | 12.4% | intergenic (+26/+16) | PP_1038 → / ← tatC‑I | hypothetical protein/Sec‑independent protein translocase protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,189,042 | 0 | C | T | 12.4% | 109.5 / 4.3 | 45 | intergenic (+26/+16) | PP_1038/tatC‑I | hypothetical protein/Sec‑independent protein translocase protein |
Reads supporting (aligned to +/- strand): ref base C (21/18); new base T (2/4); total (23/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.14e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTC > NC_002947/1188900‑1189163 | tGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCcgcg > 2:346791/1‑150 (MQ=255) gAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGACACCCGTACGAGACGCCAGCCTATGAAGTATGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGaa < 2:689949/150‑1 (MQ=255) gTGGAGCTGGTGGTGGCTGACGAGATGATTGGGCAGTTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑Aggg < 2:498619/150‑1 (MQ=255) gCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGcc > 2:659656/1‑150 (MQ=255) tGACGAGCTGTTTGTGCAGCTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTTTGGCGGCTCGCCGAGTACTAGAACGCATCACGGCTCCTCGCGGGCATTACCGCGAA‑AGGGCCCTACAATCCAGCcac < 1:381289/150‑1 (MQ=255) tGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGaa‑a < 2:274308/130‑1 (MQ=255) tGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGaa‑a > 1:274308/1‑130 (MQ=255) tGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCcac > 2:526341/1‑150 (MQ=255) gAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAAtc > 1:308401/1‑150 (MQ=255) tGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCt > 1:136567/1‑150 (MQ=255) gCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGc < 1:21634/150‑1 (MQ=255) tGCGGTCAAGCAGAGCCACGCGCACGAGACGCCAGCCTATGAAGTCTGGCGGTTCGCCGAGTTCTAGAACGCATCACTGCTATTCGGGGGCACGCCCGTGAA‑AGGGCCCTACATTCCAGCCACAATCTCCTTGGCTGGCTCCGGCCGCtt < 1:564831/150‑1 (MQ=255) tGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCAAGCA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCtt > 1:250647/1‑150 (MQ=255) gcgcTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTc < 2:75619/150‑1 (MQ=255) cTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGa < 1:712786/97‑1 (MQ=255) cTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCGTGCCCGCGa > 2:712786/1‑97 (MQ=255) cTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAAc > 2:131610/1‑150 (MQ=255) tCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACa > 1:714409/1‑150 (MQ=255) aGCAGAGCCACCCGTACGAGACGCCAGCATATGAAGTCGGGCTGCTCGCCGAGTTCTAGAACGCATCACGCCAGTTTGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTc < 1:52277/150‑1 (MQ=255) agCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACaa < 2:136567/150‑1 (MQ=255) aCCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCGCGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCa < 1:131610/150‑1 (MQ=255) cccGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAg > 1:437412/1‑150 (MQ=255) gagaCGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAACACACGTCTGAACTCCAG‑TACAATTGCTGTTCCTGGCTCCGGCCGCTTCAACATTACACAAGCCAGCAAcccc > 1:93598/1‑150 (MQ=255) gCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCACtctc > 1:134401/1‑150 (MQ=255) aGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCCCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGaa < 2:308401/150‑1 (MQ=255) cTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACac > 2:208385/1‑149 (MQ=255) aTGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCa < 2:714409/150‑1 (MQ=255) gTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAAc < 2:116969/129‑1 (MQ=255) gTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAAc > 1:116969/1‑129 (MQ=255) gcggcTCGCCGAGTTCTAGAACGCATCACGGGTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCgg < 1:100627/150‑1 (MQ=255) gcggcTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCgg > 1:684250/1‑150 (MQ=255) gcggcTCGCCGAGTTCTAGAACGCATCACGGCTATTCACGCGCACGCCCGAGAA‑AAGACACAACAATCCAACCACAATCTCCTTTCCTGTCTCCGGCCTCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCgg > 2:290981/1‑150 (MQ=255) gcggcTCGCCGAGTCCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGGCCTAGAATCCGGCCACAATCTGCTTGCCTGGGTCCGGCCGCTGCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCgg > 2:291160/1‑150 (MQ=255) cGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGc < 1:346791/150‑1 (MQ=255) gCCGAGTTCGAGACCGCATCACGGTTTTTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCTCTTCAACTGCTGACCAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGcc < 1:15522/150‑1 (MQ=255) ttCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACt < 2:742475/141‑1 (MQ=255) ttCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACt > 1:742475/1‑141 (MQ=255) ttCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATc < 1:659656/150‑1 (MQ=255) ttCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATc > 2:509818/1‑150 (MQ=255) cTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGt < 1:526341/150‑1 (MQ=255) tAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTc > 1:555929/1‑150 (MQ=255) aCGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCAccc < 2:783743/116‑1 (MQ=255) aCGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCAccc > 1:783743/1‑116 (MQ=255) gccctttCGCGGGCGTGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGaa < 2:104154/101‑1 (MQ=255) gccctttCGCGGGCGTGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGaa > 1:104154/6‑106 (MQ=255) | TGGAAGGTGGAGCTGGTGGTGGCTGACGAGCTGATTGTGCAGGTCGTCGCTGCGCTCAAGCAGAGCCACCCGTACGAGACGCCAGCCTATGAAGTCTGGCGGCTCGCCGAGTTCTAGAACGCATCACGGCTATTCGCGGGCACGCCCGCGAA‑AGGGCCCTACAATCCAGCCACAATCTCCTTGCCTGGCTCCGGCCGCTTCAACATTCCACAAGCCAGCAACCCCATCTCGAACAGCACCCACATCGGCACTGCCAGCATCGTC > NC_002947/1188900‑1189163 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |