| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 794,976 | 0 | T | C | 37.0% | 28.6 / 10.9 | 27 | intergenic (+54/‑113) | acrZ/modA | multidrug efflux pump accessory protein AcrZ/molybdate ABC transporter periplasmic binding protein |
| Reads supporting (aligned to +/- strand): ref base T (12/5); new base C (0/10); total (12/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.64e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGT > NC_000913/794916‑795041 | tCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTTCTGCCg > 1:180687/1‑68 (MQ=255)tCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCg > 1:1361063/1‑68 (MQ=255)tCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCg > 1:2556276/1‑68 (MQ=255) tGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGatat > 1:1185162/1‑68 (MQ=255) gagtGCCCGTTTAGTCGGGCTTTTTTTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTa < 1:2276366/65‑1 (MQ=255) gAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCATCTTCCTGCTTTTCCTGCCGATATTTTTTCTTa < 1:1382788/68‑1 (MQ=255) cccccGCTTAGTCGGGCATTTTTTTTTTCCCACCTTCCTGTTTTCCCTGCCGATATTTTTTCTTATCt < 1:3135158/66‑1 (MQ=255) gCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTa > 1:3386738/1‑68 (MQ=255) gCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTa > 1:1860583/1‑68 (MQ=255) gCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTa > 1:2303801/1‑68 (MQ=255) gCTTAGTGGGGCATTTTTTTTTTCTCACCTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTc < 1:841469/68‑1 (MQ=255) gCTTAGTCGGGCTTTTTTTTTTTCTCAACTTCCTGTTTTCCCTGCCGATATTGTTTCTTATCTACCTc < 1:3122480/68‑1 (MQ=255) gCTTAGTCGGGCTTTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTc < 1:1130579/68‑1 (MQ=255) gCTTAGTCGGGCTTTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTc < 1:3991842/68‑1 (MQ=255) gCTTAGTCGGGCTTTTTCTTTTTCTCAACTTCCTGCTTTCCCTGCCGATATTTTTTCGTATCTACCTc < 1:3826796/68‑1 (MQ=255) tttAGTCGGGCTTTTTCTTTTTCCAATCTACCTTTTTTTCCTGCCGATATTTTTTCTTATCTACCTca < 1:658491/67‑1 (MQ=255) gTCGGTCTTTTTCTTTTTCTCATCTTCCTCCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACaaa < 1:1152363/68‑1 (MQ=255) gTCGGGCATTTTTTTTTTCCCAACTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACaaa < 1:1920752/68‑1 (MQ=255) gTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACaaa > 1:1116362/1‑68 (MQ=255) gTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACaaa > 1:3810133/1‑68 (MQ=255) cGGGCTTTTTCTTTTTCTCACCTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg < 1:869478/68‑1 (MQ=255) cGGGCATTTTCTTTTTCTCAACTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg < 1:2163852/68‑1 (MQ=255) cGGGCATTTTCTTTTTCTCAACTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg < 1:370699/68‑1 (MQ=255) cGGGCATTTTCTTTTTCTCAACTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg < 1:1479387/68‑1 (MQ=255) cGGGCATTTTCTTTTTCTCAACTTCCTGCTTTCCCTGCCGATATTTTTTCTTACCTACCTCACAAAgg < 1:488928/68‑1 (MQ=255) gCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTa > 1:1663471/1‑68 (MQ=255) ctttttttttttctCATCTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAg < 1:1815462/66‑1 (MQ=255) ctttttctttttctCAACTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAg < 1:391681/66‑1 (MQ=255) ctttttctttttctCAACTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAg < 1:2041485/66‑1 (MQ=255) cATTTTTTTTTTGTCAACTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAg < 1:2532039/68‑1 (MQ=255) cATTTTCTTTTTCTCAACTTCCTTTTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAg < 1:2618537/68‑1 (MQ=255) cATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAg < 1:3471739/68‑1 (MQ=255) aTTTTTTTTTTCTCATCTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGc < 1:237202/68‑1 (MQ=255) aTTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGc > 1:2421028/1‑68 (MQ=255) aTTTTCTTTTCCTCAACTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGc < 1:90934/68‑1 (MQ=255) ttttctttttctCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCa < 1:3833730/68‑1 (MQ=255) ttttctttttccCAACTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCa < 1:495257/68‑1 (MQ=255) ttgtctttttctCAACTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCa < 1:3437232/65‑1 (MQ=255) tttttttttcgCATCTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCaa < 1:663238/68‑1 (MQ=255) ttttctCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAAct < 1:2998355/68‑1 (MQ=255) tttactCACCTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAAct < 1:2978613/68‑1 (MQ=255) tttctCAACTTCCTGCTTTTCCTGCCGACATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAActg > 1:2599484/1‑68 (MQ=255) tctcAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAActgct > 1:2021336/1‑68 (MQ=255) caccTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTggg < 1:3233424/65‑1 (MQ=255) gCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCg > 1:199709/1‑68 (MQ=255) gCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCg > 1:905288/1‑68 (MQ=255) tttCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGt < 1:2645527/68‑1 (MQ=255) | TCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGAGT > NC_000913/794916‑795041 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |