| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,578,446 | 0 | A | C | 51.8% | 18.6 / 37.7 | 27 | Y580D (TAC→GAC) | ggt | gamma‑glutamyltransferase |
| Reads supporting (aligned to +/- strand): ref base A (13/0); new base C (0/14); total (13/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.99e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-06 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTT > NZ_CP009273/3578360‑3578523 | tACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAccc > 1:871869/1‑90 (MQ=255)tACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAccc > 1:587194/1‑90 (MQ=255) aGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccgc > 2:727640/1‑90 (MQ=255) ttACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGtt > 2:586834/1‑90 (MQ=255) ttACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGtt > 2:491513/1‑90 (MQ=255) tagaaaaaaaaaAGGGCTCCCTTCGGTTTGCCCGAAAAAAATACCCCTCTTCCCACGAAGGGGGCCGCTACCCTTAGTCCCCCGCCGTTa < 2:51616/86‑1 (MQ=255) caaaaaaaaaaaGGCTCCCTTCGGCTTGCCCTAAAAAAATACCCCTCTCCCCAAAAAGGGGCCCGCAACCTTTAGTCCCCCGCCATTaaa < 2:500102/88‑1 (MQ=255) aaaaaaaaaaaaCCCTCCCTCCGGCTTCCCCACCAAAAATCCCCCCCTCCCCACGAGGGGGCCCGCTAACCTTAGCCCCCCGCCGTTaaa < 1:465666/88‑1 (MQ=255) aaaaaaaaaaAGGGCCCCCTGCGCTTTCCCTTGAAAAAACACCCCTCCCCCCCGGAAGAGGGCCGATATCCTTAACCCCCCGCCGTTaaa < 1:629603/88‑1 (MQ=255) aTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAatcat > 2:509644/1‑90 (MQ=255) aaaaaagGGCTCCCTTCGGTTTGCCCAAAAAAAATCCCCCTCTCCCCACGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAatcatc < 1:98251/89‑1 (MQ=255) aaaaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCa > 1:51616/1‑90 (MQ=255) tcccTCCGGCTTCCCCTGAAAAAATACCCCTCTCCCCACGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAgcgc < 2:231876/88‑1 (MQ=255) gCTTCCCCTAAAAAAATCCCCCTCCCCCCACGAAGGGGCCCGCTACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGa < 1:186094/90‑1 (MQ=255) tccccTAAAAAAATACCCCTCTCCCCACGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc < 2:423131/88‑1 (MQ=255) tGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc > 1:187046/1‑90 (MQ=255) gccccTAACAAAACCCCCCTCCCCCCACGAGGGGGCCCGCAACCTTTACCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgc < 1:586834/88‑1 (MQ=255) cccTAAAAAAATACCCCTCTCCCCACGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcgc < 2:404408/90‑1 (MQ=255) gACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTa > 1:404408/1‑90 (MQ=255) aaaaTCGCCCTCTCCCCACGAGGGGGCCCCCAACCCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACaa < 1:509644/90‑1 (MQ=255) tctTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTcc > 1:569629/1‑90 (MQ=255) ctTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg > 1:819035/1‑90 (MQ=255) cACGAAGGGGCCCGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa < 2:386325/90‑1 (MQ=255) cACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCa > 1:165257/1‑90 (MQ=255) agagGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCAt > 1:385854/1‑90 (MQ=255) gggcccGCTACCCTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATaa < 2:654867/86‑1 (MQ=255) ccTAACCTTAGCCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt < 2:518292/89‑1 (MQ=255) | TACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTT > NZ_CP009273/3578360‑3578523 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |