Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 541,443 | C→T | 8.3% | P1273L (CCG→CTG) | rpoB → | DNA‑directed RNA polymerase subunit beta |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 541,443 | 0 | C | T | 8.3% | 131.1 / 3.2 | 48 | P1273L (CCG→CTG) | rpoB | DNA‑directed RNA polymerase subunit beta |
Reads supporting (aligned to +/- strand): ref base C (22/22); new base T (2/2); total (24/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01 |
TGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAA > NC_002947/541299‑541587 | tGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTgg < 2:420324/150‑1 (MQ=255) gTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGggt > 1:735598/1‑150 (MQ=255) tGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGggtg > 1:582034/1‑150 (MQ=255) aTGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGgcgc > 2:427927/1‑150 (MQ=255) gTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGGCGGTTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑Cg > 1:571964/1‑150 (MQ=255) aCCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑Cggt < 2:546693/150‑1 (MQ=255) ccGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑Cggtg < 1:36999/150‑1 (MQ=255) cGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCgg < 1:281699/150‑1 (MQ=255) gTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGgaga > 2:58101/1‑150 (MQ=255) ccTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATg > 2:405886/1‑150 (MQ=255) tgtgACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCGGCTGCTGggtaa < 1:697392/111‑3 (MQ=255) tgtgACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCGGCTGCTGggtaa > 2:697392/1‑109 (MQ=255) gTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCAAGCGAGTACGACAGGTTCCTACGCTCGGATGACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGgcgc > 1:467238/1‑150 (MQ=255) atgtacatgCTCAAGCTGAACCACTTGGTGGCAACCAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTCGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa < 2:435306/150‑1 (MQ=255) atgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGACGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa < 1:408486/150‑1 (MQ=255) tgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAg > 2:132283/1‑150 (MQ=255) cTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGgcg < 2:511753/93‑1 (MQ=255) cTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGgcg > 1:511753/1‑93 (MQ=255) cTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCa < 1:132283/150‑1 (MQ=255) aCCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCCGCCATTCCG‑GATTCGATGGGTTTGAGGTGGTCAGCGTTTCGGGCAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATaca < 2:25908/150‑1 (MQ=255) ccACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑Cgg < 2:670626/90‑1 (MQ=255) ccACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑Cgg > 1:670626/1‑90 (MQ=255) tggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAg > 1:162100/1‑150 (MQ=255) gacgacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGaa > 1:506667/1‑150 (MQ=255) gATGCACGCGCGTTCCAATGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTg > 1:389956/1‑52 (MQ=255) gATGCACGCGCGTTCCAATGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTg < 2:389956/52‑1 (MQ=255) ttCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGATATGGAAGTGTGGGCGCTGGAAGCa > 1:167108/1‑103 (MQ=255) ttCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCa < 2:167108/103‑1 (MQ=255) cACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgacga > 1:504658/1‑150 (MQ=255) cACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgacga > 1:223681/1‑150 (MQ=255) aCTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgacgac < 2:735598/150‑1 (MQ=255) ttCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGaa < 2:467238/150‑1 (MQ=255) tCCTACAGCCTGGTTACGCAGCATCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTAAAc < 2:621035/150‑1 (MQ=255) aCAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGcc > 2:517015/1‑150 (MQ=255) aGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCG‑TTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAATTCGGACGACGTGAACGGCCGTa < 1:492152/150‑1 (MQ=255) cTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTAcc > 2:113796/1‑150 (MQ=255) gTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAg < 1:574729/150‑1 (MQ=255) atctccgagcagcCGCTTGGTTGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAg < 2:232114/143‑1 (MQ=255) ccAGCGGCTGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGaa < 2:223628/117‑1 (MQ=255) ccAGCGGCTGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGaa > 1:223628/1‑117 (MQ=255) ccAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGaagat > 1:232114/1‑146 (MQ=255) ccAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTa > 2:180006/1‑150 (MQ=255) cagcagCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTAc < 2:223681/150‑1 (MQ=255) gcagcCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACCGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACaa < 1:427927/150‑1 (MQ=255) agcCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGaa > 2:227238/1‑112 (MQ=255) agcCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGaa < 1:227238/112‑1 (MQ=255) gcCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTc > 2:198580/1‑103 (MQ=255) gcCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTc < 1:198580/103‑1 (MQ=255) | TGGTGCTGTTCGATGGCCGTACCGGCAACAAGTTCGAGCGTCCTGTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTT‑CGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAA > NC_002947/541299‑541587 |
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: ‑ |