Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 821,440 | 0 | C | T | 34.6% | 23.6 / 14.4 | 26 | F216F (TTC→TTT) | ybhM | Bax1‑I family protein YbhM |
Reads supporting (aligned to +/- strand): ref base C (4/13); new base T (9/0); total (13/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.58e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGTTGTGGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAG > NC_000913/821373‑821506 | agTTGTGGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGcttc < 1:2206792/68‑1 (MQ=255) tgtgGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCttt < 1:997163/68‑1 (MQ=255) ggAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCttttttttg > 1:2963950/1‑65 (MQ=255) ggAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTtgg > 1:1116959/1‑68 (MQ=255) ggAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTtgg > 1:857755/1‑68 (MQ=255) ggAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTtgg > 1:2400213/1‑68 (MQ=255) aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAAGACTCTATTACTATTTTATCGGCTTGTTTtgggt > 1:1727917/1‑66 (MQ=255) aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTtgggt > 1:2353997/1‑66 (MQ=255) aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTTGGAt > 1:3017413/1‑68 (MQ=255) ccGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTTGGGTatt > 1:3435969/1‑66 (MQ=255) gTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTTGGGTAAtttt > 1:3536612/1‑66 (MQ=255) tCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTGGGTTAATTGCcg > 1:1545352/1‑68 (MQ=255) gcTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGc < 1:2259719/68‑1 (MQ=255) gcTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGc < 1:97417/68‑1 (MQ=255) gcTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGc < 1:2837355/68‑1 (MQ=255) gcTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGc < 1:2843581/68‑1 (MQ=255) cTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt < 1:1662134/68‑1 (MQ=255) cTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt < 1:754298/68‑1 (MQ=255) cTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt < 1:3279498/68‑1 (MQ=255) cTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt < 1:1406486/68‑1 (MQ=255) cTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt < 1:2886589/68‑1 (MQ=255) cTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt < 1:1769571/68‑1 (MQ=255) cAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACg < 1:1365959/68‑1 (MQ=255) gTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTTGGATACTTTCCGCGATAGCTATAACGCtt > 1:1794523/1‑68 (MQ=255) gTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTTGGATACTTGCCGCGATAGCTATAACGCtt > 1:2962169/1‑68 (MQ=255) gTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGGTACTTTCCGCCCCAGCTATAACGGtt > 1:335013/1‑68 (MQ=255) tGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTg < 1:676187/68‑1 (MQ=255) tAACACTCTATTACTATTTTATCGGCTTTTTTTGGATACTTTCCCCGATAACTATAAAGCTTGTCTgg > 1:2906532/1‑68 (MQ=255) tAACACTCTATTACTATTTCATCGGCTTTTTTTGGGTACTTTCCGCGATAGCTATAACGCTTGTCTgg > 1:550411/1‑68 (MQ=255) aaCACTCTATTAATATTTTATCGGCTTTTTTGGTATACCTGTCGCGATAGCTATTAAGCATGTCTggg > 1:2192239/1‑67 (MQ=255) cacTCTATTACTATTTTATCGGCTTTTTTTGGGTACTTGCCGCGCTAGCTATTACGCTTGTCTTGCaa > 1:1727890/1‑68 (MQ=255) acTCTATTACTATTTTATCGGCTTTTTTTGGATACTTTCCGCGATAGCTATAACGCTTGTCTGGCaaa > 1:647956/1‑68 (MQ=255) acTCTATTACTATTTTATCGGCTTCTTTTGGTTTGTTTCCGCCGTAGCTATAAAGCTTTTCTGGCaaa > 1:3239775/1‑68 (MQ=255) acTCTATTACTATTTTATCGGCTTCTTTTGGGTACTTTCCGCGCTAGCTATAACGCTTGTCTGGCaaa > 1:1008464/1‑68 (MQ=255) tctATTACTATTTTATCGGCTTTTTTTGGGTACTTGCCGCGCTAGCTATAACGCTTGTCTGGCAAAgg > 1:499148/1‑68 (MQ=255) ttttATCGGCTTTTTTTGGGTACCTTCCGCGATAACTTTAACCCTTGTTTGGCAAAAGGATAcgcggt > 1:1533156/1‑66 (MQ=255) tcttTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAg < 1:704393/68‑1 (MQ=255) | AGTTGTGGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTTTTTCATTAG > NC_000913/821373‑821506 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |