Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,197,867 | 0 | T | G | 33.3% | 65.6 / 13.6 | 36 | A47A (GCA→GCC) | cirA | ferric iron‑catecholate outer membrane transporter |
Reads supporting (aligned to +/- strand): ref base T (2/22); new base G (12/0); total (14/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.27e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTTGAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAAC > REL606/2197792‑2197939 | tttGAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGt < 2:262318/76‑1 (MQ=255) tttGAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGt < 2:1764286/76‑1 (MQ=255) tGAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTgg > 2:372066/1‑75 (MQ=255) gAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCa > 2:122863/1‑76 (MQ=255) ggACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATc < 1:441648/76‑1 (MQ=255) aCATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGGATCtt > 1:1748267/1‑76 (MQ=255) aCATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCtt < 2:735566/76‑1 (MQ=255) aCATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCtt < 1:1707310/76‑1 (MQ=255) cATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCttt > 2:1295683/1‑76 (MQ=255) aTCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTa > 2:1618493/1‑76 (MQ=255) aTCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCtttt > 1:1746948/1‑75 (MQ=255) tCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTaa < 2:238251/76‑1 (MQ=255) ccTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAg < 1:589173/76‑1 (MQ=255) cTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGa < 1:212593/76‑1 (MQ=255) cTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGa < 1:1374182/76‑1 (MQ=255) cTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTTAGa > 1:613969/1‑76 (MQ=255) ttCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGAt < 2:600607/76‑1 (MQ=255) cAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTAAGAttt > 2:7922/1‑76 (MQ=255) ttCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTTAGATtttttt > 1:1424743/1‑76 (MQ=255) tCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTc < 1:233419/76‑1 (MQ=255) cTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTcn < 1:1595968/76‑2 (MQ=255) cTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTcc < 2:609716/76‑1 (MQ=255) tGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTTTTTCa > 1:1021444/1‑76 (MQ=255) tACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTTAGATTTTTTTTCACg > 2:1204883/1‑76 (MQ=255) tACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTTTTTCACg > 1:1728976/1‑76 (MQ=255) aCCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTTTTCCACgg > 2:19817/1‑76 (MQ=255) ccGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGa < 2:1147336/76‑1 (MQ=255) cGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCAGATACTGGCAGGGGGATCTTTTAGATTTTTTTTCACGGaa > 2:1259363/1‑76 (MQ=255) ggTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCACCTTTTAAATTTTTTTTCACGGAAg > 2:401949/1‑76 (MQ=255) ttttCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTTTTTCACCGCAGaa > 1:89482/1‑75 (MQ=255) tttCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTTAGATTTTGTTTCACGGAAGATg > 1:92074/1‑76 (MQ=255) tttCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTTTTTCACGGAAGaag > 1:1274420/1‑74 (MQ=255) gCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCGTCTTTTAGATTTTGTTTCACGGGAGATGCAgg > 2:393494/1‑75 (MQ=255) tGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGa < 2:833184/76‑1 (MQ=255) tGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTTTTTCACCGAAGAAGCAGTGa > 1:936044/1‑76 (MQ=255) tGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTGTTCCACGGAAGATGCAGGGa > 2:826306/1‑76 (MQ=255) gCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGGATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGaa > 1:1930002/1‑75 (MQ=255) gCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGGGCATCTTTTAGATTTTGTTTCACGGAAGAAGCAGTGaa > 1:692492/1‑75 (MQ=255) ggTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAATATTTTGTTCCACGGAAGATGCAGTGACAAc < 1:401949/76‑1 (MQ=255) ttCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATc < 2:1930001/76‑1 (MQ=255) tGGGTAATGACGCTGATACTGGCAGGTGCATCTTTTAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGttt > 2:38971/1‑76 (MQ=255) gggTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTc < 1:304277/76‑1 (MQ=255) ggTAATGACTCTGATACTGGCAGGCGGGATTTTTAGATTTTTTTCCCCCGCAGATGGAGGGAAAACCCTCGGTTCg > 2:384198/1‑76 (MQ=255) tAATGACGCTGATACTGGCAGGTGCATCTTTATGATTTTGTTTCACGGAAGATGCAGTGACAACCATCGTTTCGcc > 2:227409/1‑76 (MQ=255) aaTGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGATGATGCAGTGACAACCATCGTTTCGCca < 1:1687219/76‑1 (MQ=255) aCGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCcatcat < 1:252001/76‑1 (MQ=255) cGCTGATACTGGCAGGGGGATCTTTTAGATTTTTTTTCACGGGAGATGCAGGGACAAACATCGTTTTGCcctcatc > 2:943500/1‑76 (MQ=255) cGCTGATACTGGCAGGGGCATATTTAAGATTTTTGTTCACCGAAGAGGGAGTGGCAAACAACCTTTCGCcctcatc > 1:1901095/1‑76 (MQ=255) gCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCcatcatca < 2:533749/76‑1 (MQ=255) gATACTGGCAGGGGCATCTTTAAGATTTTTTTTCACGGAAGATGCAGGGACAACCATCGTTTCCCCCTCATCATCg > 2:914904/1‑76 (MQ=255) aCTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGAcc < 1:1848704/76‑1 (MQ=255) cTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCg < 1:939773/76‑1 (MQ=255) ngCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCt > 1:391516/2‑76 (MQ=255) aGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAac < 1:935569/76‑1 (MQ=255) | TTTGAGGACATCCTTCAGATTCTGTACCGGTTTTCGCTGCAGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAGATTTTGTTCCACGGAAGATGCAGTGACAACCATCGTTTCGCCATCATCATCGACCGCTAAC > REL606/2197792‑2197939 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |