Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,354,816 | 0 | A | C | 57.9% | 10.0 / 17.7 | 19 | G26G (GGT→GGG) | rhmD | L‑rhamnonate dehydratase |
Reads supporting (aligned to +/- strand): ref base A (4/4); new base C (1/10); total (5/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.11e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAG > NZ_CP009273/2354569‑2355054 | tGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCCCCCCCTGGTCGTGATAATCAcc > 2:117969/1‑250 (MQ=255) ccTGAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctgc < 2:321147/250‑1 (MQ=255) cAAAAAAACAGCCCATTCCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAAc < 1:266093/250‑1 (MQ=255) tccccGGCTGTCGAACCGGCGAACCCGGTCTGGCCTTTTTCTCTTTCGCCTTCAAAAACCACCGTCCCAAGAACGTTAATGCCAAAGGACGGGCGTGACTGCTCGTAATCGCGGTTTTTCCTCATCGGGGGGGAAATAGGATCGTCTTCCCAGTGATTCGCCCCCTGGTCGGGATAACCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGt < 2:468683/246‑1 (MQ=255) aaaCGGCAAACCCGTTCTGGCCTTTTTCTGCTTCGACTTCAACACCCAACGTGCCAAAAACTTAAATGCAAAATGACTGGCGTGCCTGCTCGAAATCGCGGTATTTCCTCATCGGGGTGGAAATAGGATCGTCTATCAAGTGTTCCGCCCCCTGTTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGtt < 2:308639/250‑1 (MQ=255) gAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCAtt > 1:497163/1‑250 (MQ=255) ccTGGCCGTTTTCTGCTCCGACTTCACCAACCAACGTGCAAGGAACGTTAATGCCAAATGATTGGCGGGACGGCTCGTAACCGCGGTATTTCCGCATCGGGGTGGCAATATGTTCGTTTCCCCAGTGATTCCCCCCCTGGTCTTGATAACCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAg < 2:389250/249‑1 (MQ=255) ccAAAAACGTTAATGCCAAAGGACGGGGGTGCCTGCTCATAATCGCGGTATTTACTCACGGGGGTGGAAATATGATCGTCTATCCAGTGATCCCCCCCCTGGTCTGAAAAATCACCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACg < 2:238533/250‑1 (MQ=255) aaaaaCTTTAATGCCAAATGACTGGGGTGGCTTCTCGAAGTCCGGGTTTTTCCCCGGGGGGGTGTCAATGTGATCTGCTTCCCGGTGATGCGCCCCCTGGTCGTAATAATCCCCGCCGCCAGCGCTTTTTTCTGCGGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTg < 2:50643/247‑1 (MQ=255) tAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATcc > 2:149761/1‑250 (MQ=255) aCTGGCGTTTTTGCTCGCAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCACTGTTTCCCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCATTAAACCAGGCGCGCACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCCAAGGGACTTTTCCATCCCAACAGACGTGAAATACCCCTCGCGCAGGc > 2:86110/1‑249 (MQ=255) cTCGTAATCGCGTGTTTTCCTCATCGGGGTGGCAATATGATCGTCTCTCCGTTGATTCGCCCCCTGGTCGTAATAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCg < 2:88603/250‑1 (MQ=255) ggTATTTACTCATCGGGGGGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGCCGTAATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCa < 1:439694/250‑1 (MQ=255) ccccagcGGGGTGGCAATATGACCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCCCCGCCGCCGGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTAt < 1:373082/243‑1 (MQ=255) gTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGAt < 1:78479/250‑1 (MQ=255) gATCTTCTATCCATTGATCCGCCCCCTGGTCGGAATAACCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCgg < 2:186172/250‑1 (MQ=255) tCCTGTGTTCCGCCCCCGGGTCGTGAAAACCCCCGCCGCCGGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCg < 2:153090/247‑1 (MQ=255) ctgTGATTCCCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCCCCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGaccac < 2:70858/248‑1 (MQ=255) tCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGAACGCAGCGGTCAg > 2:208625/1‑250 (MQ=255) | TGATATCACTGACACATTTACCCTCAATGAAACGGTTAAGATGTTTTTCGACAATAAAACAGCCCATTTCACCGGCTGTCGAAACGGCGAATCCGGTCTGGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATATCCCTGGCGCAGGCAATGGCCTTGCCCGCCAGATAATCACGGTATTCTTCATTGATTTGTAAGCGGGTACCGACCACCGAGATCGCAGCGGTAAG > NZ_CP009273/2354569‑2355054 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |