Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,193,661 | 0 | T | G | 53.3% | 14.8 / 16.7 | 15 | pseudogene (3709/3796 nt) | BW25113_RS25680 | WGR and DUF4132 domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (3/4); new base G (8/0); total (11/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.56e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAG > NZ_CP009273/2193572‑2193749 | ggCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGt < 1:292369/90‑1 (MQ=255) aTGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGGGGAAt > 2:71265/1‑90 (MQ=255) aGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAAcaca < 2:58771/90‑1 (MQ=255) gcTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTtc < 1:84364/90‑1 (MQ=255) gAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGGGt > 2:304033/1‑63 (MQ=255) gAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGGGTAATTTAACACAGGAAGAGCAGTTTTcc > 2:579443/1‑90 (MQ=255) gAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGt < 1:304033/63‑1 (MQ=255) gAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTcc < 1:21478/90‑1 (MQ=255) gcTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGGAATTCAAAACAGGAAGCGCAGTTCTCCGGAAt > 2:127883/1‑90 (MQ=255) tAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGGGGAATTTAAAACAGGAAGCGCAGTTCTCCGTAATCGATg > 2:37595/1‑90 (MQ=255) tAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGGGGAATTTAAAACAGGAAGCGCAGGTCTCCGTAATCGATg > 2:21478/1‑90 (MQ=255) aagcTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGGGGAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATc > 2:180834/1‑90 (MQ=255) aagcTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGGGGAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATc > 2:341762/1‑90 (MQ=255) agcTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGGAATTCAAAACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCa > 1:529423/1‑90 (MQ=255) aGGCAAAGCAGAAAGATATGGCTGGGGGAGGAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAAt > 2:56677/1‑90 (MQ=255) cTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGtt < 2:529423/90‑1 (MQ=255) gTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAg > 1:583671/1‑90 (MQ=255) gTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAg > 1:70983/1‑90 (MQ=255) | GGCATGTCGCCAGCCGAACTCAGCGCTGAGCAGCTCTTAAGCAAGCTGTGGCTATGGGAAGGCAAAGCAGAAAGATATGGCTGGGGGAGTAATTCAACACAGGAAGCGCAGTTCTCCGTAATCGATGCCATCACCGCCAGCGAGCTAATTAACGATATTGAAGCCCTGTTTGAATAAG > NZ_CP009273/2193572‑2193749 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |