| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,300,628 | 0 | T | G | 52.9% | 31.1 / 42.6 | 34 | T198P (ACC→CCC) | deaD | ATP‑dependent RNA helicase DeaD |
| Reads supporting (aligned to +/- strand): ref base T (6/10); new base G (18/0); total (24/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.11e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-06 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGTCCAGTAGCTCTGGCTGATGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGTAATGCGACGAATCGCTTCCGGCATGGTTGCAGAGAACAGAGCGGTCTGATGACCTTCCGGGATCTGCGCCATAATGGTTTCAACGTCTT > NZ_CP009273/3300542‑3300717 | aGTCCAGTAGCTCTGGCTGATGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGGAAt > 1:10958/1‑90 (MQ=255) gTCCAGTAGCTCTGGCTGATGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGTAATg > 1:14817/1‑90 (MQ=255) gTCCAGTAGCTCTGGCTGATGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGGaaag > 2:526180/1‑88 (MQ=255) gTCCAGTAGCTCTGGCTGATGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGGaaag > 2:412466/1‑88 (MQ=255) gTCCAGTAGCTCTGGCTGATGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGGAATg > 1:240314/1‑90 (MQ=255) gTCCAGTAGCTCTGGCTGATGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGGAATg > 1:392347/1‑90 (MQ=255) gTCCAGTAGCTCTGGCTGATGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGGAATg > 2:105553/1‑90 (MQ=255) gTAGCTCTGGCTGATGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGGAATGcgacg > 2:259410/1‑90 (MQ=255) tctGGCTGATGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGTAATGGGACGAAGCg > 1:364061/1‑90 (MQ=255) tctGGCTGATGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGTAATGCGACGAATCg > 1:184877/1‑90 (MQ=255) tctGGCTGATGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGGAATGGGACGAATCg > 2:553910/1‑90 (MQ=255) tctGGCTGATGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGGAATGCGACGAATCg > 2:566537/1‑90 (MQ=255) tctGGCTGATGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGGAATGCGACGAATCg > 2:53751/1‑90 (MQ=255) tctGGCTGATGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTCTCATAAAGCGGCGGGGAATGGGAAGGATCg > 1:518686/1‑90 (MQ=255) ggCTGATGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGGAATGCGACGAATCGCtt > 1:609922/1‑90 (MQ=255) tGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGTAATGCGACGAATCGCTTCCGGCa > 2:11495/1‑90 (MQ=255) gTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGTAATGCGACGAATCGCTTCCGGCAt < 2:14817/90‑1 (MQ=255) aGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGGAAAGCGACGAATCGCTTCCGGCATTGt > 2:358087/1‑90 (MQ=255) gggTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGTAATGCGACGAATCGCTTCCGGCATGGTTGCag < 2:485843/90‑1 (MQ=255) gggTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGTAATGCGACGAATCGCTTCCGGCATGGTTGCag < 2:184877/90‑1 (MQ=255) ggTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGGAATGCGACGAATCGCTTCCGGCCTGGTTGCaga > 1:4815/1‑90 (MQ=255) ggTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGGAATGCGACGAAACGCTTCCGGCATGGTTGCaga > 1:93613/1‑90 (MQ=255) ggTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGGAAAGCGACGAATCGCTTCCGGCATGGGTGCagg > 1:358646/1‑89 (MQ=255) gTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGTAATGCGACGAATCGCTTCCGGCATGGTTGCAGAGAAc < 2:494520/90‑1 (MQ=255) cTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGTAATGCGACGAATCGCTTCCGGCATGGTTGCAGAGAACAGAGCg < 1:358087/90‑1 (MQ=255) gAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGGAATGCGACGAATCGCTTCCGGCATGGTTGCAGAGAACAGAGCGGTCtgat > 1:74200/1‑90 (MQ=255) aTGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGTAATGCGACGAATCGCTTCCGGCATGGTTGCAGAGAACAGAGCGGTCtgatga < 2:392347/90‑1 (MQ=255) ttCCTGCGGCTCTTTCATAAAGCGGCGGGGAATGCGACGAATCGCTTCCGGCATGGTTGCAGAGAACAGAGCGGTCTGATGACCTTCCgg > 1:81565/1‑90 (MQ=255) tCATAAAGCGGCGGGTAATGCGACGAATCGCTTCCGGCATGGTTGCAGAGAACAGAGCGGTCTGATGACCTTCCGGGATCTGCGCCATaa < 1:317553/90‑1 (MQ=255) aaaGCGGCGGGTAATGCGACGAATCGCTTCCGGCATGGTTGCag < 2:379886/44‑1 (MQ=255) aaaGCGGCGGGTAATGCGACGAATCGCTTCCGGCATGGTTGCag > 1:379886/1‑44 (MQ=255) gcggcgGGTAATGCGACGAATCGCTTCCGGCATGGTTGCAGAGAACAGAGCg > 1:67087/1‑52 (MQ=255) gcggcgGGTAATGCGACGAATCGCTTCCGGCATGGTTGCAGAGAACAGAGCg < 2:67087/52‑1 (MQ=255) gcggGTAATGCGACGAATCGCTTCCGGCATGGTTGCAGAGAACAGAGCGGTCTGATGACCTTCCGGGATCTGCGCCATAATGGTTTCAAc < 1:53751/90‑1 (MQ=255) tAATGCGACGAATCGCTTCCGGCATGGTTGCAGAGAACAGAGCGGTCTGATGACCTTCCGGGATCTGCGCCATAATGGTTTCAACGTCtt < 2:57035/90‑1 (MQ=255) | AGTCCAGTAGCTCTGGCTGATGTCAGGACGGGTAGTCACGCTGGACTGAATGCGCACTTCCTGCGGCTCTTTCATAAAGCGGCGGGTAATGCGACGAATCGCTTCCGGCATGGTTGCAGAGAACAGAGCGGTCTGATGACCTTCCGGGATCTGCGCCATAATGGTTTCAACGTCTT > NZ_CP009273/3300542‑3300717 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |