Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,578,446 | 0 | A | C | 64.3% | 0.7 / 11.8 | 14 | Y580D (TAC→GAC) | ggt | gamma‑glutamyltransferase |
Reads supporting (aligned to +/- strand): ref base A (5/0); new base C (0/9); total (5/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.00e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTT > NZ_CP009273/3578363‑3578523 | aGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccgc > 1:51606/1‑90 (MQ=255) tcctttCGGTTTCCCCTAAAAAAATACCCCTCTCCCCAGGAAGGGGGCCGCAAACCTTAGTCCCCCGCCGTTAAATCATCCACCGAgcgc < 2:56514/86‑1 (MQ=255) tccccTAAAAAAATACCCCTCTCCCCACGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc < 2:161371/88‑1 (MQ=255) ccttGAAAAAATCCCCCTCTCCCCACGAAGGGGCCCGCTACCCTTACCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcgc < 2:321934/87‑1 (MQ=255) aaTAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACt > 1:226827/1‑90 (MQ=255) cccTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCg > 1:31600/1‑90 (MQ=255) tccccccAAGAGGGGGCCCGCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc < 1:152136/86‑1 (MQ=255) ctTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg > 1:305359/1‑90 (MQ=255) caggAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa < 2:305359/87‑1 (MQ=255) cacaAAGGGGCCCGCTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa < 2:192363/90‑1 (MQ=255) cACGAGGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa < 2:100500/90‑1 (MQ=255) ggggCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATaa < 2:55128/89‑1 (MQ=255) ggCCCCTACCCTTACCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATg < 1:21503/90‑1 (MQ=255) gCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt > 2:385236/1‑90 (MQ=255) | AGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTT > NZ_CP009273/3578363‑3578523 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |