Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,812,229 | A→T | 12.5% | intergenic (+104/‑289) | PP_4238 → / → PP_5665 | hypothetical protein/McbA family microcin B17‑like processing protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,812,229 | 0 | A | T | 12.5% | 89.2 / 4.2 | 32 | intergenic (+104/‑289) | PP_4238/PP_5665 | hypothetical protein/McbA family microcin B17‑like processing protein |
Reads supporting (aligned to +/- strand): ref base A (13/15); new base T (2/2); total (15/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01 |
CACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAAC > NC_002947/4812083‑4812368 | cacaTTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAAc < 2:484442/150‑1 (MQ=255) ggTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCCCCTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGc < 2:104949/150‑1 (MQ=255) aGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGcccc < 2:156809/150‑1 (MQ=255) gCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACTCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCtt > 1:601237/1‑150 (MQ=255) ccGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCtta > 1:169433/1‑150 (MQ=255) ttCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCttatta > 1:245926/1‑150 (MQ=255) ccGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaa > 2:699427/1‑150 (MQ=255) cGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACaaa < 2:169433/150‑1 (MQ=255) cGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAAGTCTTCGGAGGATGCTTATATCTAACAACAGTATGTTCGGGCCCTTAATTAAAATTTCGCTCCGCGCATAACaaa > 2:39527/1‑150 (MQ=255) atcTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTc < 2:599160/150‑1 (MQ=255) tcACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTGCACACCGCGCATAACAAAGAACACTcaca < 1:318130/150‑1 (MQ=255) aCCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTc < 1:495129/150‑1 (MQ=255) cAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACcaca < 2:245926/150‑1 (MQ=255) aTCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGCAt > 2:443535/1‑80 (MQ=38) aTCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTTTGACATCAGCAt < 1:443535/80‑1 (MQ=38) tCCCCACGCAGGTCTAACCCTGCAACTNAAAATATCTAGANACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTc > 1:750278/1‑150 (MQ=255) ccccACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTcacca < 1:639622/139‑1 (MQ=255) ccccACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTcacca > 2:639622/1‑139 (MQ=255) ccACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCa > 1:703334/1‑150 (MQ=255) cACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCAc < 2:406290/150‑1 (MQ=255) tAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCa < 1:463326/150‑1 (MQ=255) aaCCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAg > 2:54152/1‑150 (MQ=255) aCCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGa < 1:54152/150‑1 (MQ=255) ccTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAg > 2:311666/1‑150 (MQ=255) cTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACAGCGCATGCGGTCTGGCACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGc > 2:209341/1‑150 (MQ=255) aCTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGcttctt < 2:422709/150‑1 (MQ=255) tAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCAAGTGGCGGTTAGAAGCTTCTTACATAAATATaa < 1:209341/150‑1 (MQ=255) aGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCa < 1:74655/76‑1 (MQ=255) aGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCa > 2:74655/1‑76 (MQ=255) aaCAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATaagaag > 1:546499/1‑150 (MQ=255) aGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTaa < 2:601237/150‑1 (MQ=255) tGCTGATGTCAAACATCAGCATGCTCGCGCCCCGTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAAc > 1:696978/1‑150 (MQ=255) aaaCATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGcttc > 2:269311/1‑105 (MQ=255) aaaCATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGcttc < 1:269311/105‑1 (MQ=255) | CACATTGGTGATTGGACGCAGGCCGTTCTTTCGCAAGCAGTAAATCTCCCGAAATCAATCTCACGACCCGGCAAATCTATCCCCACGCAGGTCTAACCCTGCAACTCAAAATATCTAGAAACAATTCTTCGGAGCATGCTGATGTCAAACATCAGCATGCTCGCGCCCCTTATTAAAACTTCACACCGCGCATAACAAAGAACACTCACACGCTCACCACAGAGCTTTCGCACGCCAGAATCAGAAGCTTCTTACATAAATATAAGAAGCACTGAAAACTAAAAAC > NC_002947/4812083‑4812368 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |