Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,578,446 | 0 | A | C | 61.5% | ‑0.7 / 10.0 | 13 | Y580D (TAC→GAC) | ggt | gamma‑glutamyltransferase |
Reads supporting (aligned to +/- strand): ref base A (5/0); new base C (0/8); total (5/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.77e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCC > NZ_CP009273/3578360‑3578533 | tACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAccc > 1:331644/1‑90 (MQ=255) tACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAccc > 2:632583/1‑90 (MQ=255) aCAGGATTACAAAAAAAAAAAGGCTCCCTTCGGCTTGCCCTAAAAAAATACCCCTCTCCCCACGAAGGGGGCCGCTAACCTTAGTCcccc < 2:203040/90‑1 (MQ=255) agagAAAAAAAAAAGGCTCCCTTCGGCTTCCCCTGAAAAAATACCCCTCTTCCCACGAAGGGGGCCGCTAACTTTAGTCCCCCGCCGTTa < 1:629630/87‑1 (MQ=255) aaaaGGCCCCCTCCGCCTTCCCCTAAAAAAACCCCCCTCTCCCCACGAAGGGGGCCGCTACCCTTACTCCCCCGCCGTTAAATCATCCAc < 2:331644/90‑1 (MQ=255) cTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAgcg > 2:364937/1‑90 (MQ=255) caggAAGGGGGCCGCAACCCTAAGACCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa < 1:306379/87‑1 (MQ=255) cacaAAGAGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa < 1:49815/90‑1 (MQ=255) cGAAGGGGGCCGCTACCCTAATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc < 1:538334/90‑1 (MQ=255) gggCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAAt > 1:18346/1‑90 (MQ=255) ctccccTTTGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCCCCGTCCGGCCCAACCATAATGCTTtg < 1:632583/86‑1 (MQ=255) cTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtg < 1:580022/90‑1 (MQ=255) gTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGcc > 2:450771/1‑90 (MQ=255) | TACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCC > NZ_CP009273/3578360‑3578533 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |