| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | 1 | 593,572 | 0 | A | G | 48.5% | ‑0.6 / 31.0 | 33 | intergenic (+96019/‑105771) | unknown/unknown | tRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA |
| Reads supporting (aligned to +/- strand): ref base A (9/8); new base G (6/10); total (15/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.91e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.34e-01 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGCTCTGCACCCCAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGA > 1/593349‑593712 | cacCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAGGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGtt > 1:275174/1‑251 (MQ=21)cacCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGtt > 1:215795/1‑251 (MQ=21) gCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTAt > 1:341581/1‑251 (MQ=21) aaaTGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGa < 1:56730/251‑1 (MQ=18) ggCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGcc < 1:130955/249‑1 (MQ=18) ggCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCtt > 1:858432/1‑251 (MQ=18) cacaCCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGcc > 1:68607/1‑243 (MQ=18) acacCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGaa < 1:580042/251‑1 (MQ=18) tCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAAt > 1:26391/1‑251 (MQ=17) cTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAcc < 1:254574/251‑1 (MQ=14) cTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAcc < 1:427356/251‑1 (MQ=14) cgctttccgaaacaaaaaagtgatctGCACCCCAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAAt < 1:768512/228‑1 (MQ=11) cgctttccgaaacaaaaaagtgatctGCACCCCAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAAt < 1:456693/228‑1 (MQ=11) ctttccgaaacaaaaaagtgatctGCACCCCAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAg < 1:731158/230‑1 (MQ=11) ctttccgaaacaaaaaagtgatctGCACCCCAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAg < 1:327616/230‑1 (MQ=11) aatttcaaCCGCACTTTGCTCTGCACCCCAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTg < 1:688941/244‑1 (MQ=2) aatttcaaCCGCACTTTGCTCTGCACCCCAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTg < 1:336185/244‑1 (MQ=2) ttCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAAt < 1:828780/251‑1 (MQ=11) ttCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAAt < 1:185103/251‑1 (MQ=11) acaaaaaagtgatctGCACCCCAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAAc < 1:73785/218‑1 (MQ=11) gCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat > 1:503364/1‑222 (MQ=11) cACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAccc < 1:9422/227‑1 (MQ=1) cACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAg > 1:516002/1‑244 (MQ=1) cACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACAAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCAGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGCAAGAGGGCTTAAt > 1:190524/1‑201 (MQ=255) ttGCTCTGCACCCCAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATgg < 1:325746/193‑1 (MQ=255) gtgatctGCACCCCAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAAt > 1:402714/5‑191 (MQ=11) gtgatctGCACCCCAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGaaa > 1:846203/5‑251 (MQ=9) gCTCTGCACCCCAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAg > 1:256004/1‑232 (MQ=255) gCACCCCAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAAt > 1:210998/1‑225 (MQ=255) ccccAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATgg > 1:216402/1‑183 (MQ=255) ccccAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat < 1:205035/201‑1 (MQ=9) ccccAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGa > 1:684359/1‑251 (MQ=9) aaaaCCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAGAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTGAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGa < 1:511197/251‑1 (MQ=9) | CACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGCTCTGCACCCCAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGA > 1/593349‑593712 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |