| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,890 | 0 | C | A | 44.4% | 26.3 / 10.4 | 27 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base C (10/5); new base A (0/12); total (10/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.64e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG > NC_000913/357744‑358023 | tAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaacaa > 1:787609/1‑149 (MQ=255) gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAccg < 1:14162/151‑1 (MQ=255) cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTcc > 1:449837/1‑151 (MQ=255) aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTcc > 2:588037/1‑150 (MQ=255) attcccAGCCGCGAACTTGTGTAGGCGGCTCCCGTTATTCATTCCCCCCACCACCCGCGCATTTACCAACCCCCGTGTAAAAAAATCCCTCC‑GGGGGGGCT‑TCCGCCCGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagta < 1:408698/146‑1 (MQ=255) ttCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTg > 2:143974/1‑151 (MQ=255) cAGACCTGTTTCCGGGGCTCCCGTTTTTCTTTCCCCCCAACACCCGCCCTTTTACCCACCGCCATGTGAAAAAATCCTTTC‑GGGGCGCCTGTTTGCCAGTTTTTTTC‑CCCCCGCAACAAACCCGCCGTTCTCTCCAGTAGTg < 2:356698/142‑1 (MQ=255) cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg > 2:134824/1‑148 (MQ=255) ccTGTGCCGGCGCCCCCCGTGTTTCTTTCCCCCCAACACCCGCCCATTTACCCAACGCCTTTTGCAAAAATGCTTCC‑GGGGGGGCTTTCTCCCAGTTTTTTTT‑CCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAaga < 2:583896/149‑1 (MQ=255) ggCTACCGTGTTTCTTTCCCGCCAACACCCGCGATTTTACCCAACGCCATGTGAAAAAATGCTTTC‑GCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAaga < 2:251349/138‑1 (MQ=255) aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGgaaaaaggt > 1:137068/1‑138 (MQ=255) cggtaaTTCTTTCCCCCCAACACCCCCGCTTTTACCCCACCCCATGTAAAAAATCGCCTCC‑GGGGCGGCTGTCCGCAACTTTTTTTC‑CCGCACAAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCg < 1:453983/146‑1 (MQ=255) gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGAGaaaacgatttagacc > 1:831579/1‑136 (MQ=25) aacaccccccattttACCCACCCCCATGTGCAAAAGTGCCTGC‑GGGGGGGTTGTCTCCCAGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc < 2:221086/139‑1 (MQ=255) caacCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCGAGAGAAAAAAATTTTAACCCCTCCCGTTTTGTgggg > 1:749801/1‑144 (MQ=25) acccccatttttccccACCGCCATGTG‑AAAAATGCCTTCGGGGGCGGTTGTCTGCCAGTTTTTTTT‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc < 1:576338/134‑1 (MQ=255) cccccgaTTTCTCCAACCCCGTGTGCAAAATTCCCCTC‑GGGGCGGTTGTCCCCCTCTTTTTTTC‑CCCCCTCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt < 2:538368/139‑1 (MQ=255) ccGCTCTTTTACCAAACCCCATGTGAAAAAATCCCTCC‑GCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg < 2:698637/147‑1 (MQ=255) cATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGAAAAAAATTTAAACACCTCATGTTTTGTGGGGATCgc > 2:625742/1‑143 (MQ=37) ccaaccccgtgaaaaaaaaTGCCTCC‑GGGGCGGCTCTCGCCCAGCTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAAtg < 1:788210/136‑1 (MQ=255) caccGCCATGTAAAAAAATTCCTGC‑GCGGCGTTTGTCTCCCAGTTATTTTC‑CCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAAtg‑tg < 2:12388/148‑1 (MQ=255) aaTCCCTCC‑GGGGCGGCTGTCTGCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATG‑TGGAAAGGGAAGTgcggc < 2:686928/151‑1 (MQ=255) cTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAGAAGGAAAAAATTTAAACACGTCATATGTTGTGGGAATAGCTGC‑TCGGGACAAAGAGGGAAAGGGag > 2:679201/1‑139 (MQ=255) ggcgtcTCTCTGCCGTTTTTTTTT‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCt < 2:281501/149‑1 (MQ=255) ccattttttttc‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgcg < 2:754415/141‑1 (MQ=255) gCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAAAAAAATTTAACCCCTCCTCTTTTTTGGCGATCCCTGCTTGGTAAAAAAT‑GGGAAAGGG‑AGTGGGGCGGATCAACTCCCAAACCCCGCCaa > 2:226734/1‑149 (MQ=255) ttttg‑CCGCCGAAACCAAGCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTcc < 1:725895/132‑1 (MQ=255) ttttc‑CCGCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg < 1:69390/148‑1 (MQ=255) | TAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTC‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGC‑TGGTAACAATG‑TGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG > NC_000913/357744‑358023 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |