| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,842 | 0 | C | A | 37.1% | 51.1 / 13.0 | 35 | L283F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base C (8/14); new base A (0/13); total (8/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.52e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGA > NC_000913/357696‑357988 | gTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATGCATGTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑Caaa < 1:455048/150‑1 (MQ=255) gAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑Caaaaa < 1:410865/150‑1 (MQ=40) gAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAAt < 2:553476/151‑1 (MQ=40) aaCGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATg < 1:512414/151‑1 (MQ=40) ccTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCtt < 2:363612/151‑1 (MQ=40) ccTTATCCGCCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCtt > 1:293181/1‑151 (MQ=255) aCGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGt < 1:533121/150‑1 (MQ=255) aCGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGt < 2:434754/150‑1 (MQ=255) 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) | GTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGA > NC_000913/357696‑357988 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |