Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,873 | 0 | C | T | 37.9% | 26.5 / 11.3 | 30 | S273N (AGC→AAC) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (7/11); new base T (0/11); total (8/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.59e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCG > NC_000913/357724‑358014 | cACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACTGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCt < 1:258469/151‑1 (MQ=255) gAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTAt < 2:393915/149‑1 (MQ=255) gAACTCGGAGGCCTGCTAAGCGTAGCGCATCAGCCCATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTa < 1:282063/148‑1 (MQ=255) cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTc < 1:309652/149‑1 (MQ=255) ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc > 2:391089/1‑151 (MQ=255) aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc < 1:661220/151‑1 (MQ=255) cttcgcacATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc < 1:40172/142‑1 (MQ=255) cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc > 1:75540/1‑151 (MQ=255) actcccccccGCAGACTTCGGTCAGCGCCTCCCGGGTTTCTTTACCCCACACACCCGCG‑TTTTTTCCAACCCCCTGGGAAAAAAATCCCTT‑CGGGGGGGTTGCCTCCCAGTTTTTTTTCCGCGCAAACAAAACCGCCGTTCTCTCCagt < 1:197917/141‑1 (MQ=255) aaTCCCAGCCGCAGCCCTGTGTCAGCGGCTACCGTGTTTCTTTCCCGCCACCACCCGCGCTTTTACCCAACGCCATTGG‑CAAAAATCCCTC‑CGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCagta < 2:561137/150‑1 (MQ=255) ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTTATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGATATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg > 1:617442/1‑151 (MQ=255) ctgacttgtacaaGCGGCTACCGGGATCCATTCCCGCCAACACCCGCGCTTTTATCCAACGCCATGTA‑AAAAACTGCCTG‑CGCGGCGGCTGTCTGCGAGCTATTGTGCCGCCGCAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCa < 1:147332/138‑1 (MQ=255) gCGGCTACCGTGATTCATTCCCCCCAACAACCGCGCTTTTACCCACCGCCATGTC‑AAAAAATGCCTC‑CGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt < 2:124995/151‑1 (MQ=255) ggCTCCCGTGTTTCTTTCCCGCCAACACCCGCCCTTTTACCCAACCCCATGT‑‑AAAAAAAGCCTTGCGCGGGGCCTGTCTCCCAGTTTTTTTTCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTa < 2:9295/150‑1 (MQ=255) ggCTCCCGTGTTTCTTTCCCCCAAACACCCGCGCATTTACCCACCCCCTTGTG‑AAAAAATGCCTT‑CGCGGCGGCTTTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa < 2:557842/151‑1 (MQ=255) cTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaaaaaattttaa > 2:28866/1‑140 (MQ=255) ccgtatttCATTCCCCCCAACACCCCCGCTTTTACCCACCCCCATGTG‑AAAAAATGCTTT‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc < 2:298204/143‑1 (MQ=255) ttCTTTCCCGCCAACACCCGCGCTTTTACCCACCCCCATGTC‑AAAAAATCCCTC‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTc < 2:312514/147‑1 (MQ=255) ttCATTCCCCCCAACACCCGCGCTTTTACCCACCGCCATGTG‑AAAAAATGCCTC‑CGGGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtg < 2:152515/150‑1 (MQ=255) aTTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAAAAACATTTAACCCGTCCtgttgt > 2:47267/1‑151 (MQ=37) tcccGCCACCACCCGCGCATTTACCCACCCCTGTGTG‑AAAAAATGCTTT‑GCCGGGGCCTGTCTCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgtt < 2:536539/145‑1 (MQ=255) cccgcctttttACCCACCCCCTTGTG‑AAAAAATGCCTC‑GGCGGCGGTTGTCTCCCAGTTATTTTCCCGCCCCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg < 2:513376/143‑1 (MQ=255) ccGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACATTTAAACCGTCATGTTTTGTggggaac > 2:206246/1‑141 (MQ=25) ccGCCGTGGG‑AAAAAATGCCTT‑CGGGGCGGTTGTCTCCCTTTTTTTTCGCCGCCAAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa < 1:126716/150‑1 (MQ=255) ccttgtg‑AAAAAATGCCTC‑CGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaa < 2:516042/146‑1 (MQ=255) tgtg‑AAAAAATGCCTT‑CGGGGCGGTTGTCTGCCAGCTTTTTTCCCCCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg < 1:215525/150‑1 (MQ=255) g‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAGGGGAAAAAAAATTTATACCCCCCTCTTGTGTGGGGATCGCTGCTGGGGACAATATGGGAggggg > 1:598561/1‑147 (MQ=255) gcggcgGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa < 2:167817/148‑1 (MQ=255) gTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGGGAAAAAGATTTTATCCCGCCCTGTTGTGTGGGGAATGCCGCCGGTGAAAAAATGGGAAGGGGAGGGGGGCGGATCAGCTCCCGaa > 2:258974/1‑148 (MQ=255) aGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg < 1:163822/144‑1 (MQ=255) | CACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCG > NC_000913/357724‑358014 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |