Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,571,380 | 0 | A | C | 55.6% | 13.4 / 19.4 | 18 | V166G (GTA→GGA) | yddB | TonB‑dependent receptor plug domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base A (7/1); new base C (0/10); total (7/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.51e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.80e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCC > NZ_CP009273/1571292‑1571468 | cGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAc > 2:622942/1‑90 (MQ=255) tGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCcccc < 1:295021/90‑1 (MQ=255) tGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCcccc < 2:611870/90‑1 (MQ=255) ttAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGc > 2:321343/1‑90 (MQ=255) gTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGaa > 2:387894/1‑90 (MQ=255) aCGCGTAGTGCGATACCCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGAcc < 1:283513/90‑1 (MQ=255) aCGCGTAGTGCGATACCCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGAcc < 2:60883/90‑1 (MQ=255) gcgTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCaa > 2:295791/1‑90 (MQ=255) gATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCg > 1:444192/1‑90 (MQ=255) cccAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGcaca < 1:295791/90‑1 (MQ=255) caGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACAt < 2:287512/90‑1 (MQ=255) ttGATTTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTg < 1:101196/90‑1 (MQ=255) ttGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTg < 2:402740/90‑1 (MQ=255) ttGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTg < 1:55037/90‑1 (MQ=255) tttGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCa < 2:382220/90‑1 (MQ=255) tCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataa < 1:92120/90‑1 (MQ=255) aaTTACCCCGCCATTGAAGCGACCAAATTCAACATGCGCAAAACTGTCATAAACCATCACAATGGCCACGAAACTAACCTCAAGacaatc > 1:380657/1‑89 (MQ=255) tACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATAccc > 1:272731/1‑90 (MQ=255) | CGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCC > NZ_CP009273/1571292‑1571468 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |