Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,915,892 | T→C | 10.5% | Y1623Y (TAT→TAC) | PP_2561 → | hemolysin‑type calcium‑binding bacteriocin |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,915,892 | 0 | T | C | 10.5% | 99.7 / 3.0 | 38 | Y1623Y (TAT→TAC) | PP_2561 | hemolysin‑type calcium‑binding bacteriocin |
Reads supporting (aligned to +/- strand): ref base T (16/18); new base C (2/2); total (18/20) | |||||||||||
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.00e-01 |
GCCAGCCGGTTTCTTCGGCCCAGGTTCCCCCGCGATCCCGGGCTCGTCGTACCAGCAGACCAGCGGCCCGGTGTTCGATTCTCAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCATCGTCGCCAGCCCGGGCATGGACGGCCAGTTCGGCACCACCGACGACCATGACG > NC_002947/2915751‑2916031 | gcCAGCCGGTTTCTTCGGCCCAGGTTCCCCCGCGATCCCGGGCTCGTCGTACCAGCAGACCAGCGGCCCGGTGTTCGATTCTCAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACcgc < 2:230247/150‑1 (MQ=255) gcCAGCCGGTTTCTTCGGCCCAGGTTCCCCCGCGATCCCGGGCTCGTCGTACCAGCAGACCAGCGGCCCGGTGTTCGATTCTCAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTACGCCACcgc > 2:76676/1‑150 (MQ=255) gcCAGCCGGTTTCTTCGGCCCAGGTTCCCCCGCGATCCCGGGCTCGTCGTACCAGCAGACCAGCGGCCCGGTGTTCGATTCTCAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTACGCCACcgc > 2:77177/1‑150 (MQ=255) gCCGGTTTCTTCGGCGCTGGTTCCCGCGCGATCCCGTGCTCGTCTTACCAGCAGACCAGCGGCCCGGTGTTCGATTCTCAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTAt < 2:560817/150‑1 (MQ=255) cGGTTTACTCGGCACAGGTACCCCCGCGATCCCTGTCTTGTCGTACCAGCAGACCAGCGGCCCGGTGTTCGATTCTCAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGa < 2:561011/150‑1 (MQ=255) ttcttcGGCCCAGGTTCCCCCGCGATCCCGGGCTCGTCGTACCAGCAGACCAGCGGCCCGGTGTTCGATTCTCAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCg < 1:1070359/150‑1 (MQ=255) gaaGGTT‑CCCCGGGATCCCGGGCTCGTCGTACCAGCAGACCAGCGGCCCGGTGTTCGATTCTCAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCTGGGTGCCGATg > 2:539978/3‑150 (MQ=255) ccAGGTTCCCCCGCGATCCCGGGCTCGTCGCACCAGCAGACCAGCGGGCCGGTGTTCAATTCTCAGCCACGTACCATCAACAACGTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGATACGGGTACCGAt > 2:480273/1‑150 (MQ=255) ccccGCGATCCCGGGCTCGTCGTACCAGCAGACCAGCGGCCCGGTGTTCGATTCTCAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTACGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCt < 1:76676/150‑1 (MQ=255) ccccGCGATCCCGGGCTCGTCGTACCAGCAGACCAGCGGCCCGGTGTTCGATTCTCAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTACGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCt < 1:77177/150‑1 (MQ=255) tCCCGGGCTCGTCGTACCAGCAGACCAGCGGCCCGGTGTTCGATTCTCAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCg < 1:1139278/150‑1 (MQ=255) agcagACCAGCGGCCCGGTATTCGATTCTCAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAacgacg > 2:889650/1‑150 (MQ=255) cagACCAGCGGCCCGGTGTTCGATTCTCAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAacgacgac < 1:413021/150‑1 (MQ=255) gCCCGGTGTTCGATTCTCAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAgg > 1:932717/1‑150 (MQ=255) tCGATTCTCAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTgc < 1:904728/150‑1 (MQ=255) tctcAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCAtc < 1:277058/150‑1 (MQ=255) tcAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCAtcgt > 2:1110695/1‑150 (MQ=255) cAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCAtcgtc < 1:1162534/150‑1 (MQ=255) aGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCAtcgtcg > 2:664413/1‑150 (MQ=255) cACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCATCGTCGCCa < 2:932717/150‑1 (MQ=255) cACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCATCGTCGCCa < 1:1110695/150‑1 (MQ=255) cGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCATCGTCGCCAGc > 1:897531/1‑150 (MQ=255) gTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGAt > 2:15774/1‑128 (MQ=255) gTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGAt < 1:15774/128‑1 (MQ=255) cATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCAtcgtc < 2:167344/139‑1 (MQ=255) cATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCAtcgtc > 1:167344/1‑139 (MQ=255) gCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTc > 2:419873/1‑86 (MQ=255) gCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTc < 1:419873/86‑1 (MQ=255) ccTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCATCGTCGCCAGCCCGGGCATGGACg > 1:726393/1‑149 (MQ=255) ccTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCATCGTCGCCAGCCCGGGCATGGACg < 2:726393/149‑1 (MQ=255) aCCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCATCGTCGCCAGCCCGGGCATGGACGGCCAGTTCGGca > 1:830084/1‑150 (MQ=255) aCCTCCAACAACCCGGAAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCATCGTCGCCAGCCCGGGCATGGACGGCCAGTTCGGCaccacc > 2:357915/1‑150 (MQ=255) ccTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCATCGTCGCCAGCCCGGGCATGGACGGCCAGTTCGGCACCACcg < 1:889650/150‑1 (MQ=255) cTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCATCGTCGCCAGCCCGGGCATGGACGGc > 2:490593/1‑133 (MQ=255) cTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCATCGTCGCCAGCCCGGGCATGGACGGc < 1:490593/133‑1 (MQ=255) aacaacCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCaa > 1:25872/1‑89 (MQ=255) aacaacCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCaa < 2:25872/89‑1 (MQ=255) cgtcAGCCTATGCCACCGCGTATGAGCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGTCGTGCTCAAGGATGGCGTGCGCATCGGCGCCAGCCCGGGCCTGGGCGGCCGGTTCGAGACCACCGACGACCATGACg > 2:21528/4‑150 (MQ=255) | GCCAGCCGGTTTCTTCGGCCCAGGTTCCCCCGCGATCCCGGGCTCGTCGTACCAGCAGACCAGCGGCCCGGTGTTCGATTCTCAGCCACGTACCATCAGCAACCTGATTGTCGACCAGACCTCCAACAACCCGGCAGCCTATGCCACCGCGTATGACCCGGGTGCCGATGGCGTGCTCAACTTCGGCGCCGCAGGCAACGACGACGTGCTCAAGGATGGCGTGCGCATCGTCGCCAGCCCGGGCATGGACGGCCAGTTCGGCACCACCGACGACCATGACG > NC_002947/2915751‑2916031 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |