Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 794,976 | 0 | T | C | 32.4% | 48.6 / 12.2 | 34 | 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 (16/7); new base C (0/11); total (16/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.27e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGA > NC_000913/794916‑795039 | tCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTg > 1:590485/1‑65 (MQ=255) tCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTg > 1:566875/1‑65 (MQ=255) tGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGa > 1:1448417/1‑65 (MQ=255) tGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGa > 1:1677245/1‑65 (MQ=255) aTTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGata > 1:493301/1‑65 (MQ=255) aTTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGata > 1:697210/1‑65 (MQ=255) tCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATAtt > 1:2036618/1‑65 (MQ=255) aaTGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCt < 1:549943/65‑1 (MQ=255) aaTGCCCGCTTAGTCGGGCATTTTCTTTTTCAAATCTTCCTGCTTTTCCTGCCGATATTTTTTCt < 1:1597850/65‑1 (MQ=255) aaTCCCCGCTTAGTCGGGCTTTTTTTTTTTCTCACCTTCCTGCTTTCCCTGCCGATATTTTTTCt < 1:1601407/65‑1 (MQ=255) tGCCCGCTTTTTCGGACTTTTTCTTTTTCTCAACTTCTCGCTTTCCCTGCCGATATTTTTTCTTa < 1:1680492/65‑1 (MQ=255) gCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTAt > 1:1276510/1‑65 (MQ=255) ccccGTTTGGTCGGGTTTTTTCTTTTCCACAACTCCCTGTTTTCCCTGCCGATATTTTTTCTTAt < 1:1761106/64‑1 (MQ=255) gCTTAGTCGGGCTTTTTCTTTTTCTCAACTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTAc < 1:424367/65‑1 (MQ=255) gCTTAGTCGGGCATTTTTTTTTTCCCAACTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTAc < 1:79290/65‑1 (MQ=255) gCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTAc < 1:1867111/65‑1 (MQ=255) tttAGTCGGGCTTTTTCTTTTTCTCACCTTCCTGTTTTTCCTGCCGATATTTTTTCTTATCTAcc < 1:673379/64‑1 (MQ=255) ttAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCt > 1:823741/1‑65 (MQ=255) ttAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCt > 1:1536302/1‑65 (MQ=255) ttgTCGGGCATTTTTTTTTTTTCAAATTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTc < 1:699311/63‑1 (MQ=255) tcgTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTc < 1:800514/63‑1 (MQ=255) tAGTCGTGCTTTTTCTTTTTCTCACCTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTc < 1:3648927/65‑1 (MQ=255) tAGTCGGGCTTTTTTTTTTTCTCACCTTCCTCTTTTTCCTGCCGATATTTTTTCTTATCTACCTc < 1:1006738/65‑1 (MQ=255) tAGTCGGGCTTTTTCTTTTTCTCACCTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTc < 1:1405087/65‑1 (MQ=255) tAGTCGGGCTTTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTc < 1:2105005/65‑1 (MQ=255) tAGTCGGGCATTTTTTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTc < 1:479368/65‑1 (MQ=255) tAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTc < 1:2457434/65‑1 (MQ=255) tAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTc < 1:1483537/65‑1 (MQ=255) cGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACaa > 1:1851342/1‑65 (MQ=255) gggCTTTTTCTTTTTCTCAACTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACaaa < 1:1853152/65‑1 (MQ=255) tctttttttttttccAATCTTCCTGTTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg < 1:720074/62‑1 (MQ=255) tctttttctttttctCAACTTCCTTCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg < 1:1178109/62‑1 (MQ=255) gttttttttttttccCACCTGCCTGTTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg < 1:3196703/62‑1 (MQ=255) gtaTTTTTTTTTTCCAAACTTCCTGTTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg < 1:322599/63‑1 (MQ=255) gctttttttttttccCACCTCCCTTTTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg < 1:2262644/62‑1 (MQ=255) gctttttctttttctCATCTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg < 1:1199980/62‑1 (MQ=255) gctttttctttttctCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg < 1:3584382/62‑1 (MQ=255) gctttttctttttccCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg < 1:1647640/62‑1 (MQ=255) gctttttcttttcctCACCTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg < 1:1955711/62‑1 (MQ=255) gCATTTTCTTTTTCTCATCTTCCTATTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg < 1:2184008/65‑1 (MQ=255) gCATTTTCTTTTTCTCAACTCCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAgg < 1:3543814/65‑1 (MQ=255) cttttttttttccaCAACCTCCCTCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGt < 1:3345907/63‑1 (MQ=255) aTTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGtt > 1:2497613/1‑65 (MQ=255) aTTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGtt > 1:461288/1‑65 (MQ=255) aTTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGtt > 1:2135124/1‑65 (MQ=255) tttctttttctCACCTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAg < 1:1882814/65‑1 (MQ=255) ttctttttctCAACTGCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGc < 1:463734/65‑1 (MQ=255) tttttttctCAACTCCCTGGTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCa < 1:3077304/63‑1 (MQ=255) tctttttctCACCTTCCTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCa < 1:2378618/65‑1 (MQ=255) tctttttctCAACTTCCTTTTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCa < 1:2339592/65‑1 (MQ=255) tctttttccCATCTGCCTCCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCa < 1:3162726/65‑1 (MQ=255) tttttctCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAAt > 1:393303/1‑65 (MQ=255) ctcccTGTTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTggg < 1:1076645/62‑1 (MQ=255) cTGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTNGGaaaa < 1:2221315/65‑1 (MQ=255) tGCTTTCCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAAt < 1:3701280/65‑1 (MQ=255) gCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAAtt > 1:3308067/1‑65 (MQ=255) tttcCCTGCCGATATTTTTTCTTATCTGCCTCACACAGGTTAGCAATAACTGTTGGGAAAATTcc < 1:2324157/65‑1 (MQ=255) tttCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCg > 1:731190/1‑65 (MQ=255) ttCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGa > 1:3294681/1‑65 (MQ=255) | TCATTGATTCCCTGAATGCCCGCTTAGTCGGGCATTTTCTTTTTCTCAACTTCCTGCTTTTCCTGCCGATATTTTTTCTTATCTACCTCACAAAGGTTAGCAATAACTGCTGGGAAAATTCCGA > NC_000913/794916‑795039 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |