Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,015,516 | 0 | G | T | 25.6% | 199.7 / 31.0 | 82 | intergenic (+57/‑199) | pqiC/rmf | intermembrane transport lipoprotein PqiC/ribosome modulation factor |
Reads supporting (aligned to +/- strand): ref base G (24/37); new base T (21/0); total (45/37) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.46e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACC > NC_000913/1015468‑1015565 | aTTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGt < 1:7398721/50‑1 (MQ=255) aTTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGt > 1:1783858/1‑50 (MQ=255) tttGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGtc < 1:3598448/50‑1 (MQ=255) tGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTtctc > 1:3760317/1‑50 (MQ=255) tGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTtctc > 1:3321415/1‑50 (MQ=255) tGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGtct < 1:2536504/49‑1 (MQ=255) gTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTtctc > 1:1891695/1‑49 (MQ=255) gTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGtctc > 1:3455721/1‑49 (MQ=255) tAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTtctcc > 1:3009924/1‑48 (MQ=255) tAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTttct > 1:4200461/1‑49 (MQ=38) tAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTccct > 1:3197581/1‑49 (MQ=255) aaaTATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTTTCTaa > 1:4280231/1‑50 (MQ=39) aaaTATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTaa > 1:3437978/1‑50 (MQ=255) aaaTATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTa > 1:4200567/1‑49 (MQ=255) aaTATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTaa < 1:7821523/49‑1 (MQ=255) atatAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACg < 1:6038050/50‑1 (MQ=255) tataACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGt > 1:6169093/1‑50 (MQ=255) tataACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGt < 1:2067753/50‑1 (MQ=255) tataACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGt < 1:6665543/50‑1 (MQ=255) tataACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGt > 1:5841006/1‑50 (MQ=255) tataACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGt < 1:8555644/50‑1 (MQ=255) tataACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGt < 1:5814960/50‑1 (MQ=255) tataACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGt < 1:3557542/50‑1 (MQ=255) tataACCGTCTCCGGTATGGTGCCTGAGGCGGTTTTTTTGTCTCTAACGt < 1:797034/50‑1 (MQ=255) taACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTTACGTGc > 1:5878052/1‑50 (MQ=39) taACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGc > 1:658886/1‑50 (MQ=255) taACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTg < 1:4803370/49‑1 (MQ=255) taACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTg < 1:2254233/49‑1 (MQ=255) taACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGc < 1:7890087/50‑1 (MQ=255) taACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGc < 1:4236062/50‑1 (MQ=255) aaCCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCg > 1:3626473/1‑50 (MQ=255) aaCCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCg > 1:4919720/1‑50 (MQ=255) aaCCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTTTCTAACGTGCg > 1:8230408/1‑50 (MQ=255) aaCCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTTTCCAACGGGCg > 1:4410045/1‑50 (MQ=38) aaCCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGGGCg > 1:1696061/1‑50 (MQ=255) aCCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTTTCTTACGTGCgg > 1:2089395/1‑50 (MQ=39) aCCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCgg > 1:6152172/1‑50 (MQ=255) aCCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCgg < 1:3809666/50‑1 (MQ=255) ccGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCCAACGTGCgg > 1:785399/1‑49 (MQ=255) cGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACCTGCggg > 1:8535986/1‑48 (MQ=38) agTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTTACGTGCggg > 1:177816/2‑48 (MQ=39) gTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCggga > 1:6780555/1‑47 (MQ=255) gTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGGCTCTAACGTGCGGaaa < 1:5467209/50‑1 (MQ=255) tctcCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGaaaa < 1:5009262/50‑1 (MQ=255) tctcCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGaaaa < 1:5435160/50‑1 (MQ=255) ctcCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCGGaaaaa > 1:8326823/1‑50 (MQ=255) tcCGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCGGGAAAAt > 1:700869/1‑50 (MQ=39) ccGGTATGTTGCCTGAGGCGGTTTTTTTTTCTTCAACGGCCGGaaaaaca > 1:3362799/1‑48 (MQ=25) ccGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCGGGAAAAtt > 1:1383315/1‑50 (MQ=39) ccGGTATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCGGAAAAAtt > 1:3055019/1‑50 (MQ=255) cGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGaaaaa < 1:6120157/47‑1 (MQ=255) cGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAAttt < 1:7858694/50‑1 (MQ=255) cGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAAttt < 1:1386945/50‑1 (MQ=255) cGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAAttt < 1:3723232/50‑1 (MQ=255) tATGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCGGGAAAAtttttt > 1:3930387/1‑50 (MQ=38) aTGTTGCCTGAGGCGGTTTTTTTTTCTCTAACGTGCGGGAAAAATTGttt > 1:4569082/1‑49 (MQ=38) aTGTTGCCTGAGGCGGTTTTTTTGTCTCTAAAGTGCGGAAAAATTTGTTc > 1:2030597/1‑50 (MQ=255) tGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTcc < 1:238345/50‑1 (MQ=255) tGTTGCCTGAGGCGGTTTCTTTGTCTCTAACGTGCGGAAAAATTTGTTcc < 1:5637094/50‑1 (MQ=255) ttGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCctc < 1:3116241/50‑1 (MQ=255) gCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCctct < 1:2196382/49‑1 (MQ=255) gCCTGAGGCGGTTTTTTTGTCTCTAACCTGCGGAAAAATTTGTTCCTCtt > 1:4460431/1‑50 (MQ=255) ccTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCtt < 1:1027788/49‑1 (MQ=255) aGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTcaca < 1:6520365/49‑1 (MQ=255) aGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTcaca < 1:7885433/49‑1 (MQ=255) aGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAt < 1:8515594/50‑1 (MQ=255) aGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAt < 1:6594330/50‑1 (MQ=255) ggCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAtt > 1:5559410/1‑50 (MQ=255) ggCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAtt > 1:3871258/1‑50 (MQ=255) gCGGTTTTTTTGTCTCTAACGGGCGGAAAAATTTGTTCCTCTTCACAttt < 1:7296376/50‑1 (MQ=255) cGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAtttt > 1:6050695/1‑50 (MQ=255) cGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAtttt < 1:3365621/50‑1 (MQ=255) ggTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAttttt > 1:1043625/1‑50 (MQ=255) ggTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAttttt > 1:1481393/1‑50 (MQ=255) gTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACAtttttt > 1:2456415/1‑50 (MQ=255) tttttttGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTg > 1:3431863/1‑50 (MQ=255) ttttttGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTg > 1:5893594/1‑49 (MQ=255) tttttGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGt < 1:4619211/49‑1 (MQ=255) tttttGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGt < 1:6059817/49‑1 (MQ=255) tttGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTAc > 1:4392542/1‑49 (MQ=255) ttGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACaa > 1:5914334/1‑50 (MQ=255) ttGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACa < 1:2093078/49‑1 (MQ=255) gTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAAcc > 1:7730911/1‑50 (MQ=255) gTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAAc < 1:6615980/49‑1 (MQ=255) | ATTTGTAAATATAACCGTCTCCGGTATGTTGCCTGAGGCGGTTTTTTTGTCTCTAACGTGCGGAAAAATTTGTTCCTCTTCACATTTTTTGTACAACC > NC_000913/1015468‑1015565 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |