Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 360,762 | 0 | T | G | 62.5% | 10.8 / 16.4 | 16 | T328P (ACC→CCC) | lacI | DNA‑binding transcriptional repressor LacI |
Reads supporting (aligned to +/- strand): ref base T (1/5); new base G (10/0); total (11/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.37e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATAT > NZ_CP009273/360635‑360899 | aTGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTGTGTGGATTGGGGACCAGGGGGGGTTTTCTTTTCa > 1:238272/1‑139 (MQ=255) gTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTGGGGACCAGGGGGGGTTTTCTTTTTACCAGt > 1:264785/1‑139 (MQ=255) gTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGGATTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGt > 2:121342/1‑139 (MQ=255) cgcgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGTCCAACGCGCGGGGAGAGGCGGTTTTCGTTTTGGGCACCCGGGTGGGTTTTCTTTTCACCCGCGGGACGGGCaa > 1:190316/1‑139 (MQ=255) cgcgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGGACCAGGGGGGGTTTTTTTTTTACCAGGGAGACGGGcca > 1:183228/1‑137 (MQ=255) cgcgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGCACCAGGGTGGGTTTTCTTTTTACCAGTGAGACGGGcc > 1:491610/1‑137 (MQ=255) cgcgTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTGGGGACCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCaa > 1:442293/1‑139 (MQ=255) tCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGcc < 2:467706/138‑1 (MQ=255) ccAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGGACCAGGGGGGGTTTTTTTTTCACCCGTGAGACGGGCAACAGCTGATTGGCCTTCACCCCCCGGCCCTGAGAGAGGTgcagca > 1:326210/1‑138 (MQ=255) cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCCCCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAg < 2:342338/139‑1 (MQ=255) tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGATTGGGGACCAGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCGACAGCTGGTTTCCCTTCCCCCCCCGGGCCTGAGAGAGGTGCAGCAAGCGGt > 2:25811/1‑139 (MQ=255) tCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGATTTGGGACCCGGGGGGGTTTTCTTTTTCCCCGTGAGACGGGCAACAGCTGATTTCCCTTCACCCCCTGGCCCTGTGAGAGGTGCAGCAAGCGGTCCACGCTGGgtt > 2:437808/1‑139 (MQ=255) gCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGa < 1:196210/138‑1 (MQ=255) aCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAg < 1:298219/84‑1 (MQ=255) aCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAg > 2:298219/1‑84 (MQ=255) ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg < 2:183228/139‑1 (MQ=255) gtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTTAGAGAGTTGCAGCAAGCGGTCCGCGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat < 1:181213/139‑1 (MQ=255) gtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGata > 2:440766/1‑138 (MQ=255) | ATGCGGTCCTCAATCGCGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTGGGCACCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATAT > NZ_CP009273/360635‑360899 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |