Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | CP009974 | 1,920,754 | 1 bp→TTT | intergenic (‑383/+172) | RPPX_08560 ← / ← RPPX_08565 | 16S ribosomal RNA/CoA‑transferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009974 | 1,920,753 | 1 | . | T | 87.5% | 20.0 / 19.1 | 48 | intergenic (‑382/+173) | RPPX_08560/RPPX_08565 | 16S ribosomal RNA/CoA‑transferase |
Reads supporting (aligned to +/- strand): ref base . (1/5); new base T (19/23); total (20/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.79e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01 | |||||||||||
* | CP009974 | 1,920,753 | 2 | . | T | 87.5% | 20.0 / 19.1 | 48 | intergenic (‑382/+173) | RPPX_08560/RPPX_08565 | 16S ribosomal RNA/CoA‑transferase |
Reads supporting (aligned to +/- strand): ref base . (1/5); new base T (19/23); total (20/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.79e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01 | |||||||||||
* | CP009974 | 1,920,754 | 0 | A | T | 100.0% | 24.8 / NA | 48 | intergenic (‑383/+172) | RPPX_08560/RPPX_08565 | 16S ribosomal RNA/CoA‑transferase |
Reads supporting (aligned to +/- strand): ref base A (0/0); new base T (20/28); total (20/28) |
CTTGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTTGATCGTAGCCCTTACAGC‑‑A‑‑CCTTCCGAATCGTCTAACTCATTGAATCTCAAGGAGTTTGTCGTTCCGTTGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAA > CP009974/1920691‑1920867 ||||| cTTGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCgag < 1:199302/151‑1 (MQ=2) cTTGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCgag > 1:3927399/1‑151 (MQ=2) ttGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGCTGTCGCTGGAAGTGGGGCGCACTATAATGGGATTCgaga > 2:1465645/1‑151 (MQ=2) ttGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCgaga < 1:5821978/151‑1 (MQ=2) ttGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCgaga < 2:2855223/151‑1 (MQ=2) ttGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCgaga < 1:4873599/151‑1 (MQ=2) ttGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCgaga < 2:6126431/151‑1 (MQ=2) ttGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGC‑‑TTTCCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCgaga < 2:1913893/151‑1 (MQ=2) ggAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCTTTCCGATGTCGCTGGAAGTGGGGCGCATTATACGGCGCTTCGAGAgg > 2:232995/1‑151 (MQ=2) gAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAggg < 1:296656/151‑1 (MQ=2) gAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAggg < 1:6005577/151‑1 (MQ=2) gAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAggg > 1:1748310/1‑151 (MQ=2) aaGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGc > 1:2643515/1‑151 (MQ=2) aaGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGc > 1:638951/1‑151 (MQ=2) gcgcTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGt < 2:4032767/151‑1 (MQ=2) gcgcTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGt < 2:2918939/151‑1 (MQ=11) cgcTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTc < 2:4428118/151‑1 (MQ=2) cgcTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATGATAAGGGGATTCGAGAGGGCGTc > 2:2514959/1‑151 (MQ=2) cgcTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGC‑‑TTTCCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTc < 1:597109/151‑1 (MQ=2) gcTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCa > 2:1090577/1‑151 (MQ=2) gcTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCa < 2:1201368/151‑1 (MQ=2) cTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCTATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCaa > 2:1739571/1‑151 (MQ=2) cTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCaa > 1:611301/1‑151 (MQ=2) cTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCaa > 2:2071646/1‑151 (MQ=2) tGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCaa > 2:2521693/1‑150 (MQ=2) tGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCaa < 1:2521693/150‑1 (MQ=2) tGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAAc < 1:5317769/151‑1 (MQ=2) tGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAAc > 1:4154449/1‑151 (MQ=2) gTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGC‑‑TTTCCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAAcc < 2:638550/151‑1 (MQ=11) tCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCt > 1:173209/1‑151 (MQ=2) cAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCtt < 2:1673337/151‑1 (MQ=2) aaCCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCttt < 2:425745/151‑1 (MQ=2) aaCCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGC‑‑TTTCCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCttt > 1:5513620/1‑151 (MQ=2) aCCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTa > 1:1515934/1‑151 (MQ=2) ccATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTaa < 1:81832/151‑1 (MQ=2) ccATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTaa < 1:3548520/151‑1 (MQ=2) ccATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTaa < 2:2910332/151‑1 (MQ=2) catcatCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAAt > 2:3542561/1‑151 (MQ=2) catcatCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAAt > 1:3280911/1‑151 (MQ=2) catcatCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAAt > 2:684853/1‑151 (MQ=2) catcatCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAAt > 2:989994/1‑151 (MQ=2) atcatcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAAtt < 1:406450/151‑1 (MQ=2) atcatcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGC‑‑TTTCCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAAtt < 2:2715844/151‑1 (MQ=2) tcatcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAAttt < 2:5443525/151‑1 (MQ=2) atcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCa < 1:5552615/151‑1 (MQ=2) atcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCa < 2:613995/151‑1 (MQ=2) atcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGC‑‑TTTCCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCa < 2:3325031/151‑1 (MQ=2) cAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTC‑ACTACT‑TAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATaa < 2:905783/151‑1 (MQ=11) ||||| CTTGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTTGATCGTAGCCCTTACAGC‑‑A‑‑CCTTCCGAATCGTCTAACTCATTGAATCTCAAGGAGTTTGTCGTTCCGTTGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAA > CP009974/1920691‑1920867 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |