Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 40.6% | 29.2 / 11.2 | 32 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (10/9); new base A (0/13); total (10/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.65e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCT > NC_000913/357741‑358032 | tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac > 1:119643/1‑150 (MQ=255) gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct > 2:1001012/1‑151 (MQ=255) aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctct > 1:20022/1‑151 (MQ=255) aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctct > 2:770381/1‑151 (MQ=255) gCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt > 1:509397/1‑151 (MQ=255) ccGCAGCCTTGTGTCAGCGGCTCCCGTGTTTCTTTCCCGCCACCACCCCCGCTTTTAT‑CAACCGCAATGTCCAAAAAGCCTTT‑CCCGGCGGCTGTCGGCCACCTATTTTTCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg < 1:894907/149‑1 (MQ=255) cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg > 2:422954/1‑149 (MQ=255) aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 1:747585/1‑145 (MQ=255) ccTGTGTCAGCGGCTACCGTGATTAATTCCCGCCAACACCCGCGCATTTAC‑CCACCGCCTGGTGAAAAAATGCCTT‑CGCGGCGGCTGTCTCCCAGTTATTTTCCCCCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGc < 1:953171/145‑1 (MQ=21) cAGCGGCTCCCGTGATTCTTTCCCGCCAACACCCGCGATTTTAC‑CCAACCCCATGTGCAAAAATGCCTC‑CGGGGCGGCTTTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg < 1:333024/149‑1 (MQ=255) cAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGaaaaa > 1:674295/1‑146 (MQ=255) gCGGCTCCCGTGTTTCATTCCCCCCAACACCCGCGATTTTAC‑CCACCGCCTTTTGAAAAAATGCCTC‑CGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag < 1:996254/141‑1 (MQ=255) cGGCTCCCGTTTTTTTTTCCCCCCAACACCCGCGCTTTTCC‑CCACCGCCTTGTGAAAAAATGCTTC‑CGCGGGGGTTGTCTGCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATa < 2:898688/145‑1 (MQ=255) ggCTACGGTATTTCATTCCCGCCAACACCCGCCCATTTAC‑CCAACGCCATGTGAAAAAATGCCTC‑CGCGGCGGCTGCCTGCCAGTTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCt < 1:851906/147‑1 (MQ=255) gcctcccGTTATTCTTTCCCGCAACCACCCGC‑CATTTATCCCACCGCCATGTGAAAAAATGCCTT‑CGCGGGGGCTTTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCt < 2:612676/142‑1 (MQ=255) cGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaa > 2:588131/1‑135 (MQ=255) tgtttCTTTCCCCCCAACACCCGCCCTTTTAC‑CCACCCCCAGGTGAAAAAATCCTTT‑CGCGGGGGCTGTCGGCCGGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTc < 2:881362/147‑1 (MQ=255) tGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGgaaaaacatttatac > 2:801386/1‑134 (MQ=25) gtttCATTCCCCCCAACACCCGCGCATTTAT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa < 2:661353/148‑1 (MQ=255) ttCTTTCCCCCCAACACCCGCGTTTTTAC‑CCACCGCCATGTGCAAAAAGCCCTC‑CGCGGCGGCTGTCGGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTc < 2:283341/147‑1 (MQ=255) ttcccGCCAACACCCGCGCTTTTAT‑CCAACCCCATGTG‑AAAAATGGCTTCCGCGGCGGCTGTCGGCCACCTTTTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttgt < 1:892647/147‑1 (MQ=255) ccAACGCCATGTG‑AAAAATTCCTTGCGGGGCGGTTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAt < 2:712064/149‑1 (MQ=255) cgcattgtaCAAAAATCCCTT‑CGCGGCGGCTGTCTCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg < 1:472818/144‑1 (MQ=255) aaaacTGCCCT‑CGCGGCGTCTCTATGCCAGTTATTTCGCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa < 1:232315/146‑1 (MQ=255) aaaaTCCTTT‑CGGGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg < 1:996409/151‑1 (MQ=255) cTT‑CGGGCCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGa < 2:788290/148‑1 (MQ=255) tt‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAATTTTAAACCGTCACGTTGTTGGGGGAACGCCGCTGGGGAAAAAATGggaagggaagggggggg > 1:531121/1‑139 (MQ=255) tcgcccgctttttttcCCGCGGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 2:412223/138‑1 (MQ=255) cTGCCAGCTTTTTTTCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 1:50790/148‑1 (MQ=255) ccggCTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 1:465887/142‑1 (MQ=255) agtttttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa < 1:824212/145‑1 (MQ=255) tttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt < 2:673685/150‑1 (MQ=255) cgccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCt < 1:480016/151‑1 (MQ=255) | TGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTAT‑CCAACGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCT > NC_000913/357741‑358032 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |