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.0% | 37.3 / 51.6 | 37 | V166G (GTA→GGA) | yddB | TonB‑dependent receptor plug domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base A (15/2); new base C (0/20); total (15/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.45e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCC > NZ_CP009273/1571291‑1571468 | gCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTa > 2:108875/1‑90 (MQ=255) aTGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTAccc > 1:156022/1‑90 (MQ=255) tGTTAACCAGTCCGAACGCGTAGTGCGATACCCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATTTTTGCATAAATTCcccc < 1:896850/90‑1 (MQ=255) tGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCcccc < 1:182586/90‑1 (MQ=255) tGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCTTTTGATCTTTGCATCAATTCcccc < 1:661557/90‑1 (MQ=255) tGTTAACCAGTCAGAACGCCTGGTGCGATAACCCAATTTCTCCTTGCTTTCACCAGCGTGGAGGCGTTTGATCTTTGCATCAATTCcccc < 1:1283017/90‑1 (MQ=255) gTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGaa > 1:165013/1‑90 (MQ=255) gAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGa > 2:546808/1‑90 (MQ=255) aCGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGAcc < 2:1045715/90‑1 (MQ=255) gTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAAt > 2:927705/1‑90 (MQ=255) gTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAAt > 2:580100/1‑90 (MQ=255) cGATACCCCAATTCCCCCTTGCTATCATCAGCGTTGAAGCTTTTGATCTTTGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAAcc < 1:108875/90‑1 (MQ=255) cGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCTTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAAcc < 2:1366386/90‑1 (MQ=255) cGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAAcc < 2:1264312/90‑1 (MQ=255) gATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCg > 2:468672/1‑90 (MQ=255) aaTTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACaaaa > 2:19921/1‑90 (MQ=255) aTTTCCCCTTCCTACCACCGGCGTTGAAGCGTTTGTTCTTTCCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAAc < 1:927705/90‑1 (MQ=255) aTTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAAc < 2:1142924/90‑1 (MQ=255) ttGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATaa < 2:958589/90‑1 (MQ=255) catcaGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTca > 2:649573/1‑90 (MQ=255) caGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACAt < 2:40649/90‑1 (MQ=255) caGCGTTGAAGCGTTTGATCTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACAt < 2:477174/90‑1 (MQ=255) caGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACAt > 2:539607/1‑90 (MQ=255) gTTGAAGCTTTTGATCTTTGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTc < 2:792579/90‑1 (MQ=255) gAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAAcc < 2:289650/55‑1 (MQ=255) gAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAAcc > 1:289650/1‑55 (MQ=255) gAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAg > 1:792579/1‑90 (MQ=255) gAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAg > 1:137423/1‑90 (MQ=255) ttttGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc < 1:418160/89‑1 (MQ=255) ttttGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc < 1:548210/89‑1 (MQ=255) cTTTGCATCAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc < 2:81141/90‑1 (MQ=255) cTTTGCATAAATTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATc < 2:882307/90‑1 (MQ=255) tttGCAACATTTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCa < 1:600883/90‑1 (MQ=255) ttGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCaa > 2:851288/1‑90 (MQ=255) aaTTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGataata < 2:1346335/90‑1 (MQ=255) aTTCCCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACAt < 2:643568/61‑1 (MQ=255) aTTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACAt > 1:643568/1‑61 (MQ=255) tACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATAccc > 1:143150/1‑90 (MQ=255) | GCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGCACAAAACTGTCATAAAGCGTCACATTGTCCAGTAAGCTGACATCAAGATAATACCC > NZ_CP009273/1571291‑1571468 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |