Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,571,380 | 0 | A | C | 54.8% | 28.9 / 38.9 | 31 | V166G (GTA→GGA) | yddB | TonB‑dependent receptor plug domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base A (10/4); new base C (0/17); total (10/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.26e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAG > NZ_CP009273/1571294‑1571459 | aTGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAccc > 1:16077/1‑90 (MQ=255) tGTTAACCAGTCCGAACGCGTAGTGCGATACCCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCcccc < 2:958244/90‑1 (MQ=255) tGTTAACCAGTCCGAACGCGTAGTGCGATACCCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCcccc < 2:611492/90‑1 (MQ=255) tGTTAACCAGTCCGAACGCGTAGTGCGATACCCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCcccc < 2:513290/90‑1 (MQ=255) tGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTCCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCcccc < 2:1175617/90‑1 (MQ=255) ttAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATTTTTGCATCAATTACCCCGc < 2:710993/90‑1 (MQ=255) aaCCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCa > 1:750766/1‑90 (MQ=255) aGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGa > 2:61780/1‑90 (MQ=255) gTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGaa > 1:127593/1‑90 (MQ=255) ccGAACGCGTAGTGCGATACCCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGc < 1:410657/90‑1 (MQ=255) cgTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCaaa > 2:225733/1‑90 (MQ=255) cGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAAcc < 1:1005021/90‑1 (MQ=255) ccAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACaa > 1:211246/1‑90 (MQ=255) cAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTCCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaa < 2:404888/90‑1 (MQ=255) aTTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAAc < 1:180197/90‑1 (MQ=255) aTTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAAc < 1:61780/90‑1 (MQ=255) aTTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAAc > 2:252180/1‑90 (MQ=255) aTTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAAc < 1:252180/90‑1 (MQ=255) ttccTATCATCAGCGTTAAAGCTTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATaa < 2:16077/87‑1 (MQ=255) ttGCTATCACCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATaa < 1:645490/90‑1 (MQ=255) atcatcAGCGTTGAAGCTTTTGATTTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGt < 1:503577/90‑1 (MQ=255) caGCGTTGAAGCGTTTGATTTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACAt < 1:610431/90‑1 (MQ=255) gCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTg > 1:903609/1‑90 (MQ=255) gCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTg > 2:413107/1‑90 (MQ=255) gAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGTAGCGACCAAGGTCAACCGCCACAAAATTTTCATAAAACGTCACATTGTCCAg > 2:799740/1‑90 (MQ=255) aaGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGt < 1:413107/90‑1 (MQ=255) ttGATCTTTCCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTg < 2:903609/90‑1 (MQ=255) aTCTTTACAACAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACa < 2:750766/90‑1 (MQ=255) tttGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCa < 1:143102/90‑1 (MQ=255) tttGCATAATTTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCa < 2:403653/90‑1 (MQ=255) tGCATCAATCCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAg < 2:612518/90‑1 (MQ=255) | ATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAG > NZ_CP009273/1571294‑1571459 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |