Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,197,867 | 0 | T | G | 29.5% | 89.0 / 23.8 | 44 | A47A (GCA→GCC) | cirA | ferric iron‑catecholate outer membrane transporter |
Reads supporting (aligned to +/- strand): ref base T (5/26); new base G (13/0); total (18/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.65e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTC > REL606/2197817‑2197917 | aCCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGt < 1:3635714/51‑1 (MQ=255) cGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGc < 1:4439008/51‑1 (MQ=255) gTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCAt < 1:775578/51‑1 (MQ=255) gTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCAt < 1:5870178/51‑1 (MQ=255) gTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCAt < 1:5804802/51‑1 (MQ=255) cAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGAtt < 1:2913203/51‑1 (MQ=255) cAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCttt > 1:6315367/1‑45 (MQ=255) aGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGAttt < 1:1141196/51‑1 (MQ=255) aGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCttt > 1:6056126/1‑44 (MQ=255) ggTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGAtttt < 1:746209/51‑1 (MQ=255) ggTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGAtttt < 1:3775704/51‑1 (MQ=255) gTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTg < 1:1346146/51‑1 (MQ=255) tCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGCTCTTTTAGAtttttt > 1:4692208/1‑49 (MQ=25) tCTTCCGGGGTAATGACGCTGATACTGGCAGGTGCATCtt < 1:1757143/40‑1 (MQ=255) ttCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTc < 1:2003397/51‑1 (MQ=255) tCCTGGGTAATGACGCTGATACTGGCAGGGGGCTCttt > 1:596475/1‑38 (MQ=25) cTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTTAGATtttttt > 1:3110041/1‑47 (MQ=255) tGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACg < 1:5395361/51‑1 (MQ=255) gggTAATGACGCTGATACTGGCAGGGGGATCTTTAAAAtttt > 1:617191/1‑42 (MQ=37) gTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGaa < 1:5281870/51‑1 (MQ=255) gTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGaa < 1:4660957/51‑1 (MQ=255) gTAATGACGCTGATACTGGCAGGGGGATCTTTTAGATtttttt > 1:3323858/1‑43 (MQ=255) tAATGACGCTGATACTGGCAGGGGcttctt > 1:5917588/1‑30 (MQ=25) tAATGACGCTGATACTGGCAGGGGGCTCTTTTAGATtttttt > 1:979320/1‑42 (MQ=255) tAATGACGCTGATACTGGCAGGGGCATCTTTTAGATtttttt > 1:3222997/1‑42 (MQ=37) tAATGACGCTGATACTGGCAGGGGCAGCTTTAAGAtttt > 1:6721506/1‑39 (MQ=37) aaTGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGa < 1:699076/51‑1 (MQ=255) aTGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGAt < 1:5385028/51‑1 (MQ=255) aTGACGCTGATACTGGCAGGGGCATCTTTAAGATTTTGTTTCCCgg > 1:4497592/1‑46 (MQ=38) tGACGCTGATACTGGCAGGGGCATCTTTTAGATtttttt > 1:3765749/1‑39 (MQ=25) cTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTg < 1:2817242/51‑1 (MQ=255) tACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACaa < 1:2082449/51‑1 (MQ=255) cTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACaa > 1:3072945/1‑49 (MQ=255) cTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAAcc < 1:3448301/51‑1 (MQ=255) cTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAAcc > 1:6334480/1‑51 (MQ=255) tGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCa < 1:5103929/51‑1 (MQ=255) ggCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCAt < 1:3214444/51‑1 (MQ=255) gCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATc < 1:499490/51‑1 (MQ=255) cAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCg > 1:551215/1‑51 (MQ=255) cAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCg < 1:4028141/51‑1 (MQ=255) aGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGt > 1:7283663/1‑51 (MQ=255) aGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGt > 1:986489/1‑51 (MQ=255) ggTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGtt < 1:1982211/51‑1 (MQ=255) gTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGttt < 1:5751632/51‑1 (MQ=255) gTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGttt < 1:6539409/51‑1 (MQ=255) tGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTc > 1:6420963/1‑51 (MQ=255) | ACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTC > REL606/2197817‑2197917 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |