Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,571,380 | 0 | A | C | 57.9% | 10.8 / 21.5 | 19 | V166G (GTA→GGA) | yddB | TonB‑dependent receptor plug domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base A (8/0); new base C (0/11); total (8/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.32e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATAC > NZ_CP009273/1571292‑1571466 | cGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAc > 1:367351/1‑90 (MQ=255) tGTTACCCAGTCCGAACGCGTAGTGCGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCcccc < 2:43746/90‑1 (MQ=255) tGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCcccc < 2:541736/90‑1 (MQ=255) gTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGaa > 1:574442/1‑90 (MQ=255) gTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGaa > 1:64879/1‑90 (MQ=255) aaCGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGAc > 2:199821/1‑90 (MQ=255) gCGATAACCCAATTTCCCCTTATTATCCCCAGCGTTGACGCTTTTGTTTTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAAc < 2:502852/90‑1 (MQ=255) cGATAACCCAATTCCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATTTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAAcc < 2:440572/90‑1 (MQ=255) ttGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATaa < 2:281762/90‑1 (MQ=255) atcatcAGCGTGGAAGCTTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGt < 2:571775/90‑1 (MQ=255) atcaccAGCGTTAAAGCTTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGt < 2:23285/90‑1 (MQ=255) caGCTTGGAGGCTTTTGATCTTTCCACAATTTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACAt < 1:329268/90‑1 (MQ=255) caGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACAt > 1:239211/1‑90 (MQ=255) caGCGTTAAAGCTTTTGATCTTTCCAAAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACAt < 1:38781/90‑1 (MQ=255) ttttGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc < 2:64879/89‑1 (MQ=255) tttGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCa > 2:99999/1‑90 (MQ=255) cAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataat < 1:69779/90‑1 (MQ=255) aaTTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataata > 1:435123/1‑90 (MQ=255) aTTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATAc > 2:355829/1‑90 (MQ=255) | CGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATAC > NZ_CP009273/1571292‑1571466 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |