Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 702,766:1 | +G | 17.9% | intergenic (+11/‑121) | PP_23SD → / → PP_5SD | 23S ribosomal RNA/5S ribosomal RNA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 702,766 | 1 | . | G | 17.9% | 71.9 / 17.2 | 28 | intergenic (+11/‑121) | PP_23SD/PP_5SD | 23S ribosomal RNA/5S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base . (12/11); new base G (2/3); total (14/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01 |
CTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGC > NC_002947/702631‑702905 | cTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTA‑GCAGATTGctt < 1:159857/150‑3 (MQ=2) gggATCTGGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGGTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAACACTAATTGGCCGTGAGGCTTGACCATATAACACCCCCACAATCTGC‑GTAGCAGATTGcttgg < 2:65807/150‑5 (MQ=11) ttGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAg > 2:398903/1‑150 (MQ=32) tGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGa < 1:398903/150‑1 (MQ=32) ttGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTA‑GCAGATTGct < 2:466191/98‑2 (MQ=11) ttGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTA‑GCAGATTGct > 1:466191/1‑97 (MQ=11) ggTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAata > 2:39947/1‑150 (MQ=32) tAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGAc < 1:281808/100‑1 (MQ=34) tAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGAc > 2:584789/1‑100 (MQ=33) tAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGAc < 1:584789/100‑1 (MQ=34) tAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGAc > 2:281808/1‑100 (MQ=34) aGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCg > 1:40196/1‑150 (MQ=32) gCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAACAATCTGCGTA‑GCAGATTGctt > 1:401231/1‑75 (MQ=11) gCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAACAATCTGCGTA‑GCAGATTGctt < 2:401231/77‑3 (MQ=11) aGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTg > 2:429158/1‑150 (MQ=32) aCCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTg > 1:103030/1‑96 (MQ=34) aCCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTg < 2:103030/96‑1 (MQ=34) ccGTGAGGCTTGACCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCt < 1:540823/137‑1 (MQ=33) ccGTGAGGCTTGACCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCt > 2:540823/1‑137 (MQ=33) tGAGGCTTGACCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACa < 2:40196/150‑1 (MQ=32) gCTTGACCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAAt > 2:235348/1‑150 (MQ=32) aCCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAAt < 2:429069/107‑1 (MQ=33) aCCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCtt < 1:429158/150‑1 (MQ=32) aCCATATAACACACAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAAt > 1:429069/1‑107 (MQ=17) ccATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTg > 2:512825/1‑150 (MQ=32) acCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATa > 1:577929/1‑150 (MQ=32) acCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATa < 1:512825/150‑1 (MQ=32) aaGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGc < 1:235348/150‑1 (MQ=32) | CTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGC‑TACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGC > NC_002947/702631‑702905 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |