Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,176,010 | 0 | T | G | 28.8% | 138.1 / 13.2 | 59 | A37A (GCA→GCC) | fhuE | ferric‑rhodotorulic acid outer membrane transporter |
Reads supporting (aligned to +/- strand): ref base T (4/38); new base G (17/0); total (21/38) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.16e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTGA > REL606/1175935‑1176085 | tACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGt < 1:930145/76‑1 (MQ=255) cGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgc < 2:1073961/76‑1 (MQ=255) cGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgc < 2:2598/76‑1 (MQ=255) cGCTATAATCATTTTCGCCATCATCTGTAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgc < 1:344111/76‑1 (MQ=255) ttAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCGCCGTTTCTTCAGTGGCTGGTgcagca < 1:394993/75‑1 (MQ=255) gTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgcagca < 1:2859/76‑1 (MQ=255) tAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCaa < 2:1473608/76‑1 (MQ=255) tAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCaa > 2:851645/1‑76 (MQ=255) cATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGc < 2:464799/76‑1 (MQ=255) cATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGc < 2:38362/76‑1 (MQ=255) aTTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAACg > 2:1612795/1‑76 (MQ=255) tttCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGc < 1:1519554/76‑1 (MQ=255) cGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAGAGGgaga > 2:1295649/1‑76 (MQ=255) cGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGa > 2:1207466/1‑76 (MQ=255) cGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGa > 1:366017/1‑76 (MQ=255) ccATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAg < 2:1311014/76‑1 (MQ=255) ccATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAg > 1:1640406/1‑76 (MQ=255) catcatCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAGCGGCAGAAgg > 2:641043/1‑76 (MQ=255) atcatcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGt < 1:388537/76‑1 (MQ=255) atcatcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAGCGGCCGAAggg > 2:1129707/1‑75 (MQ=255) tcatcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTa < 1:1421623/76‑1 (MQ=255) tcatcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTa < 1:1146997/76‑1 (MQ=255) catcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGaa > 2:124466/1‑76 (MQ=255) atcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGGACAAGCGGGAGAGGGGAAt > 1:857950/1‑76 (MQ=255) atcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAGCGGCAGAAGGGAAt > 1:799033/1‑76 (MQ=255) atcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAAt > 1:857958/1‑76 (MQ=255) tcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATa < 1:1139380/76‑1 (MQ=255) tcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATa > 2:1486655/1‑76 (MQ=255) cTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAg > 2:1475372/1‑76 (MQ=255) tGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCCGCAAAAACGGCAGAAAGGAATaag > 1:658298/1‑74 (MQ=255) tGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGt > 1:16431/1‑76 (MQ=255) ggAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGCGCAAAAGCGGCCGAAGGGAATAATg > 1:742849/1‑76 (MQ=255) gCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCa < 1:1422606/76‑1 (MQ=255) tgtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGt < 1:8232/76‑1 (MQ=255) tgtgGCTGAACCCTCAACAATCACCATTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGt < 1:1017610/76‑1 (MQ=255) gCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTa < 2:206007/76‑1 (MQ=255) cTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTAt < 1:280984/76‑1 (MQ=255) tGAACCCTCAACAATCACCGTTTCTTCAGTGGTTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATg < 1:1612777/76‑1 (MQ=255) gAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGc < 1:1038314/76‑1 (MQ=255) gAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGc < 1:816823/76‑1 (MQ=255) gAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGc < 1:632557/76‑1 (MQ=255) gAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGc > 2:94170/1‑76 (MQ=255) aaCCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCa < 1:1062365/76‑1 (MQ=255) aaCCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCa > 2:80131/1‑76 (MQ=255) aCCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCaa > 2:1521506/1‑76 (MQ=255) ccTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAAcc > 2:917975/1‑76 (MQ=255) cTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCg < 1:1431391/76‑1 (MQ=255) cTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCg < 1:298398/76‑1 (MQ=255) caacaaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCGGc > 1:1149492/1‑76 (MQ=255) caacaaTCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCCGTGCTATGCAACCGGc > 2:410459/1‑76 (MQ=255) acaaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGGAGAAGGGAATAGTGCCAGTGCTATGCAACCGGCaa > 2:4320/1‑76 (MQ=255) caaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAag < 2:934166/76‑1 (MQ=255) aaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAagt < 2:1185009/76‑1 (MQ=255) tcttcAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAAACGGCAAGTAGTGACgggt > 2:565977/1‑74 (MQ=255) cttcAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGttt < 1:1559388/76‑1 (MQ=255) tcAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTgg < 1:1490722/76‑1 (MQ=255) cAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGt < 1:510838/76‑1 (MQ=255) gTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGa < 2:485050/76‑1 (MQ=255) tGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGAt < 2:1612029/76‑1 (MQ=255) tGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGAt > 2:776208/1‑76 (MQ=255) ggCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATg < 2:1434701/76‑1 (MQ=255) cTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGc > 2:9270/1‑76 (MQ=255) ggTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCtt < 1:208431/76‑1 (MQ=255) gTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTg < 1:1298932/76‑1 (MQ=255) gTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTg < 2:253193/76‑1 (MQ=255) gTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTg < 1:1523079/76‑1 (MQ=255) tGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTGa > 1:944041/1‑76 (MQ=255) | TACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTGA > REL606/1175935‑1176085 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |