Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,176,010 | 0 | T | G | 35.6% | 84.5 / 11.5 | 45 | A37A (GCA→GCC) | fhuE | ferric‑rhodotorulic acid outer membrane transporter |
Reads supporting (aligned to +/- strand): ref base T (4/25); new base G (16/0); total (20/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.49e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCT > REL606/1175937‑1176082 | cGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGgc > 2:173356/1‑76 (MQ=255) cGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGgc > 1:1286701/1‑76 (MQ=255) gCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgca < 1:486864/76‑1 (MQ=255) cTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGgcag > 1:1343375/1‑76 (MQ=255) tGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGgcagg > 1:560201/1‑75 (MQ=255) tGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGgcagc > 2:575804/1‑76 (MQ=255) gTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGgcggca > 2:931378/1‑76 (MQ=255) gTAATCATTTTCGCCATCATCTGGAGCTGCGGCTGAGCCCTCAACAATCACCGTTTCNTCAGTGGCTGGTgcagca < 1:1535580/76‑1 (MQ=255) tAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCaa > 1:491910/1‑76 (MQ=255) aTCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCaaaa < 2:117848/76‑1 (MQ=255) aTCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCaaaa < 2:1139637/76‑1 (MQ=255) tCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAg < 1:965718/76‑1 (MQ=255) cATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGc < 1:829543/76‑1 (MQ=255) tttCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAACGGc > 2:641871/1‑76 (MQ=255) tttCGCCATCATCTGGAGCNGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCAGGGGCAGCAAAAGCGGc > 1:1507979/1‑76 (MQ=255) ttCGCCATCATCTGGAGCTGTGGATGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAACGGCa > 2:1304091/1‑76 (MQ=255) ttCGCCATCATCTGGAGCTGTGGATGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGGAAAAGCggga > 2:1304108/1‑74 (MQ=255) tCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAg < 1:519743/76‑1 (MQ=255) gcgcatcatCTGGNGCTGTGGCTTCACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGa < 1:212833/73‑1 (MQ=255) gCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAANGCGGCAGaa < 1:40688/76‑1 (MQ=255) ccATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAGCGGCAGGAg > 1:879158/1‑76 (MQ=255) catcatCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAgg > 1:333152/1‑76 (MQ=255) catcTGGAGCTGTGGCTGNACCCTCAACAATCACCGTTTCTTCAGTAGCTGGGGCAGCAAAAGCGGCCGAAGGGaa > 1:1549471/1‑76 (MQ=255) catcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGaa > 2:118401/1‑76 (MQ=255) catcTGGAGCTGTGGCTGAACCCTCAACAATAACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTaa < 1:840207/76‑1 (MQ=255) atcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAAt < 1:252914/76‑1 (MQ=255) gAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAGCGGCGGAAGGGAATAGTGc > 1:344220/1‑76 (MQ=255) gCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCa < 1:252963/76‑1 (MQ=255) gtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTg < 2:1496671/76‑1 (MQ=255) tgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGc > 2:720968/1‑76 (MQ=255) ggCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGGAGGGAAAAGTGCCCGCGCt > 1:606708/1‑76 (MQ=255) cTGATCCCTCAACAATCACCGTTTCTTCAGTGGCTGATGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTAt < 1:1585726/76‑1 (MQ=255) cTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATANTGCCAGTGCTAt < 1:15141/76‑1 (MQ=255) cTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCCGTGCTAt > 2:878488/1‑76 (MQ=255) aaCCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCa < 2:198355/76‑1 (MQ=255) aaCCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCa < 1:308740/76‑1 (MQ=255) cccTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAAc < 1:1236621/76‑1 (MQ=255) ccTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAAcc < 1:243648/76‑1 (MQ=255) caacaaTCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCGGc > 2:91897/1‑76 (MQ=255) aacaaTCACCGTTTCTTCAGTGGCTGGTGCAACAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCa < 2:250698/76‑1 (MQ=255) aacaaTCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCGGCa > 1:812462/1‑76 (MQ=255) aacaaTCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCGGCa > 2:1072928/1‑76 (MQ=255) caaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGGAATAATGCCAGTGCTATTCAACCGGCCag > 2:1014633/1‑76 (MQ=255) tCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAagtag < 1:711808/76‑1 (MQ=255) cACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAagtagt < 2:430151/76‑1 (MQ=255) cACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGGGCTATGCAACCGGCAagtagg > 2:909587/1‑75 (MQ=255) ccGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGa < 1:1175015/76‑1 (MQ=255) tttCTTCAGTGGCTGGGGCGGCAAAAGCGGCAGCAGGTGATATTGCCCATGCTTTGCAACCCGCCAATCGTGaacg > 2:1305909/1‑73 (MQ=255) tttCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACgg > 2:1305889/1‑76 (MQ=255) tcAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTgg < 2:455476/76‑1 (MQ=255) cAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGGGGCAGGGCTATTCAACCGGGCAGGAGTGGCGGGTTggg > 2:1549662/1‑75 (MQ=255) cAGTGGCTGGGGAAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTggg > 2:955083/1‑75 (MQ=255) tGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGGGAt > 1:762851/1‑76 (MQ=255) ggCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATg < 1:167896/76‑1 (MQ=255) tGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCt > 1:1297808/1‑76 (MQ=255) tGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCt > 1:1021527/1‑76 (MQ=255) | CGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCT > REL606/1175937‑1176082 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |