Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,453,630 | 0 | C | G | 60.0% | 1.2 / 10.2 | 10 | A107G (GCG→GGG) | paaC | 3‑hydroxyacyl‑CoA dehydrogenase PaaC |
Reads supporting (aligned to +/- strand): ref base C (4/0); new base G (0/6); total (4/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.76e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTGGCAGCTGCGGACCTGGTCATTGAAGCGGCGTCTGAACGTCTGGAAGTCAAAAAAG‑CGCTCTTTGCACAGCTGGCGGAAGTTTGCCCGCCACAAACGCTATTGACCACTAACACTTCGTCAATCTCTATAA > NZ_CP009273/1453553‑1453686 | gcTGGCAGCTGCGGACCTGGTCATTGAAGCGGCGTCTGAACGTCTGGAAGTCAAAAAAG‑CGCTCTTTGCACAGCTGGCGGAAGTTTGccc > 1:134873/1‑90 (MQ=255) tgttCATGGAGGCGGCGTCTGACCGTC‑GGAAAAAAAAAAAGCCGCTCTTCCCACGGGGGGGGGATTTTTCCCCGCCACAAACGCTATTGa < 2:341063/87‑1 (MQ=255) tGGTCATTGAAGCGGCGTCTGAACGTCTGGAAGTCAAAAAAG‑CGCTCTTTGCACAGCTGGCGGAAGTTTGCCCGCCACAAACGCTATTGa > 1:386004/1‑90 (MQ=255) tCATTGAAGCGGCGTCTGAACGTCTGGAAGTCAAAAAAG‑CGCTCTTTGCACAGCTGGCGGAAGTTTGCCCGCCACAAACGCTATTGACCa > 2:103546/1‑90 (MQ=255) cggcgTCGGAACGTCAGAAAAAAAAAAAAG‑CCCTTTTTGCACGGGGGGGGGATTTTTCCCCGCCACAAACGCTATTGACCACTAACACtt < 2:134873/90‑1 (MQ=255) ctcaaCGTGTGAAAGAAAAAAACG‑CCCTCTTCCCACGGGTGGGGGAATTTTGCCCGCCACAAACGCTATTGACCACTAACACTTCGTCaa < 1:49932/87‑1 (MQ=255) cTGAACGTCTGAAAAAAAAAAAGG‑CTCTTTTTGCCCGGGGGGGGGATGTTTGCCCGCCACAAACGCTATTGACCACTAACACTTCGTCaa < 1:276249/90‑1 (MQ=255) cTGAACGTCGGAAAGAAAAAAAAG‑CCCTTTTTGCACGGCGGGGGGAATTTTGCCCGCCACAAACGCTATTGACCACTAACACTTCGTCaa < 1:164452/90‑1 (MQ=255) gTCTGGAAGTCAAAAAAG‑CGCTCTTTGCACAGCTGGCGGAAGTTTGCCCGCCACAAACGCTATTGACCACTAACACTTCGTCAATCTCta > 1:196687/1‑90 (MQ=255) gggAAGAAAAAAAAG‑CGCTCTTTCCACAGCGGGGGGAAGTTTGCCCGCCACAAACGCTATTGACCACTAACACTTCGTCAATCTCTATaa < 1:103546/89‑1 (MQ=255) | GCTGGCAGCTGCGGACCTGGTCATTGAAGCGGCGTCTGAACGTCTGGAAGTCAAAAAAG‑CGCTCTTTGCACAGCTGGCGGAAGTTTGCCCGCCACAAACGCTATTGACCACTAACACTTCGTCAATCTCTATAA > NZ_CP009273/1453553‑1453686 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |