Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,866,753 | A→C | 9.9% | intergenic (+20/+37) | PP_5140 → / ← thyA | MerR family transcriptional regulator/thymidylate synthase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,866,753 | 0 | A | C | 9.9% | 131.0 / 6.4 | 50 | intergenic (+20/+37) | PP_5140/thyA | MerR family transcriptional regulator/thymidylate synthase |
Reads supporting (aligned to +/- strand): ref base A (22/23); new base C (2/3); total (24/26) | |||||||||||
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.82e-01 |
GTGGAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCGG > NC_002947/5866605‑5866899 | gtggAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGaa > 1:1109194/1‑150 (MQ=255) tGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGa < 2:637067/150‑1 (MQ=255) gcggcggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCCGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCCGCGaa < 1:872197/150‑1 (MQ=255) cggcggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAag < 2:1109194/150‑1 (MQ=255) ggcggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAaga < 1:1075002/150‑1 (MQ=255) ggcggcgCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAaga < 1:52545/150‑1 (MQ=255) aCGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTaa < 1:679899/150‑1 (MQ=255) gaGGGGAGTGAATGCGCGATCTTGCAGCCACTGGAGACAAACGGGCCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATa < 2:1102563/149‑1 (MQ=255) aGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACTGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCTGTTATAACAATat > 2:404202/1‑149 (MQ=255) ggTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAgagagc < 2:1029658/149‑3 (MQ=255) gAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCTGCACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCACAAGCTGCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTAt < 1:633859/150‑1 (MQ=255) aaaTTCGCGATCTCGCAGCCACTGGAGACCAACGGGGCGGTATCGTTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTGTAACAATAGAGATCTAt < 2:615342/149‑1 (MQ=255) aTGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCa < 2:57995/150‑1 (MQ=255) gATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGc > 1:429175/1‑150 (MQ=255) tGGAGACAAACGGAGCGGTATCGGTACCGGATACCGAACATTCCCATGCGGGGCGGAGTAACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCAGCGACGAGGTCAGGTCTAACAATAGAGATCTATCAGACCGGCAT‑CGTCGCTGTCa > 2:99443/1‑150 (MQ=255) aaaCGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTATAACAATAGAGAGCTATCAGACCGCCAT‑CGGCGCTGTCATCGgcgc > 2:70434/1‑150 (MQ=255) aaaCGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTAAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAg < 1:1151206/124‑1 (MQ=255) aaaCGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTAAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAg > 2:1151206/1‑124 (MQ=255) aCGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCAt > 2:875071/1‑150 (MQ=255) aCGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCAt < 1:474302/150‑1 (MQ=255) ggATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTccag > 2:111690/1‑150 (MQ=255) gATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCCGCGaa < 2:475146/71‑1 (MQ=39) gATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTTCACCCGCGaa > 1:475146/1‑71 (MQ=255) gATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTccagc < 2:509790/150‑1 (MQ=255) ccGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGa > 1:583968/1‑150 (MQ=255) cGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCt < 1:951547/112‑1 (MQ=255) cGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCt > 2:951547/1‑112 (MQ=255) gttCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGGCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCAc > 2:826057/2‑150 (MQ=255) tCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATaga < 2:3121/80‑1 (MQ=255) tCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATaga > 1:3121/1‑80 (MQ=255) ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATa > 1:451596/1‑76 (MQ=255) ccATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATa < 2:451596/76‑1 (MQ=255) ggggCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAAt > 2:109678/1‑150 (MQ=255) gCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTg > 2:726918/1‑111 (MQ=255) gCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTg < 1:726918/111‑1 (MQ=255) gAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTtctc < 2:429175/150‑1 (MQ=255) gTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCa < 1:111690/150‑1 (MQ=255) cACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGc < 1:875071/150‑1 (MQ=255) aCGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCTGACGGCCATACTGCGC‑CAGATGGGCGCAGGGTTCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATTTTCTCCAGcc > 1:390648/1‑150 (MQ=255) aTTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTATCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCg > 2:20912/1‑150 (MQ=255) ccAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCAt < 1:109678/150‑1 (MQ=255) aaaaGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTAt < 1:1135307/54‑1 (MQ=255) aaaaGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTAt > 2:1135307/1‑54 (MQ=255) aaaaGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATAc > 1:343239/1‑150 (MQ=255) aaaaGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATAc > 1:596717/1‑150 (MQ=255) aaaGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACt > 1:535340/1‑150 (MQ=255) gCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGc < 2:343239/150‑1 (MQ=255) cTTCGCGGGTTCACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTccagcca < 2:837136/102‑1 (MQ=255) cTTCGCGGGTTCACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTccagcca > 1:837136/1‑102 (MQ=255) ttCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCa > 1:848514/1‑150 (MQ=255) gTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCgg > 2:142305/1‑150 (MQ=255) | GTGGAGCTGCGGCGGCGCTGCAACGCGCAGGGGAGTGAATGCGCGATCTTGCAGCAACTGGAGACAAACGGGGCGGTATCGGTACCGGATACCGAACATTCCCATGTGGGGCGGAGTCACGGGCATTGAGCCAAAAGCTTCGCGGGTGAACCCGCGAAGAGGTCAGTTCTAACAATAGAGATCTATCAGACCGCCAT‑CGGCGCTGTCATCGGCGCATGGTGCTCATAGCCTTCCAGCCAGAAATCACTCGGCTCAATCTTCTCCAGCCACTCCGGCTCATACTTGCCAGTTTCGG > NC_002947/5866605‑5866899 |
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: ‑ |