Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,731,312 | G→A | 100% | intergenic (‑155/+288) | rrsG ← / ← clpB | 16S ribosomal RNA of rrnG operon/protein disaggregation chaperone |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,731,312 | 0 | G | A | 92.1% | 88.6 / 1.1 | 38 | intergenic (‑155/+288) | rrsG/clpB | 16S ribosomal RNA of rrnG operon/protein disaggregation chaperone |
Reads supporting (aligned to +/- strand): ref base G (2/1); new base A (11/24); total (13/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.65e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01 |
ACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAC > NC_000913/2731241‑2731406 | aCGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTTGCTTTTTTTTTTTGTTGCGGGGTGTGCa > 1:1288144/1‑99 (MQ=11) aacttagcgTGACCACACAGATAGTCTGATAAATTGTTAAAGCGCCGTGCCGCTTCACTTTTTATAAGCGGCGCGGGGTGCGCATACTACGACTTCCCg > 1:1731988/9‑99 (MQ=0) aGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTc > 1:942919/1‑100 (MQ=11) aGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTc < 1:1437894/100‑1 (MQ=11) gATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCt < 1:1878959/100‑1 (MQ=14) ttGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCttt < 1:1412165/100‑1 (MQ=17) tGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg < 2:1979054/100‑1 (MQ=18) tGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg > 2:1464215/1‑100 (MQ=25) tGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg < 2:2904552/100‑1 (MQ=18) tGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg < 1:128391/100‑1 (MQ=18) tGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg < 1:902308/100‑1 (MQ=18) tGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCACCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg < 1:1348108/100‑1 (MQ=16) gTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGc < 2:567748/100‑1 (MQ=25) gTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGc < 2:324447/100‑1 (MQ=25) cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGAGGGGTGTGCATAATACGCCTTCTCGCTACAGAGGCAAGCATTTCTTTGCtt > 1:1699136/1‑100 (MQ=11) tGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCttt > 1:2361937/1‑100 (MQ=255) tGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTGTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCttt < 2:715033/100‑1 (MQ=255) gATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtttt > 1:500315/1‑100 (MQ=255) aTAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtc < 2:1208985/100‑1 (MQ=255) aTAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtc < 2:975602/100‑1 (MQ=255) aaaTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtctc < 1:622026/100‑1 (MQ=255) aaaTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtcac < 2:1973647/100‑3 (MQ=255) aTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTg < 2:125045/100‑1 (MQ=255) ttGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGACGCGGGGTGGGCATAATACGTCTTCCCGCTAGAGAGGCAAGCATTTATTTGCTTTTCTCGGt > 1:2682459/1‑100 (MQ=11) tgtgAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGt < 1:2661696/99‑1 (MQ=255) tGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGtt < 1:103182/100‑1 (MQ=255) gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg < 2:2560870/100‑1 (MQ=255) ttAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTga < 2:669936/100‑1 (MQ=255) tAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgag > 2:485552/1‑100 (MQ=255) tAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgag > 2:2123618/1‑100 (MQ=255) gagCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtc > 1:769474/1‑100 (MQ=255) agCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtct < 2:2877533/100‑1 (MQ=255) agCAGTGCCGCTTCACTTTTTCTCAGCGACGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtct > 2:1668536/1‑100 (MQ=255) cAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCa < 2:445497/100‑1 (MQ=255) cAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCa < 2:1601112/100‑1 (MQ=255) aGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg > 1:603295/1‑100 (MQ=255) aGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg < 1:1149133/100‑1 (MQ=255) gTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg < 1:2876971/100‑1 (MQ=255) cgcttcacttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGCGaa > 2:1930783/1‑100 (MQ=255) cgcttcacttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGaa < 2:870531/100‑1 (MQ=255) cgcttcacttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGaa > 2:915510/1‑100 (MQ=255) gcttcacttTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAc < 1:2084614/100‑1 (MQ=255) | ACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAC > NC_000913/2731241‑2731406 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |