| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 4,189,513 | +C | coding (4164/4224 nt) | rpoC → | RNA polymerase, beta prime subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,189,512 | 1 | . | C | 97.5% | 121.8 / ‑2.9 | 40 | coding (4163/4224 nt) | rpoC | RNA polymerase, beta prime subunit |
| Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (21/18); total (21/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.75e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01 | |||||||||||
ACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑A‑CGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTG > NC_000913/4189374‑4189649 | tccaaCTCCTCGGCCTGGCAGCGATCGTTATTCTGTGTCGGCTGTGCCCGGCCGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCTGCTGCACCGCAGGTGACTGCAGAAG‑A‑ccgc < 2:34783/140‑4 (MQ=255) cgcgGCCTGAAAGGGAACGTTATCGTGGGGCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCa < 2:489861/143‑1 (MQ=255) gcgGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAg > 1:610072/1‑143 (MQ=255) cgGCCTGAAAGAGAACGTTATCGTGGGTCGTCGGATCCCGGCAGGTACCGGGTAGGCGTACCACCTGGTTCGTATGCGTCGCCGTGCGGCGGGTGAAGCTCCGGCTGCACCGCGGGTGACTGCCGAAG‑ACCGCATCTGccggg > 1:404972/1‑140 (MQ=255) cagAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGCGAAGCTCCGTTTTTACCGCAGGTTACTGCAGAAG‑ACCGCATCTGCCAGCCTgg < 1:862894/141‑1 (MQ=255) tGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGc > 1:158835/1‑143 (MQ=255) tGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCGCCGGCTGCACCAGAGGTGACCGACGCAG‑ACCGCAACCGCCAGCATGGc > 2:327261/1‑143 (MQ=255) agagAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCGGaa > 1:306374/1‑143 (MQ=255) aCGTTATCGGGGGTCGTCTGGTCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAActgct < 1:526970/143‑1 (MQ=255) ttATCCGGGGTCGTCTGATGCGGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGaa < 2:748040/143‑1 (MQ=255) tCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGACCTGCTGAACGc > 1:497678/1‑143 (MQ=255) gTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCCGGTGAAGCTCCGGCGGCACCGCAGGTGACTGCAGAAG‑ACCTCATCTGCCAACCTGGCAGAAATGCTGAACGCAg > 1:58934/1‑143 (MQ=255) cgtcTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGCATCGTATGCGTCGCCGTGCTGCGCGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCCGGCAGAACTGCTGAACGCAGGTCTgg < 2:863589/143‑1 (MQ=255) ctcTGCTCCGGGCAGGGACCGGTCACGCGTACCACCAGGAGCGTATGCGTCGCCGTGCTGCGGGTGACGGCCCGGCTGCTCCGCAGGTGACTGCAGACGCC‑CGCAGCTGCCAGCCTGGCAGAACTGCTGAGCGCAGGTCTggg < 2:468216/142‑1 (MQ=255) tGATCCCGGCGGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCgg < 2:548246/143‑1 (MQ=255) agcccGGCGGGGACCGGTTACGCGGACCACGAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGCGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGtt < 1:270044/141‑1 (MQ=255) ccggggaGGTACCGGGTACGCGTACCACCAGGAGCGTATGCGTCGCTGTGGTGCGGACGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCt < 1:416376/137‑1 (MQ=255) cGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCCGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCCGa > 2:545821/1‑143 (MQ=255) cAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATaa > 2:768630/1‑143 (MQ=255) ttACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGTGTGAAGCTCCGGTTGTACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGa < 1:662498/142‑1 (MQ=255) ccGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCCACGGCTGCGCAGGCGGTGACCGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGAGGTTCTGATaaagcata > 2:789639/1‑137 (MQ=255) caGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGAATGGAAAGTGTTTAGGACGAGTAATCGATAATgcgcccat > 1:164847/1‑139 (MQ=255) caGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCGCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAAt < 1:836874/143‑1 (MQ=255) aGGATCGTAGGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATa < 2:523221/143‑1 (MQ=255) gcctccGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAAccc < 2:704609/139‑1 (MQ=255) gTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGACAACGAGTAATCGTTAATCCGCAACTAACGTAAAAAccc > 1:439365/1‑143 (MQ=255) gCCGTTCTGCGGGTCAAGCTTCGCCTGCACCCCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCt < 1:615342/143‑1 (MQ=255) ccGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTGAAAGCCCGCtt > 1:433527/1‑143 (MQ=255) gtgcGGGTTAAGCTCCGGCTGCACAGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTcggcg < 1:301669/142‑1 (MQ=255) cGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTCAAAACCCGCTTCGGCgggg > 2:483346/1‑142 (MQ=255) cGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGt < 1:888890/143‑1 (MQ=255) aagcctcGGCTGCATCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGCGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGttttttt < 1:239325/137‑1 (MQ=255) gCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTAt > 2:521512/1‑143 (MQ=255) ccGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATggg > 1:91664/1‑143 (MQ=255) ccGGCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGCTTTTTTATggg > 1:389035/1‑143 (MQ=255) ggCTGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGCTCCTGATAACGAGTACTCGTTAATCCGCAAATAACGTAAACACCCGCTTCGGACGGTTTTTTTATgggac > 2:680132/1‑141 (MQ=255) tGCACCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGACTTTTTATGGGGGGAg > 2:272403/1‑143 (MQ=255) aCCGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAACAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGttt > 2:267239/1‑143 (MQ=255) ccGCAGGTGACTGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTa > 1:525419/1‑143 (MQ=255) cGCAGGTGACGGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATTACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAg < 1:789639/143‑1 (MQ=255) tGCAGAAG‑ACCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCa < 1:768630/143‑1 (MQ=255) gaag‑acCGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGCGCATTTg > 1:455170/1‑143 (MQ=255) | ACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑A‑CGCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTG > NC_000913/4189374‑4189649 |
| 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: ‑ |