Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,340,075 | 0 | A | T | 31.7% | 38.5 / 25.8 | 45 | intergenic (+60/‑271) | ribAB‑II/PP_3814 | 3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein |
Reads supporting (aligned to +/- strand): ref base A (12/16); new base T (7/6); total (19/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.37e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGAGGGTGACACCACCGACGCATTCGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTG > NC_002947/4339929‑4340223 | gggCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGAgg < 6:34117/149‑1 (MQ=255) gggCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGAgg < 5:161188/149‑1 (MQ=255) gggCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGAgg < 5:147786/149‑1 (MQ=255) cgccgcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGAGGGTGa < 6:360564/149‑1 (MQ=255) cgcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGTCGGTggtgt > 7:167559/1‑145 (MQ=255) cgcTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGAGGGTGacac < 4:229230/149‑1 (MQ=255) tGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGTCGGTggtgtcac > 8:320788/1‑142 (MQ=255) tACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGAGGGTGACACCACCGACg < 3:232480/149‑1 (MQ=255) ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGTTATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc < 1:313695/149‑1 (MQ=25) ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc < 8:167559/149‑1 (MQ=25) cTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc > 6:300433/1‑139 (MQ=25) cTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc < 5:300433/139‑1 (MQ=25) tGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc < 6:246887/136‑1 (MQ=25) tGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGTCGGTGGTGTCACCCTCGCATTcgc > 5:246887/1‑136 (MQ=25) tGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCCGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGAGGGTGACACCACCGACGCATTcgc‑g < 3:10259/137‑1 (MQ=255) tGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCCGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGAGGGTGACACCACCGACGCATTcgc‑g > 4:10259/1‑137 (MQ=255) aGCTGGCGGGGGGGCAGAGCATCCCGTGCGACCCGGGTTTGGATGGCTGATAACACACAAGCAATTGATTTAGGGGGGCCCTGTTGGAGGGGGTTTACCC‑GCGAATGCGAGGGTGACACCACCTACGCATTCGC‑GGGT‑‑AAACCCGCTc < 2:88243/148‑1 (MQ=37) cATCCCGTGCGACCGGGGGATGTATGGCTGAGAAAACACACGAAATTTATTTTGCGGGGCCCTGTAGGAGGGGGTTTACCC‑GCTAAGGCGAGGGTGACACCACCGACGCTTTCGC‑GGTT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCtt < 2:291259/149‑1 (MQ=25) tCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGAGGGTGACACCACCGACGCATTAGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCtttt > 1:300460/1‑149 (MQ=255) ccGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGAGGGTGTCACCCCCGCCGCATTAGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGa > 1:24356/1‑149 (MQ=255) cGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑CCGAATGCGAGGGGGACACCACCCCCGCCTTAGC‑GGGT‑‑AAACCCGCTCCTACACGGCCCGAGGTGCGTGTGACAggg > 8:283493/1‑149 (MQ=25) gACCGGGGTATGGATGGCTGATAACACACAAGCAATTTTTTTCGCGGGGCCCTGTAGGAGCGGTTTTACCC‑GCTAATGCGAGGGTTACACCACCGACGCATTCGC‑GGTT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAggg < 2:24356/148‑1 (MQ=37) aCCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTACC‑CCCAAAGCGAGGGTGTGCCCCCCGCCCCAATAGC‑GGGG‑‑AAACCCCCGGCGAAACGGCCCTTGGGGGTATTTACACGGCg > 6:160442/1‑149 (MQ=255) ggatggCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTCCC‑CCGAAATCGAGGGTGACAACACCGCCCCATTAGCGGGGT‑‑AAACCC‑CACCTACAGGGCCCATGGTGCTTTTGACAgg > 1:73706/1‑136 (MQ=255) gatggCTGACAACACACACCCAATTGATTTCGCGGGGCCCTGTCGGCGGGGGTTTCCCC‑GCTAATGCGCGGGTTACACCACCGACGCATTCGG‑GGTT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGc < 2:300460/148‑1 (MQ=37) acaacaaacaAGCATTTTTTTTCGCGGCGCCGTGT‑GGGGGGGGTTTACCCGGCTATGGCGCGGTTCACCCCCCCGCCG‑ATT‑GC‑GGGTTAAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGccc < 8:173545/147‑1 (MQ=255) acacaAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGAGGGTGACACCACCGACGCATTCGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATAc > 8:283039/1‑149 (MQ=255) aaGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGAGGGTGACACCACCGACGCATTCGC‑GGGT‑‑AAACCCGc > 5:357792/1‑88 (MQ=255) aaGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGAGGGTGACACCACCGACGCATTCGC‑GGGT‑‑AAACCCGc < 6:357792/88‑1 (MQ=255) gCAATTTATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCTAATGCTAGGGTGACACCACCTACGCATGCGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTg < 7:299668/149‑1 (MQ=38) aTTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGAGGGTGACACCACCGACGCATTCGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTg < 5:362005/117‑1 (MQ=255) gCGGGGCCCTGTAGGCGCGGGTTAACCC‑GCGAATGCGAGGGTGACACCACCGACGCATTCGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg < 7:283493/148‑1 (MQ=255) gCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGAGGGTGACACCACCGACGCATTCGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg < 7:283039/148‑1 (MQ=255) gCGGGCCCCTGTTGGTGGGGTTATACCG‑GGTAATGCGCGGGTGACACCACCGCCTCTTTCGG‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg < 1:311827/148‑1 (MQ=255) gCGGGCCCCTGTCGGTGCGGTTTTACCC‑GCTAATCCGCGGGTGACACCACCGACGCATTCGC‑GGTT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg < 5:12948/148‑1 (MQ=38) cggcgCGGGTTTACCC‑GCTAATGCGTCGGTGGTGCCACCCTCGCATTCGC‑GGTT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGaa < 2:298490/135‑1 (MQ=21) aGGAGCGGGTTTACCC‑GCGAATGCGTCGGTGGTGTCACCCTCGCATTCGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGGTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGaa > 1:298490/1‑139 (MQ=18) gTTTACCC‑GCGAATGCGAGGGTGACACCACCGACGCATTCGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATcc > 3:241188/1‑149 (MQ=255) cc‑GCGAATGCGTCGGTGGTGTCACCCTCGCATTCGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa > 3:339078/1‑149 (MQ=25) cc‑GCGAATGCGAGGGTGACACCACCGACGCATTCGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa > 1:76081/1‑149 (MQ=255) cgtaattGCGTGGGTCACACCACCGACGCATTGGG‑GTGT‑‑AACCCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc < 5:207819/143‑1 (MQ=255) gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc > 8:213613/1‑149 (MQ=25) gcgAATGCGAGGGTGACACCACCGACGCATTCGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc > 5:349898/1‑149 (MQ=255) gcgAATGCGAGGGTGACACCACCGACGCATTCGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc > 2:148170/1‑149 (MQ=255) cgAATGCGAGGGTGACACCACCGACGCATTCGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt < 2:76081/149‑1 (MQ=255) aGGGTGACACCACCGACGCATTCGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTg > 3:351310/1‑149 (MQ=255) | GGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTACCC‑GCGAATGCGAGGGTGACACCACCGACGCATTCGC‑GGGT‑‑AAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTG > NC_002947/4339929‑4340223 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |