Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 899,438 | 0 | A | C | 73.7% | ‑2.5 / 10.8 | 19 | G162G (GGT→GGG) | ybjP | lipoprotein |
Reads supporting (aligned to +/- strand): ref base A (5/0); new base C (0/14); total (5/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.60e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTG > NZ_CP009273/899356‑899527 | aTTCCTCACAAAATGCCGGAAATTTTGGGGGTTTTTTTGTTTGATAAGGGCGTTAGCGTTTCTCGATGGACTGAGGGAGCGTCCccgccg < 2:8500/90‑1 (MQ=255) tCACAAAATGCCGGAAATTTTAGGGGGTTTATTGGTTGATCAAGGCGTTAGCGTTTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCa < 1:410769/90‑1 (MQ=255) tGCCGGACATTTTAGGGGTTTTATTGGTTGATCAAGGGGTTAGCGTTTCTCGATGGACGGAGGGGGCGTCCCCGCCGTGGCATGAACGCt < 2:606354/90‑1 (MQ=255) gCCGGACATTTTCGGGGTTTTATTGGTTGATCAAGGCTTTAGCGATTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTa < 1:486651/90‑1 (MQ=255) ttATTGGTTAATAAGGGCGTTAGCTATTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACAt < 2:232367/90‑1 (MQ=255) ttGATAAAGGCGTTAGCGATTCTCGAGGGACTGAGGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGAt < 1:488928/90‑1 (MQ=255) ttGACCAGGGCGTTAGCGTTTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGAt < 2:159717/90‑1 (MQ=255) tGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATg > 1:550318/1‑90 (MQ=255) aaaGGCGTTAGCGTTTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGAccc < 1:269752/89‑1 (MQ=255) ttAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTg > 2:785426/1‑90 (MQ=255) ttAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTg > 2:173410/1‑90 (MQ=255) gCGTTTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGcc < 1:292861/90‑1 (MQ=255) cttttCCCGAGGGACGGAGGGGGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCt < 2:159511/87‑1 (MQ=255) aTTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTc > 1:278144/1‑90 (MQ=255) ctcGATGGACTGCCGGAGCCCCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTg < 1:142874/90‑1 (MQ=255) cGATGGACGGACGGAGCGCCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGaa < 2:581737/90‑1 (MQ=255) gACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAaca > 2:406461/1‑90 (MQ=255) ggAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACtt < 1:24458/90‑1 (MQ=255) cGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATc < 1:596944/90‑1 (MQ=255) aCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTg > 1:470588/1‑90 (MQ=255) | ATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTG > NZ_CP009273/899356‑899527 |
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: ‑ |