Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,454,629 | 0 | C | A | 33.3% | 19.8 / 14.1 | 25 | R430S (CGC→AGC) | PP_4791 | transposase |
Reads supporting (aligned to +/- strand): ref base C (6/10); new base A (5/3); total (12/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.90e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTGAACAAGGTCTTGAGGAGTTCGAAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGG > NC_002947/5454567‑5454775 | ttGAACAAGGTCTTGAGGAGTTCGAAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCggg < 1:16057/148‑1 (MQ=32) cAAGGTCTTGAGGAGTTCGAAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTAcc > 5:169336/1‑149 (MQ=34) aaGGTCTTGAGGAGTTCGAAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTAccc < 7:26118/149‑1 (MQ=34) ggTCTTGAGGTGTCCTAACCCAGACGCCAACACCGGCAGGCCCGGGCTTTGCTTCGCAAACCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCcgc < 2:363460/149‑1 (MQ=11) aggagTTCGAAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTGGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACg > 7:227217/1‑149 (MQ=38) gTTCGAAAACATAATCCAAGCCCCGCAGGTCCGGGCTTTGGTTCGCAAATCCCGGAAACTGGGGTATCAAATGGTGGCCGCTTTAGTGCTTAGAAAACAATTGGTTTCATTTTTTTTGTTGATATCGGGTTTACCCGCGAATACGGTa < 6:438386/148‑1 (MQ=255) aaCCAGAAGCCAAGACCGGCAGGTCCGGGGTTTGCTTCGCAAAGCCCGGAAACTGGGGTATCAACTGGTGTCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCgg < 2:352808/149‑1 (MQ=255) aaCCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCgg < 1:12086/149‑1 (MQ=255) cAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCGGCAACGGTGaa < 5:25537/149‑1 (MQ=255) aaGACCGGCAGGTCCGGGCATTTCTTGGCAAAGCCCCGAAACTGGGGTATAAACAGGTGGCCCCTAGAGTGGTTATAAAAAAATGGGTTTGATTTTTTTTGTTTATAGCCGGTTTTCCCGCGAAAAAGGTAATCGCGGGCACGGTGAAc > 4:369726/1‑149 (MQ=9) ccGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCGGCAACGGTGAAcggcg < 8:227217/149‑1 (MQ=255) cGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCGGCAACGGTGAAcggcgg < 5:214390/149‑1 (MQ=255) ggCAGGTCCGGGCTTTGCTTAGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCGGCAACGGTGAACGGCggg < 3:369726/149‑1 (MQ=255) gCAGGTCCGGTCTTTGCTTCGCAAAGCCCGTAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCGGCAACGGTGAACGGCGGGt < 5:6536/149‑1 (MQ=255) ccGGGCTTTGCGAAGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCAGCGAATACGGTAGTCGCGGCAACGGTGAACGGCGGGTGGAGAt > 5:393652/1‑149 (MQ=255) gggCTTTGCGAAGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCGGCAACGGTGAACGGCGGGTg > 7:289852/1‑142 (MQ=255) gggCTTTGCGAAGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCGGCAACGGTGAACGGCGGGTg < 8:289852/142‑1 (MQ=255) gggCTTTGCGAAGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCGGCAACGGTGAACGGCGGGTGGAGATTg > 6:258645/1‑149 (MQ=255) gggCTTTGCGAAGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCGGCAACGGTGAACGGCGGGTGGAGATTg > 2:35599/1‑149 (MQ=255) ggCTTTGCTTCGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGTGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCGGCAACGGTGAACGGCGGGTGGAGATTgg > 7:324414/1‑149 (MQ=255) ggCTTTGCGAAGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCGGCAACgg < 4:202449/128‑1 (MQ=255) ggCTTTGCGAAGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCGGCAACgg > 3:202449/1‑128 (MQ=255) cTTCGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCGGCAACGGTGAACGGCGGGTGGAGATTGGCcagcag > 6:39159/1‑149 (MQ=255) cTTCGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCGGCAACGGTGAACGGCGGGTGGAGATTGGCcagcag > 7:298559/1‑149 (MQ=255) ttCGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAgcagg > 5:192830/1‑149 (MQ=255) | TTGAACAAGGTCTTGAGGAGTTCGAAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTATCAACTGGTGGCCGCTTGAGTGCTTAGAAAACAATGGGTTGCATTTTGTTTGATGAGAGCGGGTTTACCCGCGAATACGGTAGTCGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGG > NC_002947/5454567‑5454775 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |