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.5% | 35.6 / 59.9 | 44 | V166G (GTA→GGA) | yddB | TonB‑dependent receptor plug domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base A (16/4); new base C (0/24); total (16/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.16e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATA > NZ_CP009273/1571292‑1571465 | cGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAc > 2:199389/1‑90 (MQ=255) ttttAACCAGTCCGAACGCGTAGTGCGATACCCCATTTTCCCCTTGCTATCATCAGCGTTGAAGCTTTTGATCTTTGCATCAATTCcccc < 1:666212/88‑1 (MQ=255) tGTTACCCAGCCCGAACGCGTAGTGCGATACCCCAATTCCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATTTTTGCATCAATTAcccc < 1:535661/90‑1 (MQ=255) tGTTAACCAGTCCGACCGCGTGGTGCGATACCCCATTTCCCCCTTGCTATCATCAGCGTTGAAGCTTTTGATCTTTGCATCAATTCcccc < 1:746545/90‑1 (MQ=255) tGTTAACCAGTCCGAACGCGTAGTGCGATACCCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCcccc < 1:59950/90‑1 (MQ=255) tGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAcccc < 2:925933/90‑1 (MQ=255) ttAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGc < 2:371930/90‑1 (MQ=255) gTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGaa > 2:953499/1‑90 (MQ=255) gTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGaa > 1:321265/1‑90 (MQ=255) gTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGaa > 2:502522/1‑90 (MQ=255) gTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGaa > 2:674004/1‑90 (MQ=255) gTCCGAACGCGTAGTGCGATAACCCAATTTCACCTCGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGaa > 1:479806/1‑90 (MQ=255) ccGAACGCGTAGTGCGATAACCCAATTCCCCCTTGCTACCATCAGCGTTGAAGCTTTTGATCTTTGCATAAATTCCCCCGCCATTGAAGc < 1:775216/90‑1 (MQ=255) cgcgTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCa > 2:23401/1‑90 (MQ=255) gTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAAt > 2:780751/1‑90 (MQ=255) gTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAAt > 1:787020/1‑90 (MQ=255) gTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAAt > 1:260678/1‑90 (MQ=255) gTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAAt > 2:819689/1‑90 (MQ=255) aGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGAGACCAAATTc > 2:795786/1‑90 (MQ=255) cGATACCCCATTTCCCCCTTCCTATCATCAGCGTTAAAGCTTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAAcc < 1:795786/90‑1 (MQ=255) cGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAAcc < 2:225213/90‑1 (MQ=255) aTAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCgg > 2:225713/1‑90 (MQ=255) aTTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAAc < 2:200293/90‑1 (MQ=255) ttGCTATCATCAGCGTGGAAGCGTTTGATCTTTGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATaa < 2:692590/90‑1 (MQ=255) gCTATCATCAGCGTTGAAGCTTTTGTTCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAg < 1:963519/90‑1 (MQ=255) caGCGTTGAAGCGTTTGATTTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACAt < 2:744320/90‑1 (MQ=255) caGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACAt < 2:378998/90‑1 (MQ=255) gAAGCTTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAg < 2:533497/90‑1 (MQ=255) aGCTTTTTATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTa < 1:743277/90‑1 (MQ=255) gTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGc > 2:439676/1‑90 (MQ=255) ttGATTTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTg < 1:199389/90‑1 (MQ=255) ttGATCTTTGCCTAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTg < 2:787020/90‑1 (MQ=255) ttGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTg < 2:27439/90‑1 (MQ=255) ttGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTg < 2:69959/90‑1 (MQ=255) tttGCATCATTTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCa < 1:101913/90‑1 (MQ=255) tttGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCa < 2:177625/90‑1 (MQ=255) tttGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCa < 2:321265/90‑1 (MQ=255) tttGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCa < 2:241426/90‑1 (MQ=255) tttACATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCa < 1:421414/90‑1 (MQ=255) ttGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc < 2:14650/88‑1 (MQ=255) ttGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc > 1:14650/1‑88 (MQ=255) cAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataat < 2:437652/90‑1 (MQ=255) aaTTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataata < 1:225713/90‑1 (MQ=255) aaTTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataata > 2:357827/1‑90 (MQ=255) | CGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATA > NZ_CP009273/1571292‑1571465 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |