Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 42.3% | 13.3 / 11.3 | 26 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (8/7); new base A (0/11); total (8/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.39e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01 |
AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCC > NC_000913/357745‑358031 | 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) cgccgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAAAAAATTTAAGACCCGCCTCTTTTTTGGGAACTGCTGCTGTTTACAATGTGGAAAGGGGAGAGGGGGGGAGAAGCTCCCGAACCGCGCCCATTTAATTCgtc > 1:64093/1‑148 (MQ=255) | AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTG‑CAAAAATGCCTT‑CGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCC > NC_000913/357745‑358031 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |