Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,578,446 | 0 | A | C | 59.1% | 3.4 / 24.0 | 22 | Y580D (TAC→GAC) | ggt | gamma‑glutamyltransferase |
Reads supporting (aligned to +/- strand): ref base A (8/1); new base C (0/13); total (8/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.81e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAA‑CCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTT > NZ_CP009273/3578372‑3578522 | aaaTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAA‑CCTTAGTACCCCGCCGTTAAatc > 1:499851/1‑90 (MQ=255) aaaaaaagGGCTCCCTTCGCTTTCCCCTGAAAAAATACCCCTCTCCCCACGAAGGGGGCCGCTAC‑CCTTAGTCCCCCGCCGTTAAatcat < 2:388758/88‑1 (MQ=255) aTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAA‑CCTTAGTACCCCGCCGTTAAatcat > 2:111575/1‑90 (MQ=255) aaaGGCTCCCTTCGGCTTGCCCTGAAAAAATACCCCTCTCCCCACGAAGGGGCCCGCTAC‑CCTTAGTCCCCCGCCGTTAAATCATCCAcc < 2:110036/90‑1 (MQ=255) ggCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAA‑CCTTAGTACCCCGCCGTTAAATCATCCACCGAg > 1:666207/1‑90 (MQ=255) cTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAA‑CCTTAGTACCCCGCCGTTAAATCATCCACCGAgcg > 1:282082/1‑90 (MQ=255) tACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAA‑CCTTAGTACCCCGCCGTTAAATCATCCACCGAgcgc > 2:352715/1‑90 (MQ=255) aCCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAA‑CCTTAGTACCCCGCCGTTAAATCATCCACCGAGcgcg > 1:110036/1‑90 (MQ=255) ggCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAA‑CCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCgg > 1:100553/1‑90 (MQ=255) tccccTAAAAAAATCCCCCTCTCCCCAAGAGGGGGCCCCCAAC‑CCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc < 1:111575/88‑1 (MQ=255) tGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAA‑CCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc > 2:178271/1‑90 (MQ=255) aaaaaTACCCCTATCCCCCGGAAGGGGGCCGCTAC‑CCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACa < 2:593528/89‑1 (MQ=255) atcccccTCTCCCCAGGAAGAGGCCCGCTAC‑CCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTc < 2:100553/86‑1 (MQ=255) tcccccTCTCCCCACGAAGGGGGCCGCTAC‑CCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCa < 1:352715/87‑1 (MQ=255) cttccccACGAGGGGGGCCCTAAA‑CTTTAGCCCCCCGCCGTTAAAACATCCTCCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc < 2:908441/86‑1 (MQ=255) caggAAGGGGGCCGCAAC‑CTTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa < 2:78864/87‑1 (MQ=255) caaaaaGGGGCCCGCTAC‑CCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa < 2:157882/86‑1 (MQ=255) cACGAGGGGGCCCGCTAC‑CCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa < 2:525740/90‑1 (MQ=255) cccaaGGGGGCCCCTTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa < 1:615676/87‑1 (MQ=255) cGAAGAGGGCCGCTAA‑CCTTAGTACCCGGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAc < 2:253048/90‑1 (MQ=255) ggggCCGCAAC‑CCTTACCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATaa < 1:166848/89‑1 (MQ=255) cGCTAC‑CTTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt < 2:920849/90‑1 (MQ=255) | AAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAA‑CCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTT > NZ_CP009273/3578372‑3578522 |
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: ‑ |