Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,313,871 | T→G | 17.8% | intergenic (+33/‑126) | hemB → / → proP | 5‑aminolevulinate dehydratase/osmosensory proline/betaine/H+ permease |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,313,871 | 0 | T | G | 17.8% | 83.3 / 19.0 | 45 | intergenic (+33/‑126) | hemB/proP | 5‑aminolevulinate dehydratase/osmosensory proline/betaine/H+ permease |
Reads supporting (aligned to +/- strand): ref base T (14/23); new base G (4/4); total (18/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.94e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01 |
TGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAAG > NC_002947/3313727‑3314016 | tGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGc < 1:550874/150‑1 (MQ=255) aTCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCGGCt > 2:122007/1‑150 (MQ=255) cgacgaAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCttt < 1:252193/150‑1 (MQ=255) gacgaAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCtttt < 1:1260800/150‑1 (MQ=255) cgaAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCGGCtttcat < 2:87310/150‑4 (MQ=255) cgaAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCAATGGCAGCGGCtttcat < 1:46844/150‑4 (MQ=255) ggCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTtgt < 2:878666/150‑1 (MQ=255) cGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTgcg > 2:198886/1‑150 (MQ=255) tggtgCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTACGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGTGCCGt < 2:23165/150‑1 (MQ=255) aaaGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGgtggtg > 2:680207/1‑150 (MQ=255) aaGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGc < 2:540123/150‑1 (MQ=255) aTTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCa > 1:491542/1‑150 (MQ=255) ggCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAAc < 1:2444/150‑1 (MQ=255) gcgcAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACg < 2:239798/150‑1 (MQ=255) cgcAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGc > 1:1251225/1‑150 (MQ=255) gCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTCCGCCGTTGGGTGGTGCCGCAAACGCACaaaaaaa > 1:514068/1‑150 (MQ=255) gTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACa < 2:514068/150‑1 (MQ=255) tGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAg < 2:568653/150‑1 (MQ=255) gCTGGTGCAGACCTGATCCTGACCTACTTCGCCATAGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGt < 1:680207/150‑1 (MQ=255) cTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTa > 2:533387/1‑150 (MQ=255) aCCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAgtgt > 2:991775/1‑150 (MQ=255) ttCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCtgtt < 1:484103/150‑1 (MQ=255) ttCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCtgtt < 2:491542/150‑1 (MQ=255) ttCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTAATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCtgtt < 2:1115661/150‑1 (MQ=255) ttgACGTGGCGCGCGAAGGTATCTGCAGCCAGGTGTTGGTCACGAAAGCCGCTGCCATTGCAGCGGCGGTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTg < 1:536778/148‑1 (MQ=255) aCCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATc > 2:508718/1‑150 (MQ=255) gcgAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGc > 1:952921/1‑150 (MQ=255) aaGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCAtt > 2:50705/1‑150 (MQ=255) aGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTc < 1:448079/150‑1 (MQ=255) gTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCAt < 1:396565/150‑1 (MQ=255) tATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATc < 1:50705/150‑1 (MQ=255) tGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACa < 1:991775/150‑1 (MQ=255) cAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATa > 2:927481/1‑150 (MQ=255) aGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATAc < 2:952921/150‑1 (MQ=255) cTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACggat < 1:533387/150‑3 (MQ=255) gTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGa > 2:1196035/1‑150 (MQ=255) aaagaaaaGCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTa > 1:934079/6‑119 (MQ=39) aaagaaaaGCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTa < 2:934079/114‑1 (MQ=39) gCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAAc > 2:519559/1‑66 (MQ=255) gCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAAc < 1:519559/66‑1 (MQ=255) gCCGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCa > 1:1088519/1‑150 (MQ=255) ccGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAggg < 2:11010/128‑1 (MQ=255) ccGCTGCAATGGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAggg > 1:11010/1‑128 (MQ=255) ccATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAAg < 1:508718/150‑1 (MQ=255) ccATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAAg > 1:26079/1‑150 (MQ=255) | TGCCATCGACGAAGGCCGGGTGGTGCGCGAAAGCATTGGCGCAATCAAGCGTGCTGGTGCAGACCTGATCCTGACCTACTTCGCCATGGACCTGGCGCGCGAAGGTATCTGAAGCCAGGCCTTGGTCATGAAAGCCGCTGCCATTGCAGCGGCTTTTCTTTGTGCGCCGTTGGGTGGTGCCGCAAACGCACAAAAAAACAGTAAGCGCTGAATTATTTAGTGTCCTGTTGTATCTGATCCCTTAGAGCATTCATCGACAGGGATACGGGTATGAAACCTCGCAAGAAAAG > NC_002947/3313727‑3314016 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |