Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,578,446 | 0 | A | C | 57.1% | 3.0 / 15.5 | 14 | Y580D (TAC→GAC) | ggt | gamma‑glutamyltransferase |
Reads supporting (aligned to +/- strand): ref base A (6/0); new base C (0/8); total (6/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.33e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAC > NZ_CP009273/3578361‑3578512 | aaaGGATAACAAAAAAAAAAGGGCCCCCTTCGGCTTCCCCTGAAAAAATACCCCTTTCCCCACGAAGGGGGCCGCTAACCTTAGcccccc < 1:17210/88‑1 (MQ=255) aTTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGt > 2:103371/1‑90 (MQ=255) ttACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGtt > 1:237347/1‑90 (MQ=255) ttACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGtt > 1:43737/1‑90 (MQ=255) ttACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGtt > 2:81291/1‑90 (MQ=255) aaaaaaaaaaGGCTCCCTTCGCCTTCCCCTAAAAAAATACCCCTCTCCCCAGGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAatc < 2:322277/90‑1 (MQ=255) ggCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAg > 2:29762/1‑90 (MQ=255) tcccTCCGCCTTCCCCTAAAAAAATCCCCCTCTCCCCACGAAGGGGGCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAgcgc < 1:29762/88‑1 (MQ=255) tCGGCTTCCCCTAAAAAAATACCCCTCTCCCCACGAAGGGGGCCGTTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTc < 2:43737/90‑1 (MQ=255) cccTAAAAAAATCCCCCTCCCCCCACGAGGGGGCCCCCTACCTTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcgc < 1:14113/90‑1 (MQ=255) cACGAGGGGGCCCCCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa < 1:264169/90‑1 (MQ=255) cACGAAGAGGGCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa < 2:46426/90‑1 (MQ=255) cGAAGGGGCCCGCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc < 2:136698/90‑1 (MQ=255) gggCCGCTAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcgc < 2:184785/59‑1 (MQ=255) gggCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcgc > 1:184785/1‑59 (MQ=255) | ACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAC > NZ_CP009273/3578361‑3578512 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |