Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,527,926 | 0 | A | C | 57.1% | 19.7 / 35.0 | 28 | G219G (GGT→GGG) | nanM | N‑acetylneuraminate epimerase |
Reads supporting (aligned to +/- strand): ref base A (10/2); new base C (0/16); total (10/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.03e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTTATCACCTTTAT‑T‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGA > NZ_CP009273/4527844‑4528012 | tattaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTTATCACCTTTAt < 2:631292/90‑1 (MQ=255) attaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTTATCACCTTTAt‑t > 1:828441/1‑90 (MQ=255) attaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑GTTTATCACCTTTAt‑t > 1:323390/1‑90 (MQ=255) taCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTTATCACCTTTAT‑T‑ca < 2:681545/90‑1 (MQ=255) gTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTTATCACCTTTAT‑T‑CACAAcc > 1:53035/1‑90 (MQ=255) gTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTTATCACCTTTAT‑T‑CACAAcc > 2:678559/1‑90 (MQ=255) aaTCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTTATCACCTTTAT‑T‑CACAAccgccg > 1:348802/1‑90 (MQ=255) ttaaaaTCCGCCACCCTTTCCAAACCCTGTTTTGG‑CTC‑CCCCTTAAAAAAGCAAGT‑TTTTACCCCTTTTAT‑C‑CACAACCGCCGCACCAGc < 2:389211/87‑1 (MQ=255) tCAAATCCGGCCCCCTTTCGCAAT‑CTGGTTTGG‑GTTCCCCCTTTAAAAAGCCGGT‑TTTTACCCCTTTTAT‑T‑CACAACCGCCGCACCAGcc < 2:525428/90‑1 (MQ=255) aaaaaTACGCCTCCCGTTCGAAACCCGGGTTTGG‑CTG‑CCCATTTAAGAAGCCAGA‑TTTTTTCCCCTTATC‑A‑AAAAACCGCCGCACCAGcc < 2:110419/88‑1 (MQ=255) aTACGCCACCCTTCCGCACCCCTGTTTTGG‑CTT‑CCCCATTAATAACCCAGG‑TTTTATCCCCTTTAT‑C‑CACAACCGCCGCACCAGCCGTTc < 2:92561/90‑1 (MQ=255) aTACGCCACCCGTTC‑CAACTCTTGTTTGGTCTC‑CGCCA‑TAAAAAACCAGGTTTTTACCCCTTTTAT‑C‑CAAAACCGCCGCACCAGCCGTTc < 2:613079/90‑1 (MQ=255) cGGCTTCCGTTCCCAACCCTGTTTTGG‑CTT‑CCCAATTAAAAACCCAGT‑TTTTATCCCTTTTAT‑A‑CACACCCGCCGCACCAGCCGTTCCGt < 2:848712/90‑1 (MQ=255) aTCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTTATCACCTTTAT‑T‑CACAACCGCCGCACCAGCCGTTCCGTAccc > 1:110419/1‑89 (MQ=255) cccTGTTTTGC‑TTT‑CCCCATAAAAAACCCGGT‑TTTTACCCCTTTTAT‑C‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc < 2:554378/89‑1 (MQ=255) ggTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTTATCACCTTTAT‑T‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg > 1:822887/1‑90 (MQ=255) ttGG‑CCT‑CCCCAATAAAAAGCCCGG‑TTTTTCCCCCTTTTTCC‑CACACCCCCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa < 1:678559/90‑1 (MQ=255) gg‑CTT‑CGCCATTAATAAGCCAGG‑TTTTATCACCTTTAT‑T‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt > 2:39852/1‑90 (MQ=255) g‑CTT‑CGCCATTAATAAGCCAGG‑TTTTATCACCTTTAT‑T‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTc > 2:318635/1‑90 (MQ=255) ccATAAAAAAGCCGGT‑TTTTATCCCTTTTTT‑T‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc < 2:348802/90‑1 (MQ=255) ccATAAAAAACCCGGT‑TTTTATCCCTTTTAT‑T‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc < 2:350223/90‑1 (MQ=255) caaTTAACACGCCGGT‑TTTTACCCCTTTTAT‑C‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc < 1:315827/88‑1 (MQ=255) tttAAAAACCCGGT‑TTTTACCCCCTTTAT‑T‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg < 2:1326/89‑1 (MQ=255) aTTAAAAAGCCAGT‑TTTTACCCCTTTTTT‑T‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg < 2:5282/90‑1 (MQ=255) aTTAAAAAGCCAGT‑TTTTACCCCTTTTAT‑C‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg < 2:53035/90‑1 (MQ=255) ataaaaaGCCCGG‑TTTTTACCCCTTTTT‑TCCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg < 1:224668/89‑1 (MQ=255) ttttCTCCCTTTTAT‑C‑CACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa < 1:330925/90‑1 (MQ=255) aTCACCTTTAT‑T‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGa > 2:243562/1‑90 (MQ=255) | TATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTT‑CGCCATTAATAAGCCAGG‑TTTTATCACCTTTAT‑T‑CACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGA > NZ_CP009273/4527844‑4528012 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |