| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,313,327 | C→G | 12.5% | intergenic (‑598/+116) | PP_16SG ← / ← msrQ | 16S ribosomal RNA/methionine sulfoxide reductase heme‑binding subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,313,327 | 0 | C | G | 12.5% | 88.9 / 4.9 | 32 | intergenic (‑598/+116) | PP_16SG/msrQ | 16S ribosomal RNA/methionine sulfoxide reductase heme‑binding subunit |
| Reads supporting (aligned to +/- strand): ref base C (14/14); new base G (2/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 = 9.23e-01 | |||||||||||
GCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGC > NC_002947/5313181‑5313464 | gCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCACCcgcgcg < 2:109916/150‑1 (MQ=255) cTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCTAGCGCCGCCCGCgcgg < 2:17138/150‑1 (MQ=255) aaGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCt < 2:931244/150‑1 (MQ=255) cATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTcaca < 2:97113/150‑1 (MQ=255) tataAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGAtct < 1:168890/143‑1 (MQ=255) tataAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGAtct > 2:168890/1‑143 (MQ=255) taAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGc > 2:1196275/1‑150 (MQ=255) aaGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTcaca < 2:511389/144‑1 (MQ=255) aaGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTcaca > 1:511389/1‑144 (MQ=255) ggggATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCag > 1:540874/1‑150 (MQ=255) agagGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCt < 1:923120/123‑1 (MQ=255) agagGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCt > 2:923120/1‑123 (MQ=255) tCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGcc < 1:653077/150‑1 (MQ=255) tttCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCgg < 1:310599/136‑1 (MQ=255) tttCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCgg > 2:310599/1‑136 (MQ=255) ggggAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGAc > 2:466679/1‑146 (MQ=255) ggggAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGAc < 1:466679/146‑1 (MQ=255) gCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGa > 2:578727/1‑53 (MQ=11) gCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGa < 1:578727/53‑1 (MQ=11) gCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCaaa > 2:1031170/1‑144 (MQ=255) gCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCaaa < 1:1031170/144‑1 (MQ=255) aacaGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCcaca < 2:1044426/89‑1 (MQ=255) aacaGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCcaca > 1:1044426/1‑89 (MQ=255) acaGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCt > 2:1361116/1‑150 (MQ=255) caGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCGCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCGGCCCCACCCCACTCTCAAACCGCCCTACCCTg > 1:129731/1‑150 (MQ=255) gTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCt > 1:1440769/1‑150 (MQ=255) ttCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTAc < 2:1440769/150‑1 (MQ=255) aGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAAACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTg > 1:765511/1‑150 (MQ=255) tGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCtt > 2:1087370/1‑150 (MQ=255) gagaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAAc > 2:445564/1‑51 (MQ=255) gagaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAAc < 1:445564/51‑1 (MQ=255) tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGaa < 1:1196275/150‑1 (MQ=255) gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAgc < 1:1087370/150‑1 (MQ=255) | GCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAAGAAAGTGAAATATAAGGGGAAATACCGTATCGGGGCTGCCAAGTGTTCGCCACAACAGTTTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGACGCCAGAGAAACAAAGGCGGGCCCTCTCCACAGCCATTGCTATCTCAAGACAGGCCCTTGGCAGCCCCACCCCACTCTCAAACCGCCCTACCCTGCTACCTGGCAACCTTGGGAAGC > NC_002947/5313181‑5313464 |
| 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: ‑ |