Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 61.8% | ‑6.5 / 33.6 | 34 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (7/6); new base A (0/21); total (7/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.19e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-03 |
CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATT > NC_000913/357740‑358022 | cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaac > 2:276571/1‑151 (MQ=255) aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGcaacaa < 1:611284/150‑1 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAAcc < 2:514164/151‑1 (MQ=255) gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtct > 2:493693/1‑151 (MQ=255) cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTtct > 1:414649/1‑150 (MQ=255) atttCCGGCCGCAGCCTTGTTTCAGCGGCTCCCGTGTTTCTTTCCCCCAAACACCCGCGCTTTTA‑CCCACCCCATGGGGAAAAAATCCTTT‑C‑GGGGCGCCTTTCTCCCAGTTTTTTTC‑CCGCCGAAACAAAACCGCCGTTCTCTCCagtag < 2:598886/149‑1 (MQ=255) attccccgCCGCAGCCTTGTGACAGCGGCACCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTCCCCACCGCCAT‑GTAAAAAAATGCCTC‑C‑GGGGCGCTGGTCGGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCagtag < 1:776056/144‑1 (MQ=255) aGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgg > 2:106950/1‑151 (MQ=255) ccGCAGACCGGTGACAGCGTCTCCCGGGTTTCATTCCCGCCAACAACCCCGCATTTACCAACCGCCTT‑GAGCAAAATCGCCTT‑C‑GGGGCGGCTGCCGGCGAGTTTTTTCC‑CCCCCGAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCggg < 2:198595/150‑1 (MQ=255) aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 1:142892/1‑147 (MQ=255) cccggtttCAGCGGCCACCGTGTTTCATTCCCCCCAACACCCGCCCTTTTCCCCAACGCCAT‑GGGAAAAAATCCCTT‑C‑GGGGGGGCTGTCTCCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga < 1:180955/143‑1 (MQ=255) cttgtgTCAGCGTGTCCGGTTATTCATTCCCACAAAAACCCGCTCTTTTCCCCACCGCCTT‑GT‑CAAAAAGCCTTTCC‑GCGGCGTCTTTC‑GCCACTTTTTTTTCCCGCGGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAagag < 2:836393/148‑1 (MQ=255) tgtgTCGGCGCCTCCCGTGATTCTTTCCCCCCAACACCCGCGCTTTTCCCCACCGCCAG‑TTGCAAAAGTGCCTG‑C‑GGGGGGCCTTTCGGCCAGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg < 1:161348/139‑1 (MQ=255) tcgggtgtccccttGATTCTTTCCCCCCAACACCCGCCCATTTTCCCAACGCCTG‑GTGAAAAAACCCTTG‑C‑GCGGCGGTTGTCTGCCGGTTTTTTTC‑CCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg < 2:273589/137‑1 (MQ=255) cAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGaaaaaa > 1:716750/1‑146 (MQ=255) aGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAagagggaaaaaat > 1:337391/1‑142 (MQ=255) tccggTGATTCATTCCCGCAAACACCCGCGTTTTTACCCAACCCCAT‑GTGAAAAAAGGCTTT‑C‑GCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt < 2:940887/142‑1 (MQ=255) gtgtTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCAAT‑GTGAAAAAATGCCTC‑C‑GGGGCGGCTGTCGGCCAGTTTTTTTC‑CCGCCGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTc < 1:180233/151‑1 (MQ=255) ttattcattcCCCCAAACACCCGCTCTTTTACCAACCGCCTT‑GAAAAAAAATGCCTT‑G‑GCGGCGCCTGTCGGCCAGTTTTTTTC‑CCGCCACAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGt < 1:101915/147‑1 (MQ=255) caccccctcATTTACCCACCGCCTG‑GTGAAAAAATGCTTC‑C‑GCGGGGGCTGTCGGCCGGTTTTTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt < 2:933172/140‑1 (MQ=255) ccGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGGAAAAAATTTTAACCCCTCCTTGTGTGTggggagcgccg > 2:85402/1‑141 (MQ=25) cgcgCTTTTACCCACCCCCAT‑GTGAAAAAGTGCCTC‑C‑GCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct < 2:953333/151‑1 (MQ=255) gtttttACCCACCGCAAT‑GTAAAAAAAT‑CCTT‑CGGGGGGGGTTGTCTGCCAGCTTTTTTC‑CCCCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctg < 2:557435/143‑1 (MQ=255) tttACCCACCCCCAT‑GTGAAAAAATCCCTC‑C‑GCGGCGGTTGTCTCCCAGCTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa < 2:796034/150‑1 (MQ=255) aaCGCCAT‑GTAAAAAAATGCCTC‑C‑GCGGCGGTTTCCTGCCAGTTTTTTTC‑CCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgtg < 2:99660/150‑1 (MQ=255) ccAT‑GTGAAAAAATGCCTC‑C‑GCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaa < 1:286114/149‑1 (MQ=255) caagggaaaaaaaTGCTTC‑C‑GGGGCGCCTGCCTGCAGTCTTTTTCG‑CCGCCAAACAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaa < 1:39020/141‑1 (MQ=255) gtgAAAAAAGCCCTT‑C‑GCGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg < 2:165168/148‑1 (MQ=255) aaaaaTGCTTC‑C‑GGGGGGTCTGTCTGCCAGTTTTTTTC‑CCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt < 2:167136/149‑1 (MQ=255) aaaaaTGCCTC‑G‑GGGGCGGCTGTCTGCCAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc < 2:519086/151‑1 (MQ=255) aaaaaTGCCTC‑C‑GGGGGGGCTGTCGGCCAGTTTTTTTC‑CCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt < 2:958628/149‑1 (MQ=255) tttCTGCCAGCTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa < 2:705637/147‑1 (MQ=255) ctcccAGCTTTTTTC‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 2:304902/145‑1 (MQ=255) cTGCCATTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa < 1:393878/146‑1 (MQ=255) ccAGTTTTTTTC‑CCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 1:140/145‑1 (MQ=255) ttttc‑CCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt < 1:562473/147‑1 (MQ=255) | CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCAT‑GTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTT‑CCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATT > NC_000913/357740‑358022 |
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: ‑ |