Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,527,926 | 0 | A | C | 61.1% | 6.4 / 18.7 | 18 | G219G (GGT→GGG) | nanM | N‑acetylneuraminate epimerase |
Reads supporting (aligned to +/- strand): ref base A (7/0); new base C (0/11); total (7/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.14e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGA > NZ_CP009273/4527845‑4528012 | attaCCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTAtt > 2:397499/1‑90 (MQ=255) ttCAAATACGGCACCCGTTCGCAATCCTGTTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCCCC‑TTTATTCACAACCGCCGCACCAGc < 2:197130/90‑1 (MQ=255) tcaaaaTCCGGCACCCTTCCCAAACCCTGTTTTGG‑CTTCGCCATAAAAAAGCCGGT‑TTTTATCCCT‑TTTATTCACAACCGCCGCACCAGc < 2:647677/87‑1 (MQ=255) aTCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCa > 1:548750/1‑90 (MQ=255) cctttCGCAACCCTGGTTGGG‑TTTCGCCATTAATAAGCCGGT‑TTTTATCCCC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAgg < 2:162575/87‑1 (MQ=255) tCCTGGTTTGG‑CTTCGCCATTAAAAAGCCAGG‑TTTTATCCCT‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc < 2:292875/90‑1 (MQ=255) ccTGGTTTGGTCTCCCCCATAAAAAAGCCGGT‑TTTTACCCCT‑TTTTTTAACACCCGCCGCACCAGCCGTTCTGTACCAGGGCGATTCGcc < 2:327924/90‑1 (MQ=255) cTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAg > 2:302466/1‑90 (MQ=255) ttggg‑CTCCCCCATTAAAAAGCCAGG‑TTTTATCCCT‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa < 2:180386/87‑1 (MQ=255) ttGG‑GTTCCCCATTAAAAAACCCGGTTTTTATCCCT‑TTTATTCACACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa < 2:433229/90‑1 (MQ=255) ttGG‑CTTCCCCAATAAAAAACCCGG‑TTTTTCCCCCTTTTATCCCCAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTa < 2:618773/90‑1 (MQ=255) gg‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt > 1:292875/1‑90 (MQ=255) ccATTAATAAGCCAGT‑TTTTATCCCT‑TTTTTTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc < 2:135169/90‑1 (MQ=255) aaacccAGT‑TTTTACCCCT‑TTTATCCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTg < 2:217213/86‑1 (MQ=255) aGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAg > 1:108550/1‑90 (MQ=255) gCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAgg > 1:503135/1‑90 (MQ=255) ttttATCCCT‑TTTTTTCAAACCCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCaa < 2:700226/90‑1 (MQ=255) aTCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGa > 2:186996/1‑90 (MQ=255) | ATTACCGGTGAAATCAAGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGG‑CTTCGCCATTAATAAGCCAGG‑TTTTATCACC‑TTTATTCACAACCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATCAAAAGA > NZ_CP009273/4527845‑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: ‑ |