Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,170,226 | A→T | 16.6% | intergenic (‑290/‑97) | katG ← / → PP_3669 | catalase‑peroxidase/LysR family transcriptional regulator |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,170,226 | 0 | A | T | 16.6% | 55.6 / 2.7 | 24 | intergenic (‑290/‑97) | katG/PP_3669 | catalase‑peroxidase/LysR family transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base A (11/9); new base T (2/2); total (13/11) | |||||||||||
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.83e-01 |
CAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCC‑GCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGC > NC_002947/4170080‑4170368 | cAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGc < 1:408728/150‑1 (MQ=255) aGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTAACGGGCAGGCCTGTATCGGCCTCATCGCGGTCAAGcc < 1:508314/150‑1 (MQ=255) gTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCc‑g > 2:311806/1‑150 (MQ=255) cgcTGCCCTGTCGCGGAATTCGCGGGGCCATGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACC‑GCGATGAGGc < 1:571942/150‑1 (MQ=255) gcgGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCC‑GCGATGAGGCCGGCACATCCAt > 2:37070/1‑150 (MQ=255) aaTTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCC‑GCGATGAGGCCGGCACATCCATacaa < 1:311806/150‑1 (MQ=255) gcgGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCCAACGCT‑AGGCCGGCACATCCATacaaggcac > 1:308255/1‑150 (MQ=255) cACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACC‑GCGATGAg > 2:503880/1‑120 (MQ=39) cACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGGCTTGACC‑GCGATGAg < 1:503880/120‑1 (MQ=39) aCGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAACCC‑GCGCCGAGCCCGGCACATCCATACAAGGCACAAGGcccc > 1:277341/1‑150 (MQ=255) aCGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGACGACCTGGATCGGCCTCATCGCGGTCAAGCCC‑GCGATGAGGGCGGCACAGCCATCCAAGGCACAAGGcccc > 1:288940/1‑150 (MQ=255) tAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGGCGGCCTGTATCGGCCTCATCGCGGTCTAGTCC‑GCGATGAGGCCGGCACATACATACAAGGCAAAAATCTTCGCTCCATGCTTCAATTCACAGCCATCCCCtttt > 1:347536/1‑149 (MQ=255) tggctCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCC‑GCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCaa < 1:37070/150‑1 (MQ=255) ccAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCC‑GCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTgc < 1:392877/149‑1 (MQ=255) gcAGGCCTGTATCGGCCTCATCGCGGTCAAGCCC‑GCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCcgcg > 1:242044/1‑150 (MQ=255) gCCTGTATCGGCCTCATCGCGGTCAAGCCC‑GCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCtt > 1:21932/1‑80 (MQ=255) gCCTGTATCGGCCTCATCGCGGTCAAGCCC‑GCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCtt < 2:21932/80‑1 (MQ=255) tATCGGCCTCATCGCGGTCAAGCCC‑GCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGa > 1:17559/1‑150 (MQ=255) tCGGCCTCATCGCGGTCAAGCCC‑GCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAAc > 1:69862/1‑150 (MQ=255) cTCATCGCGGTCAAGCCC‑GCGATGAGGCCGGCACATCCATGCAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACctgcg < 2:288940/150‑2 (MQ=255) tCGCGGTCAAGCCC‑GCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTCCTGGTg > 2:53456/1‑150 (MQ=255) cgcgGTTAATCCC‑GCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGc < 2:308255/150‑1 (MQ=255) cgGTCAAGCCC‑GCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCtt > 1:447889/1‑81 (MQ=255) cgGTCAAGCCC‑GCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCtt < 2:447889/81‑1 (MQ=255) | CAGGTGCTGACGCGCTGCCCTGTCGCGGAATTCGCGGGGCCACGAGAAAGACTAGCCGAGCTTGGGCAGAGGTCTAATAGAGTGGCTGTTGGAGCCTGATAGGCTGGCTCTGTCAGGCGCAGGCCTGTATCGGCCTCATCGCGGTCAAGCCC‑GCGATGAGGCCGGCACATCCATACAAGGCACAAGGCCCCGCGCCATGCTTCAATACACAGCCATACCCTTTGCCCAGGACGCCGAGCCCTCATGAACCGCAACGACCTGCGCCGCGTCGACATGAACCTGCTGGTGC > NC_002947/4170080‑4170368 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |