| 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 | 96.0% | 146.3 / ‑0.7 | 50 | coding (4163/4224 nt) | rpoC | RNA polymerase, beta prime subunit |
| Reads supporting (aligned to +/- strand): ref base . (0/2); new base C (28/20); total (28/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.89e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑A‑C‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATT > NC_000913/4189371‑4189647 | gcgACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑A‑c > 1:934731/1‑143 (MQ=255) aCTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCAGTTACGCGTACCACCAGGATCGCATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCt > 1:930417/1‑143 (MQ=255) gcggcgCGAAAGGGCGCGTTAGCGTGGGTCGTCTGGTCCCGGCAGGTCCCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAg < 2:877950/143‑1 (MQ=255) ccTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCCGAAG‑ACC‑GCCTCTGCCAGccgg > 1:408492/1‑141 (MQ=255) ccTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTg > 1:754637/1‑143 (MQ=255) acgAAGGAGTACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGGACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTgg < 1:717552/141‑1 (MQ=255) tGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTgcc > 1:523407/1‑141 (MQ=255) gaACGTGACCGTGGGTCGGCTGAGCACGGCAGGTACCGGTTGGGCGGACCAGCAGGATCGTATGCGTCGCCGTGCGGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAActg < 2:676254/143‑1 (MQ=255) accGTTCTCGTGGGTCGTCTGCTCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGCAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAActgc < 2:844996/141‑1 (MQ=255) ttgtctggggTCGTCCGTCCCCGGCGGGTACCGGTTGCGCGTGCCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGaa < 2:742852/136‑1 (MQ=255) tCGTGGGTCGTCTGAGCCCGGCAGGTACCGGTGACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGc < 1:317313/143‑1 (MQ=255) gTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAg < 2:592580/143‑1 (MQ=255) gTGGGGCGTCTGAGCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGGGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAg < 1:696026/143‑1 (MQ=255) tGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGGCGTGCTGCGGGTGAAGCTGCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGACCTGCTGCAAGCAgg > 1:867546/1‑143 (MQ=255) cgtgtcaTCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTgg < 1:434674/137‑1 (MQ=255) tcTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGc > 1:266537/1‑143 (MQ=255) tcTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCACAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGc > 1:474696/1‑143 (MQ=255) tgggtCCCGGCGGGTAGGGGTTAGGCGTCCCGCCAGGAGCGTGTGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGGAGCGG‑A‑CGGCATGGATAATCTTGGCAGAACTGCGGGACGCAGGTCTGGGCg < 2:202157/139‑1 (MQ=255) cTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGGGACTGCAGGAG‑ACC‑TCATCTGCCAGCCTGGCAGAACTGCTGAACGCACGTcccggcc > 2:152078/1‑137 (MQ=255) aTCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGCCAAGCTCCGGCTGCACGGCAGGTGACTCCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCACGCCTGGGCGCtt > 2:373575/1‑143 (MQ=255) ccgtacgTACCGGGTAGGCGTACCACCAGGGTCGTATGCGGCGCCGTGCTGCGGGTGGAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGa < 2:628931/137‑1 (MQ=255) gtccGGCTGCGCGTACCACCAGGATCGTAGGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAg < 2:833611/141‑1 (MQ=255) gTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGACGCTCCGGCTGCACCGCAGCTGACTGCCGAAG‑ACC‑TCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTCCTAATAACGAGtcagc > 2:35801/1‑139 (MQ=255) gTTACGCGTACCACCAGGATCGTAGGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATc < 1:914332/143‑1 (MQ=255) ttACGCGTACCACCAGGATCGTATGCGTCGCCGTGATGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCGGCCAGCCTGGCAGGACAGCTGAACGCAGGGCTGGGCGGCACTGAAAACGAGAAATCg > 2:865385/1‑143 (MQ=255) gcgTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGCCTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTAGACGCAGGTCTGGGCGGATCTGCTGACTAATAATCGTTac > 2:507972/1‑142 (MQ=255) accaGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCa‑a > 1:236391/1‑143 (MQ=255) caGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGCGTAATCGTTAATCCGCA‑CAt > 1:373928/1‑143 (MQ=255) caGGAGCGCATGCGGCGCCGTGCGGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AAt < 2:590805/143‑1 (MQ=255) gATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAAc > 1:101824/1‑143 (MQ=255) gATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCAGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCCGGCAGAACCGCTGCACGCAGGTCTGGGCGGTTCCGATACAGAGTAATCGGTAACCCGCA‑aacaac > 2:401477/1‑143 (MQ=255) cctCGCCCTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACGGCAGAAG‑ACC‑GCATCTGCCAGCCGGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAAcc < 1:699563/141‑1 (MQ=255) gTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGCAG‑ACC‑GCATCTGCCAACCTGGCAGAACTGCTGAACGCAGGCCTGGGCGGTTCTGATAACGCGGCAACGTTAATCCCCA‑AATAAGGTAAAGAccc > 2:321679/1‑143 (MQ=255) gTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATCAGGAGTAATCGTTCATGCGCA‑AATAACGTAAAAAccc > 2:880488/1‑143 (MQ=255) tgctgcGGGTGAAGCACCGGGTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTcgg < 1:101173/143‑1 (MQ=255) gctgcGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTcggc > 1:483874/1‑143 (MQ=255) gcggcGGGTGAAGCTCCGGCGGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCAGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTcggc < 1:483554/140‑1 (MQ=255) gAAGCTCCGGCTGGACGGCGGGTGACGGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGtttttt < 2:101781/143‑1 (MQ=255) tCCGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATgg > 1:153421/1‑143 (MQ=255) ccGGCTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAAGCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATggg > 1:147790/1‑143 (MQ=255) ctcgatcGCAGGTTACTGGAGAAGCC‑C‑GCATCTGCCCGACTGGCATAACTGCTGAACGCAGGTGTGGGCGGTTCTGATAACGAGTAATCGTGAATCCGTA‑AATAACGTGAAAACCCTCTTCGGCGGGTTTTTTTATGGGGGGa < 2:70845/137‑1 (MQ=255) cTGCACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGa < 2:930417/143‑1 (MQ=255) cACCGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGATTTTTTAGGGGGGGAGtt > 1:825290/1‑143 (MQ=255) ccGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTa > 2:762546/1‑143 (MQ=255) ccGCAGGTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTa > 1:151801/1‑143 (MQ=255) gCAGGTGGCTGCATAAG‑ACC‑GCATCTGCCAGGCGGGTCGAACTGCTGAACGCAGGTCTGGGCGGTTTTGATAACGAGTAATCGTTAATTCGGA‑AATAACGTAATAACGCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAgg < 2:182005/143‑1 (MQ=255) gCAGGTGAGTGCAGAAG‑ACC‑GCATCTGCCAGCCGGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAgg < 1:190881/143‑1 (MQ=255) gTGACTGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGATTAGGGaaa > 1:522879/1‑143 (MQ=255) tGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCa > 1:454528/1‑143 (MQ=255) tGCAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGAAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCC‑AATAACGTAAAAACACGCTTCGGCGGGATTTTCTACGGGGGGAGGTTCGGGAAAGGGCa > 2:653954/1‑143 (MQ=255) cAGAAG‑ACC‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCACCATTAC‑TACGAACCTGCTTCGGCGGGTTGTTTTATGGGGGGAGTTTAGGGAAAGAGCAtt > 1:926296/1‑143 (MQ=255) | GCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTGCACCGCAGGTGACTGCAGAAG‑A‑C‑GCATCTGCCAGCCTGGCAGAACTGCTGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCA‑AATAACGTAAAAACCCGCTTCGGCGGGTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATT > NC_000913/4189371‑4189647 |
| 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: ‑ |