Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,620,602 | 0 | G | C | 53.8% | 2.5 / 13.4 | 13 | R346G (CGG→GGG) | dcp | peptidyl‑dipeptidase Dcp |
Reads supporting (aligned to +/- strand): ref base G (5/1); new base C (7/0); total (12/1) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.62e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.70e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGA > NZ_CP009273/1620530‑1620670 | cacCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGGTGGGCCTCTTCAAGAACAAATTTTTCCCCCCGTACCTCTTCCGCAt > 1:779040/1‑90 (MQ=255) cacCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTTTCCCCCCGTACCTGTTCCGCAt > 1:19845/1‑90 (MQ=255) cacCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTTTCCCCCCGTACCTGTTCCGCAt > 1:489612/1‑90 (MQ=255) cGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCCTCAAGATCAAATTTTTCCCCCCGTACCTTTTCCGCATAAAATGCCCCGTcccc > 2:629920/1‑89 (MQ=255) gtttttAATTAAAAATATGGTTTGAGCTGCCCCTCAAAAAAATAAATTTCCCCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCAc < 2:433921/87‑1 (MQ=255) gtgtTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCAc > 2:1231935/1‑90 (MQ=255) gtgtTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCCCCGTACCTGTTCCGCATAAAATGCCCAGTCCCAc > 2:1025202/1‑90 (MQ=255) tgtTTAATTCAAAATATGGCTTGGGCTGCGCCTCATCAAGATCAAATTTTTCCCCCCGTACCTGTTCCGCATAAAATTCCCAGTCCCACg > 1:308382/1‑90 (MQ=255) cAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTTTCCCCCCGTACCTGTTCCGCCTAAAAAGCCCCGTCCCCCCGCTgcccgc > 1:718012/1‑90 (MQ=255) tGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAAccc > 2:1197425/1‑90 (MQ=255) gAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCctg > 1:49767/1‑90 (MQ=255) gcgcCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTgctgct > 1:925447/1‑90 (MQ=255) catcaAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGa > 2:796253/1‑90 (MQ=255) | CACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAGATCAAATTTCTCCCGCCGTACCTGTTCCGCATAAAATGCCCAGTCCCACGGCTGCGCGCTAAACCCGCCCTGCTGCTTATCGA > NZ_CP009273/1620530‑1620670 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |