| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 535,822 | A→C | 15.6% | E85A (GAA→GCA) ‡ | rplA → | 50S ribosomal protein L1 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 535,822 | 0 | A | C | 15.6% | 67.6 / 6.5 | 32 | E85A (GAA→GCA) ‡ | rplA | 50S ribosomal protein L1 |
| Reads supporting (aligned to +/- strand): ref base A (13/14); new base C (3/2); total (16/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GTTCGTCGAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTGG > NC_002947/535679‑535962 | gTTCGTCGAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGaa < 2:417115/145‑1 (MQ=255)gTTCGTCGAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGaa > 1:417115/1‑145 (MQ=255)gTTCGTCGAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAgctgc > 2:311480/1‑150 (MQ=255) ttCGTCGAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAgctgct > 2:97360/1‑150 (MQ=255) tcgAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTgg < 2:10116/150‑1 (MQ=255) aGTCCTAAGACATCGCCGTTAACCGCGGTGTTGATGCGCGTAAATCCGACCAGGTCGTACGTCGCGCTACCGTGCTGCCACACGGCGGTGGCAAGGCCGTGCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTTCAGCTGctcgagcct < 2:133966/150‑7 (MQ=255) tCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGCCAAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGcc < 2:744271/150‑1 (MQ=255) tACGACATGGCCCTTAACCTCGGTGTTGACCCGCGTAAATCCAACCAGGTCGTACGTAGCGCTACCGTGCTGCCCCACGGCACTGGCAAGACCGTTCGCGTTCCTGTCTTCACCCAGAGTCCAGCAGCTGAACCTCCTCTGACTGCCgga > 2:19377/1‑149 (MQ=255) cGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGgcgc < 1:311480/150‑1 (MQ=255) tCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATc < 2:713391/150‑1 (MQ=255) ttAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAg < 2:682523/150‑1 (MQ=255) tGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCt < 2:487605/150‑1 (MQ=255) gACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCt < 2:774451/134‑1 (MQ=255) gACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCt > 1:774451/1‑134 (MQ=255) aCCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAg > 1:466969/1‑150 (MQ=255) gtacgtaGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGa > 1:35000/1‑119 (MQ=255) gtacgtaGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGa < 2:35000/119‑1 (MQ=255) gtacgtaGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATc > 1:332006/1‑124 (MQ=255) gtacgtaGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATc < 2:332006/124‑1 (MQ=255) gcgcTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACt > 1:648958/1‑150 (MQ=255) gcTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGACGATCCTCTGGCTGCCGAATAAGAGCGCGTTTGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTAt < 1:97360/150‑1 (MQ=255) gTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGAcgtc < 1:273137/150‑1 (MQ=255) ccACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAACGAAAGGCGGCGACCTGAACTATGACGTCGTTatc > 2:731879/1‑150 (MQ=255) cAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGc > 1:379048/1‑150 (MQ=255) gCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGt < 1:768442/150‑1 (MQ=255) cTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTc > 1:386253/1‑150 (MQ=255) tCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGc > 1:586246/1‑150 (MQ=255) tCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTgg < 2:682486/150‑1 (MQ=255) aGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGCAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAgg > 2:462230/1‑150 (MQ=255) gccagagcagcTTCAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGtt > 1:160293/6‑150 (MQ=255) ccagagcagcTTCAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGaa > 1:760636/5‑88 (MQ=255) ccagagcagcTTCAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGaa < 2:760636/84‑1 (MQ=255) agagcagcTTCAGCTGCCCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGCTAGGTTCTg > 2:694988/3‑150 (MQ=255) cagcagCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTgg < 1:15142/150‑1 (MQ=255) | GTTCGTCGAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTGG > NC_002947/535679‑535962 |
| 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: ‑ |