Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,731,312 | G→A | 64.0% | 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 | 64.0% | 51.0 / 59.5 | 64 | intergenic (‑155/+288) | rrsG/clpB | 16S ribosomal RNA of rrnG operon/protein disaggregation chaperone |
Reads supporting (aligned to +/- strand): ref base G (21/2); new base A (18/23); total (39/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.62e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01 |
TTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAC > NC_000913/2731229‑2731406 | cctgagatctccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAgcggcg > 1:986902‑M1/76‑100 (MQ=255) cctgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAgcggcg > 2:1322126‑M1/76‑100 (MQ=255) cctgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAgcggcg > 1:1774628‑M1/76‑100 (MQ=255) cctgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAgcggcg > 1:2382563‑M1/76‑100 (MQ=255) cctgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAgcggcg > 1:2403811‑M1/76‑100 (MQ=255) tgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGcgcg > 2:2452401‑M1/74‑100 (MQ=255) tgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGcgcg > 2:110763‑M1/74‑100 (MQ=255) tgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGcgcg > 1:2143905‑M1/74‑100 (MQ=255) tgagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAGCGGcgcg > 2:2020790‑M1/74‑100 (MQ=255) agagacctccggcagaaaaagcaaaaataaatgcttgactctgtagcgagaaggcgtattatgcacaccccgcGCCACTTCGCTTTTTCTCAGCGGcgcg > 1:1789093‑M1/74‑100 (MQ=255) gagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCgg > 2:673386‑M1/73‑100 (MQ=255) gagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCgg > 2:730851‑M1/73‑100 (MQ=255) gagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtatgatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCgg > 1:1736646‑M1/73‑100 (MQ=255) gagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAGCGGCGCgg > 1:399210‑M1/73‑100 (MQ=255) gagaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAGCGGCGCgg > 2:22915‑M1/73‑100 (MQ=255) agaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCgg > 1:1468472‑M1/72‑99 (MQ=255) agaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAGCGGCGCggg > 1:50697‑M1/72‑100 (MQ=255) gaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaggcgtattatgcacaccccgcGCCGCTTCGCTTTTTCTCAGCGGCGCgggg > 1:954581‑M1/71‑100 (MQ=255) gaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAGCGGCGCgggg > 2:309696‑M1/71‑100 (MQ=255) gaactccggcagagaaagcaaaaataaatgcttgactctgtagcgggaaagcgtattatgcacaccccgcGCCGCTTTGCTTTTTCTCAGCGGCGCgggg > 1:941735‑M1/71‑100 (MQ=255) ttGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTTTTCGGCTCGGCgtg > 1:449530/1‑100 (MQ=1) cgcgTGCCCACACCGCTTGTCTGATCAACTGCTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTAc < 2:483026/95‑1 (MQ=1) aTTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTcat < 2:1759290/100‑3 (MQ=2) ttGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCttt < 2:1373540/100‑1 (MQ=17) tGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg < 1:1012704/100‑1 (MQ=18) ggtcatgcgtattatgcacaccccgtGCCGCTTCACGGTTTCTCAGCGGCGTGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTg > 1:2388905‑M1/27‑100 (MQ=255) gTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGc < 1:384379/100‑1 (MQ=25) gTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGc > 2:2133235/1‑100 (MQ=25) tCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCt < 2:2039152/100‑1 (MQ=255) tCTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGATACAGAGTCAAGCATTTCTTTGCt < 1:385642/100‑1 (MQ=255) cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt > 1:2545089/1‑100 (MQ=255) cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt > 1:708256/1‑100 (MQ=255) cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt > 1:2099100/1‑100 (MQ=255) cTGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCtt > 2:1350484/1‑100 (MQ=255) tGATAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGAGGCGAGGGGTGTGCATAATACGCCTTCTCGCTACAGAGTCAAGCCTTTCTTTGCttt > 1:1630198/1‑100 (MQ=14) tAAATTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTtct > 2:2519061/1‑100 (MQ=255) aaTTGTTAAAGAGCAGTGCCACTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTctct < 2:680372/100‑1 (MQ=255) aTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTg > 1:637812/1‑100 (MQ=255) aTTGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTg > 2:2375262/1‑100 (MQ=255) ttGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGt < 2:756721/100‑1 (MQ=255) tGTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGtt < 2:715012/100‑1 (MQ=255) gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATTCGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg > 1:564117/1‑100 (MQ=255) gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg < 2:270618/100‑1 (MQ=255) gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg < 1:2021714/100‑1 (MQ=255) gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg > 2:2225052/1‑100 (MQ=255) gTTAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg < 2:637815/100‑1 (MQ=255) atgcacaccccgcGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg > 2:740368‑M1/14‑100 (MQ=255) ttAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCACGGGGTGTGCATAATACGCTTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTg > 1:101309/1‑99 (MQ=255) tAAAGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTgag > 2:495745/1‑100 (MQ=255) aaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAt < 2:225987/100‑1 (MQ=255) aaGAGCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGAt < 1:1434720/100‑1 (MQ=255) gagCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtc > 2:1339572/1‑100 (MQ=255) agCAGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATtct < 1:1194935/100‑1 (MQ=255) aGTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAg < 1:89197/100‑1 (MQ=255) gTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg < 2:2145913/100‑1 (MQ=255) gTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg > 2:2112032/1‑100 (MQ=255) gTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg < 1:1291183/100‑1 (MQ=255) gTGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAgg > 1:1056983/1‑100 (MQ=255) tGCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGga < 2:574366/100‑1 (MQ=255) GCCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGgag < 1:2188894‑M1/100‑1 (MQ=255) GCCGCTTCACTTTTTCTCAGCGGCGCGGGGAGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGgag < 1:393772‑M1/100‑1 (MQ=255) CCGCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGgaga < 1:2498747‑M1/100‑1 (MQ=255) GCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAc < 2:923658‑M1/100‑1 (MQ=255) GCTTCACTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAc < 2:1079234‑M1/100‑1 (MQ=255) | TTTCGTCTTGCGACGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTGCCGCTTCGCTTTTTCTCAGCGGCGCGGGGTGTGCATAATACGCCTTCCCGCTACAGAGTCAAGCATTTCTTTGCTTTTCTCTGTTGAGATTCTCAGGAGAAC > NC_000913/2731229‑2731406 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |