Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 294,664 | 0 | A | C | 56.2% | 12.8 / 19.1 | 16 | V661G (GTG→GGG) | paoC | aldehyde oxidoreductase molybdenum‑binding subunit PaoC |
Reads supporting (aligned to +/- strand): ref base A (7/0); new base C (0/9); total (7/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.74e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCG > NZ_CP009273/294577‑294739 | acacCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCAcc > 2:294747/1‑90 (MQ=255) tCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCa > 1:275657/1‑90 (MQ=255) ggggAGGATGGGGGGTCGTTATCACCCAGAAAAACCCCCTCTTTTTTTGGGATATCGGGATGACCCGCCCCCTCATCCCCCGCCATATCg < 1:76736/90‑1 (MQ=255) gggAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCa < 2:523127/69‑1 (MQ=255) gggAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCa > 1:523127/1‑69 (MQ=255) gaTATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTa > 1:160430/1‑90 (MQ=255) gaTATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTa > 1:498522/1‑90 (MQ=255) tggggggTGGGTATCACCAAGAAAAACCCCCTCCTTTTTTGGGATATCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAAc < 1:107984/87‑1 (MQ=255) gggTCGGTACCCCCAAGAAAAACCCCCCCTTTTTTTGGGATATCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAg < 1:100665/90‑1 (MQ=255) ggTGGGTACCACCCAGAAAAACCCCCCCCTGTTTTGGGATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGt < 2:82320/90‑1 (MQ=255) cGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGc > 1:324723/1‑90 (MQ=255) atcaccCAGAAAAACCCCCCCCTTTTTGGGGATACCCGAATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCaa < 2:501121/90‑1 (MQ=255) tcCAGAAAAACCCCCTCCTTTTTTGGGATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACgg < 2:403167/90‑1 (MQ=255) ggAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCa < 1:458951/42‑1 (MQ=255) ggAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCa > 2:458951/1‑42 (MQ=255) ccccctTTTTGGGGATACCCGAATGACCCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCa < 2:512485/87‑1 (MQ=255) tcctGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGc > 1:373820/1‑90 (MQ=255) gggATACCCGCATGACCCGGCCCCCCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc < 1:103007/90‑1 (MQ=255) tatCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTg > 2:527662/1‑90 (MQ=255) ccacATGACCCGCCCCCCCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcg < 1:16512/87‑1 (MQ=255) | ACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCG > NZ_CP009273/294577‑294739 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |