Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 443,152 | 0 | T | A | 54.5% | 15.0 / 32.0 | 22 | K4* (AAG→TAG) | cyoE | heme o synthase |
Reads supporting (aligned to +/- strand): ref base T (9/1); new base A (0/12); total (9/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.01e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACAACCAGGATAGCGA > NZ_CP009273/443065‑443240 | tgcCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTaa < 1:25836/90‑1 (MQ=255) cccTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAc > 1:157513/1‑90 (MQ=255) cccTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAc > 1:909996/1‑90 (MQ=255) gAGGCCAGCAGGACCCCCCCAACCCCCGAGATCGGTTTGCAAAAGATGATGCCTGTTTTCGTTTCTTGCGGGTATTGCTAAAACATCATa < 1:538454/90‑1 (MQ=255) agcagGACCCCCCCACCCCCCGAGATCGGGTTGCAAAAGATGATGCCTGTTTTCTTTTCTTGCAGTTTTTGTTAAAACATCAAACccgcc < 1:565952/90‑1 (MQ=255) aaccccccAACCACCGATTTCGGGTAACCAAAGATGGTGCTTGTTTTTCTTTCTTGTGGTTATTTTTAAAACATCATAACCGCCGCTCtt < 1:601969/87‑1 (MQ=255) cccccAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGt > 2:476793/1‑90 (MQ=255) cccAACCCCCGAGATCGGGTTGCAAAAGATGATGCCTGTTTTCTTTACTTGCGGGTTTTGCTAAAACATCATAACCGCCGCTCTTAGTGc < 1:177529/90‑1 (MQ=255) cAATCCCCGAGATCGGGTTGCAAAAGATGATCCCTGGTTTCGTTCTTTGCGGGTATTGCTAAAACATCATAACCGCCGCTCTTAGTGcat < 2:81387/90‑1 (MQ=255) ccGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGcatcatcat > 2:1036694/1‑90 (MQ=255) gagaTCGGGTTCCAAAAGTTGATCCTTGTTTTCTTTCCTTGCGGTTTTTTTTAAAAAATCAAACCCGCCGCTCTTAGTGCATCATCAtgt < 1:183143/90‑1 (MQ=255) gaaaTCAGGTTGCAAAAGATGATGCCTGTTTTCTTTTCTTGCAGTTTTTTTTAAAACAACATACCCGCCGCTCTTAGTGCATCATCAtgt < 1:199888/87‑1 (MQ=255) ggTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGtt > 2:53236/1‑90 (MQ=255) tGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAg > 1:539821/1‑90 (MQ=255) tGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAg > 2:452844/1‑90 (MQ=255) aaaGATGATCCCTGTTTTCTTTTCTTGCAGTTTTTGCTAAAACATAAAACCCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTc < 1:271474/90‑1 (MQ=255) gttgCTTGTTTTCTTTACTTGGGGGTTTTGCTAAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATa < 2:909996/88‑1 (MQ=255) gatgCCTGTTTTCTTTACTTGCAGGTTTTGCTAAAACATAATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATa < 2:157513/90‑1 (MQ=255) gatgCCTGTTTTCGTTCCTTGCAGTTTTTTCTAAAACATAAAACCCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATa < 1:84021/90‑1 (MQ=255) gatcCCTGTTTTCTTTCCTTGCAGTTTTTCTTAAAACATCAAAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATa < 2:89571/90‑1 (MQ=255) atgCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATaa > 2:188729/1‑90 (MQ=255) ttGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACAACCAGGATa > 2:793225/1‑90 (MQ=255) ttAATCATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACAACCAGGATAGCga < 2:75801/90‑1 (MQ=255) | TGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACAACCAGGATAGCGA > NZ_CP009273/443065‑443240 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |