Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 2,719,426 | T→C | 38.5% | intergenic (‑321/+2) | kgtP ← / ← rrfG | alpha‑ketoglutarate transporter/5S ribosomal RNA of rrnG operon |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 2,719,426 | 0 | T | C | 38.5% | 5.9 / 15.9 | 24 | intergenic (‑321/+2) | kgtP/rrfG | alpha‑ketoglutarate transporter/5S ribosomal RNA of rrnG operon |
Reads supporting (aligned to +/- strand): ref base T (9/4); new base C (4/5); total (13/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.84e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01 |
CGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGACCGC > CP009273/2719281‑2719568 | cGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCTTTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGAt < 1:153012/149‑1 (MQ=255) aGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGAGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGGTCCCTa > 1:203332/1‑149 (MQ=255) gtaaGCCCTATCCGTTGGTGAATTTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTTAGGATGGCCTTTCGCATAATTGGATGCCTGGCAGTTCCCTACTCTCg < 1:80336/147‑1 (MQ=25) tCCGTTGCTGAATCTGGCAAATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGACAAATGCCTGTATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGgaga > 1:241490/1‑149 (MQ=25) tGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGGGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCca < 2:21605/149‑1 (MQ=255) cATATGTGGGAGTATAGGTCGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCTATCCGTCAGGATGGCCTTTCGCATAATTCGATGTCTGGTAGTTACCTACTCTCGCATGGGGAGACCCCACACTACCATCgg < 1:242333/149‑1 (MQ=25) gtgGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATCATTTGATGCCTGGCAGTTCCCTACTCTCGCCTGGGGAGCCCCCACACTACCCTCggacgta > 2:202706/1‑144 (MQ=21) tgGGAGTATAAGACGCGCAGCGTCGCATCCGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAc > 1:199072/1‑149 (MQ=255) tgGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTTTTTTTTTTTGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAc > 2:50674/1‑149 (MQ=25) gattctcacaCGCGCAGCGTCGCATCAGGCCTTTTTTTCTCCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGc < 2:254169/140‑1 (MQ=25) aGTATAAGACGCGCAGCGTCGCATCAGGCTTGTCTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTTGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCg < 1:128485/149‑1 (MQ=37) aCGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACtt < 2:138600/149‑1 (MQ=255) gTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAAGATGGCCTTTAGCATAATACGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTca > 1:72507/1‑148 (MQ=255) tcgggCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGGCATGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg < 2:91173/146‑1 (MQ=255) tttttttCTGCGCCAATGCAACAAGGCCATCCGTCAGGATGACATTTAGCATAATTTGATGCCTGACAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTAAgg > 1:171735/1‑149 (MQ=18) tttttttCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCGGAGTTCGGCATGGGGTCAgg > 2:220144/1‑149 (MQ=255) cAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcgcg < 1:262603/149‑1 (MQ=255) gCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGCCAGGTGGGACCACCGCGCTAc > 1:268809/1‑149 (MQ=255) ggCCATCCGTCAGGATGGCCTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATAGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGgccgcc > 1:112399/1‑149 (MQ=255) ggtaTGTCCTTTCGCATAAGTTGATGCCTGGCAGTTCCATTTTCGAGCATGGGGATACCCCACACTACCACCGGCGCTACGGCGTTTCTCTTTTGAGTTCGGCATGCGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCt < 2:90671/146‑1 (MQ=11) agtggggCCTTTCGCATAAGTGGATGCCGGGCAGTTCGCTGCACGCGCATGGGGAGACCCCACACTACTAGCGGGGCTACGGCGTGTCACTTCTGAGTGCGGCATGGGGTCAGGTGGGACCCCCGCGCTACGGCCGGCAGGCTAATTCt < 2:90661/144‑1 (MQ=11) cTTTCGCATAATTCGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGGCAGGCAAATTCTGTTTTATc > 1:316256/1‑149 (MQ=255) ttCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAg > 1:35765/1‑149 (MQ=34) aTAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGACCGc > 1:11817/1‑149 (MQ=11) | CGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCAGGATGGCCTTTCGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAATTCTGTTTTATCAGACCGC > CP009273/2719281‑2719568 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |