| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,744,897 | 0 | A | C | 59.3% | 24.6 / 35.3 | 27 | V111G (GTT→GGT) | yiaT | MipA/OmpV family protein |
| Reads supporting (aligned to +/- strand): ref base A (11/0); new base C (0/16); total (11/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.67e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-06 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTT‑CCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAG > NZ_CP009273/3744808‑3744986 | gCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTT‑CCCCAACGCTCGTGGTGATACCAGGCaa < 2:586738/90‑1 (MQ=255) ccccccccCAGCCGTTGCTTTTACCCAGAACGCCCGCAGCGGCAGGGGTTTTCCCGGTC‑CCCCACCGCTGGTGGTGATACCAGGCAcccc < 2:340321/90‑1 (MQ=255) ccccACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTT‑CCCCAACGCTCGTGGTGATACCAGGCAAcccc > 2:270568/1‑90 (MQ=255) cccccccaGCCGTTGCTGTTACCCAAACCGCCCGCAGCTGCAGGGGCTTTCACGGTC‑CCCCAACGCTCGTGGTGATCCCAGGCAcccccc < 2:715030/90‑1 (MQ=255) acccaGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTT‑CCCCAACGCTCGTGGTGATACCAGGCAACCCCCGcc > 1:418154/1‑90 (MQ=255) cccaCCCGTGGCTGTTACCAAGACCGCCCCCCGCTGCAGGGGCTTTCCCGGTC‑CCCCACCGCTCGTGGTGATACCGGGCACCCCCCGCCa < 2:79649/90‑1 (MQ=255) ccaGCCTTTGCTGTTACCCAGAACGCCCGCAGCTGCAGGGGTTTTCACGGTT‑CCCCAACGCTCGTGGTGATACCAGGCACCCCCCGCCAt < 1:336294/90‑1 (MQ=255) caGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTT‑CCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATa > 1:265702/1‑90 (MQ=255) gCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTT‑CCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGc > 1:637044/1‑90 (MQ=255) ccctttGCTTTTACCCAGACCGCCCGCAGCGGCAGGGGCTTCCACGGTC‑CCCCACCGCTCGGGGTGATACCAGCCACCCCCCGCCATAGc < 2:580426/86‑1 (MQ=255) cccGTGGCTTTTACCCAGACCCCCCCCACCGGCAGAGGCTTCCCCGGTC‑CCCCAACCCTCGTGGTGATCCCAGCCACCCCCCGCCATAGc < 2:534609/89‑1 (MQ=255) gttgttgttACCCAGAACGCCCGCAGCTGCAGAGGCTTTCCCGGTC‑CCCCAACGCTCGTGGTGATCCCAGGCACCCCCCGCCATAGCCGt < 1:447969/90‑1 (MQ=255) gTTGCTTTTATCCAAAACGCCCGCAGCTGCAGGGGCTTTCGCGTTT‑CCCAACCGCTCGGGGTGATCCCAGGCACCCCCCGCCATAGCCGt < 1:322253/90‑1 (MQ=255) tGTTATCCAGAACGCCCGCAGCTGCAGAGGCTTCCCCGGTC‑CCCCACCGCGCGGGGTGATACCAGGCACCCCCCGCCATAGCCGTACTAt < 2:418154/90‑1 (MQ=255) accaaGACCGCCCGCCGCTGCAGGGGTTTTCGGGGTT‑CCCCACGGCGCGGGGTGATCCCGGCCACCCCCCGCCATAGCCGTACTATCACg < 1:542575/86‑1 (MQ=255) accaaGAACGCCCGCAGCGGCGGGGGCTTTCCCGGTC‑CCCCACCGCTGGTGGTGATCCCAGCCACCCCCCGCCATAGCCGTACTATCACg < 2:464463/86‑1 (MQ=255) gAACGTCCGCAGCTGCAGAGGCTTTCACGGTT‑CCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTAt > 2:466079/1‑90 (MQ=255) cccccGCAGCTGCAG‑GGTTTTCCCCGTTCCCCCACCGCTGGTGGTGATCCCGGCCACCCCCCGCCATAGCCGTACTATCACGCTTATCaa < 2:406470/87‑1 (MQ=255) cGTCCGCAGCTGCAGGGGCTTTCCCGGTC‑CCCCACCGCCGGTGGTAATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCaa < 2:4418/90‑1 (MQ=255) cccGCAGCTGCAGAGGTTTTCCCGGTC‑CCCCACCGCTCGTGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGc < 1:376657/89‑1 (MQ=255) cccGCAGCTGCAGAGGCTTTCACGGTC‑CCCCAACGCTCGGGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGc < 1:221172/89‑1 (MQ=255) ggCAGAGGTTTTCCCGGTC‑CCCCAACGCTCGTGGTGATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCAt < 1:527529/89‑1 (MQ=255) cTTTCACGGTT‑CCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAAt > 2:533411/1‑90 (MQ=255) cTTTCACGGTT‑CCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAAt > 1:361543/1‑90 (MQ=255) tttCACGGTT‑CCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATc > 2:290179/1‑90 (MQ=255) gtcccccACCGCTGGTGGTGATCCCGGCCCCCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTa < 2:610602/87‑1 (MQ=255) accGCTGGGGGTAATCCCAGGCACCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAg < 2:219947/88‑1 (MQ=255) gCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGt > 2:527546/1‑90 (MQ=255) gATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCa > 2:301411/1‑90 (MQ=255) gATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCa > 1:108120/1‑90 (MQ=255) aCCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAg > 2:471535/1‑90 (MQ=255) | GCTCCCCCACCCAGCCGTTGCTGTTATCCAGAACGTCCGCAGCTGCAGAGGCTTTCACGGTT‑CCCCAACGCTCGTGGTGATACCAGGCAACCCCCGCCATAGCCGTACTATCACGCTTATCAAGCTGTTGCATGGCATAATCGTCATTATCGGTAGGGTCAAATTCCAGCGGCATCCAG > NZ_CP009273/3744808‑3744986 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |