Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,341,237 | 0 | T | G | 54.5% | 34.4 / 42.1 | 33 | V255G (GTG→GGG) | rapZ | RNase adapter RapZ |
Reads supporting (aligned to +/- strand): ref base T (4/11); new base G (18/0); total (22/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.05e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.57e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTTGAGCTATGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCA > NZ_CP009273/3341149‑3341325 | cTTGAGCTATGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGgtg < 1:1112083/90‑1 (MQ=255) cTTGAGCTATGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGgtg < 1:1204168/90‑1 (MQ=255) cTTGAGCTATGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGgtg < 2:658965/90‑1 (MQ=255) gAGCTATGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGGACCGGTCGGGGtat > 2:988976/1‑90 (MQ=255) gCTATGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGGAGGACCCTTCGGGGTatat > 2:1405318/1‑90 (MQ=255) aTGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGGAGGACCCGTCGGGGTATAtttc > 2:968339/1‑88 (MQ=255) aTGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGGAGCACCCGTCCGGGTATAtttc > 2:689460/1‑88 (MQ=255) aTGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGgtg < 2:589757/82‑1 (MQ=255) aTGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGgtg > 1:589757/1‑82 (MQ=255) aTGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGGTTGGGGGGTATTGc > 1:275707/1‑90 (MQ=255) aTGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGGTCCGGGTATAtttc > 1:1597963/1‑88 (MQ=255) ccTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCaa < 2:55913/90‑1 (MQ=255) cTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGGAGCACCGTTCGGGGGTTATTTCAGAGCaaa > 1:373642/1‑89 (MQ=255) aTGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGGAGCACCGTTCGGGGTTTATTGGAGAGCAAATg > 2:753021/1‑90 (MQ=255) aTGCTGGAAACCAACAACCGTAGCTACCTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCCTTCGGGGGATATTGCAGAGCAAATg > 2:1234762/1‑90 (MQ=255) gCTGGACAGCAACCACAGTAGCTACTTGACGGTGGCAATTGGTTGTACCGGCGGGAAGGCACGTTCGCTGTATATTGCAGAGCAACTGGc < 1:476728/90‑1 (MQ=255) gCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGGAGCACCCTTCGGGGTATATTTCAGAGCAACTGGc > 2:1245113/1‑90 (MQ=255) gCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGc < 2:158430/90‑1 (MQ=255) ccAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGGAGCACCCGTCGGGGTATATTTCAGAGCAACTGGCAGactact > 2:715828/1‑90 (MQ=255) ccAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCCGTCGGGGGATATTTGAGAGCAACTGGCAGgctact > 1:404411/1‑90 (MQ=255) cagaCTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTcgcgcg < 1:891841/87‑1 (MQ=255) cTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGAACCCCTCGGGGTATATTGCAGAGCAACTGGCAGGCTACTTCCGCCCGCGCgg > 1:718611/1‑90 (MQ=255) cTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGGGGGTATTTCAGAGCAAATGGCCGACTAATTCCGCTCCCGCGGGaa > 2:1459841/1‑90 (MQ=255) cGCCATTGGTTGTACCGGCGGGAAGCACCCTTCGGGGTATATTTCAGAGCAACTGGCAGACTACTTCCCCTCCCGCGGGAAAAAAGTccc > 2:1588772/1‑89 (MQ=255) ttGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCAc < 1:1405318/90‑1 (MQ=255) ttGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCAc < 1:968339/90‑1 (MQ=255) ggTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGc < 2:1040655/90‑1 (MQ=255) tACCGGCGGGGAGCACCGCTCGGTGGATATTTCAGAGGAAATGGCCGACTACTTCCCCCCGCGCGGGAAAAACCACCAGTCACGCCATCg > 1:1397364/1‑90 (MQ=255) tACCGGCGGGAAGCACCGTTCGGGGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGGAAAAAAGTCCCGTCACGCCATCg > 2:459840/1‑90 (MQ=255) ccGGCGGGAAGCACCGTTCGGGGTATATTGCAGAGCAACTGGCAGACTACTTCCCCTCGCGCGGGAAAAAAGGCCAGTCACGCCATCGTa > 2:1292443/1‑90 (MQ=255) ggcggGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACg < 1:1245113/90‑1 (MQ=255) gcggGAAGCACCGTTCGGGGTATATTTCAGAGCAACTTGCAGACCACTTCCGCTCCCGCGGGAAAAAAATCCAGTCACGCCATCGTACGc > 1:1089854/1‑90 (MQ=255) aGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGaaa < 2:1077862/90‑1 (MQ=255) aGCACCATTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGAAAAAACGTCCAGTCACGCCATCGTACGCTGGaaa > 2:615336/1‑90 (MQ=255) cGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGt < 1:689460/90‑1 (MQ=255) gTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTa > 2:360246/1‑90 (MQ=255) gTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTa > 2:660955/1‑90 (MQ=255) ttCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTaa < 2:373642/90‑1 (MQ=255) gtgtATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCa > 2:684439/1‑90 (MQ=255) | CTTGAGCTATGGTTACCTATGCTGGAAACCAACAACCGTAGCTACTTGACGGTCGCCATTGGTTGTACCGGCGGGAAGCACCGTTCGGTGTATATTGCAGAGCAACTGGCAGACTACTTCCGCTCGCGCGGTAAAAACGTCCAGTCACGCCATCGTACGCTGGAAAAACGTAAACCA > NZ_CP009273/3341149‑3341325 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |