Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,176,010 | 0 | T | G | 32.6% | 102.1 / 11.0 | 49 | A37A (GCA→GCC) | fhuE | ferric‑rhodotorulic acid outer membrane transporter |
Reads supporting (aligned to +/- strand): ref base T (9/24); new base G (16/0); total (25/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.30e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCT > REL606/1175936‑1176082 | aCGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGt < 2:101243/75‑1 (MQ=255) aCGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGt < 2:1381566/75‑1 (MQ=255) cGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTg < 2:1364665/75‑1 (MQ=255) cTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGgcg > 1:687399/1‑74 (MQ=255) tctAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgcag < 1:1474052/73‑1 (MQ=255) gTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGggcggc > 2:938444/1‑75 (MQ=255) gTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgcagc < 2:1021755/75‑1 (MQ=255) gTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTgcagc < 1:313262/75‑1 (MQ=255) gTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGggagc > 1:182528/1‑75 (MQ=255) tAATCATTTTCGCCATCATCTGGAGCTGTGCCTGAACCCTCAACAATCACCGCTTCTTCAGTGGCTGGggcggcc > 2:1424674/1‑74 (MQ=255) aaTCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCaa < 2:658877/75‑1 (MQ=255) aaTCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCaa < 1:569855/75‑1 (MQ=255) tCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGGGCaaaa > 2:554970/1‑75 (MQ=255) ccATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCCAAAGCGGCAGaa > 2:334915/1‑75 (MQ=255) catcatCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAg < 1:607953/75‑1 (MQ=255) catcatCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAg > 1:1569062/1‑75 (MQ=255) tcatcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAGCGGCAGAAggg > 2:711844/1‑74 (MQ=255) tcatcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGGAGAAggg > 2:792920/1‑74 (MQ=255) atcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGCGCAAAAACGGCAGAAGGGaa > 1:1288436/1‑75 (MQ=255) atcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAACGGCAGAAGGGaa > 1:645122/1‑75 (MQ=255) tcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAAt < 2:1127172/75‑1 (MQ=255) tcTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAGCGGCAGAAGGGAAt > 1:864059/1‑75 (MQ=255) tGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGGAGCAAAAGCGGCAGAAGGGAAAAg > 1:1245778/1‑75 (MQ=255) ggAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAACGGCAGAAGGGAAAAgg > 1:485455/1‑74 (MQ=255) ggAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAACGGCAGAAGGGAAAAgg > 2:858523/1‑74 (MQ=255) aGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGc < 1:12750/75‑1 (MQ=255) tgtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGGAGGGAATAGTGCCAg > 1:1473977/1‑75 (MQ=255) tgtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAg < 2:1266265/75‑1 (MQ=255) tgtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAg > 1:225060/1‑75 (MQ=255) gtgGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAAAAGTGCCAgg > 1:1082276/1‑74 (MQ=255) cTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTa < 2:258678/75‑1 (MQ=255) cTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTa > 2:282166/1‑75 (MQ=255) aaCCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGc < 2:1054437/75‑1 (MQ=255) cTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAAcc > 2:33147/1‑75 (MQ=255) aacaaTCGCCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGc < 2:439389/75‑1 (MQ=255) acaaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCa < 2:475890/75‑1 (MQ=255) acaaTCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCGGCa > 1:1150605/1‑75 (MQ=255) acaaTCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCGGCa > 1:1438923/1‑75 (MQ=255) aCCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAagtagt < 1:378979/75‑1 (MQ=255) ccGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAg < 1:712001/52‑1 (MQ=255) ccGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAg > 2:712001/1‑52 (MQ=255) cGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGa < 1:1402548/75‑1 (MQ=255) gTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGAc < 1:1370970/75‑1 (MQ=255) cttcAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGtt < 2:395267/75‑1 (MQ=255) cttcAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGtt < 2:616206/75‑1 (MQ=255) cttcAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACgggt > 2:722523/1‑73 (MQ=255) tcAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTg < 1:1406253/75‑1 (MQ=255) tcAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGGTTg > 2:1471745/1‑75 (MQ=255) aGTGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCCGTGCTATGCAACCGGCAAGGAGTGACGGGTTTGt > 2:944133/1‑75 (MQ=255) tGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGa < 2:389134/75‑1 (MQ=255) tGGCTGGGGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCGGCAAGGAGTGACGGGTTGGGGa > 1:1269843/1‑75 (MQ=255) ggCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGAt > 2:933782/1‑75 (MQ=255) ggCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGGTTGGTGAt > 2:593171/1‑75 (MQ=255) gCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATg > 1:1461835/1‑75 (MQ=255) gCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATg > 1:271578/1‑75 (MQ=255) gCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATg < 2:715537/75‑1 (MQ=255) gCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATg < 1:473614/75‑1 (MQ=255) tGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGc > 1:794465/1‑75 (MQ=255) ggTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCt < 1:621502/75‑1 (MQ=255) | ACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCT > REL606/1175936‑1176082 |
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: ‑ |