Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913_3_vch_pgi | 860,154 | 0 | A | T | 37.5% | 4.5 / 14.5 | 16 | intergenic (‑126/+20) | ybiV/ybiW | sugar phosphatase; fructose‑1‑P/ribose‑5‑P/glucose‑6‑P phosphatase/putative pyruvate formate lyase |
Reads supporting (aligned to +/- strand): ref base A (7/3); new base T (2/4); total (9/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.02e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCAT > NC_000913_3_vch_pgi/860013‑860301 | tAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACgccgcc > 2:320001/1‑149 (MQ=255) tAACTTTTACGCTCATAGGAGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGAcgacgct < 1:372000/149‑2 (MQ=255) aCTTTTACGCTCATAGGTGTCAGCCTGGCGGGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCgg < 1:8006/149‑1 (MQ=255) aCTTTTACGCTCATAAGTGTTTTCCTAGCTCTAAAATGAAACCTTAGCAGTTTAATAAAAATTTTACTGAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATTTCAATGATTTAGCCCGGCGACGACGCCGCCgg < 1:127713/149‑1 (MQ=255) ttttACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGc > 1:126399/1‑149 (MQ=255) gCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGctaaa > 1:161393/1‑143 (MQ=255) gCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGctaaa < 2:161393/147‑5 (MQ=255) cTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCGGGctaaatca < 1:319321/149‑8 (MQ=255) cTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCCGAGCACAAATTATATTTCGTATGAAAGTAAGGATGAAATTCATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGCCCGAGGAGCTTACAGCAGATGtt > 2:204119/1‑149 (MQ=255) ggCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCa < 1:320001/149‑1 (MQ=255) taaCAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGAAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGtt > 2:263096/1‑148 (MQ=255) aTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCgg < 1:368190/59‑1 (MQ=37) aTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGGCGTCGTCGCCgg > 2:368190/1‑59 (MQ=37) aTGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGAGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCcg > 2:331115/1‑149 (MQ=255) cGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGGGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGGCGGCGACACGCACTACCAGATCGCGATACTGAGCGGCGTGtttg > 1:68584/1‑148 (MQ=255) cGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGACCGGGATGtttt > 1:107959/1‑149 (MQ=255) gacgCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCAt > 1:351610/1‑149 (MQ=255) | TAACTTTTACGCTCATAGGTGTGCTCCTGGCTCGAAAATGAAACCGTAACAGTGTAATAACAATGTGACGCAGAGCACAAATTATATTTCGAATGAAAGTAAGGATGAAATTGATGATGTGAATGATTTAGCCCGGCGACGACGCCGCCGGGCCGAGGAGATTACAGCATATGTTCAGTACGGGCGATGATATCGTCCTGAGCGTCTGGAGAGAGCGCGGTGAAGAACGCGGAATAGCCCGCGACACGCACTACCAGATCGCGATACTGATCGGGATGTTTTTTCGCAT > NC_000913_3_vch_pgi/860013‑860301 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |