| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,890 | 0 | C | A | 58.8% | 3.3 / 45.1 | 34 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base C (13/1); new base A (0/20); total (13/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.51e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-06 | |||||||||||
ATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC > NC_000913/357743‑358033 | aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa > 2:971084/1‑150 (MQ=255)aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa > 2:1181281/1‑150 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc < 2:1242629/151‑1 (MQ=255) cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc > 2:1127108/1‑151 (MQ=255) gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCCCCGCAACAAAAccgcc > 2:158820/1‑151 (MQ=255) cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct > 1:929417/1‑150 (MQ=255) cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt > 1:111249/1‑150 (MQ=255) ccccccGCAGCCTTGTGTCGGCGGCTCCGGTTTTTCATTCCCCCCAACACCCCCCATTTTACCCAACCCCGTGTAAAAAATTCCCTG‑C‑GCGGCGGCTTTCTCCCAGTTTTTTTCCCCCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGc < 2:470400/147‑1 (MQ=255) cggCAGCCTTGTGTCGGCGCCTACCGTGTTTCTTTCCCCCAAACACCCGCGATTTTACCCACCGCCATGTGAAAAAATGCCTT‑C‑GGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg < 2:985710/149‑1 (MQ=255) ccgaagccttgtgTCGGCGCCTCCCGTGTTTCTTTCCCCCCAACACCCGCCATTTTACCCACCGCATGGTACAAAAATGCTTT‑C‑GGGGCGGCTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg < 2:743190/140‑1 (MQ=255) cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg > 1:626346/1‑147 (MQ=255) gACCTGTGTCGGCGGCCCCGGTTATTCTTTCCCCCAAACACCCGCGCTTTTCCCCACCGCTTTGT‑AAAAAATGGCTTCC‑GGGGCGTTTTCCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTg < 1:161335/140‑1 (MQ=255) ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggtgg > 2:1069296/1‑142 (MQ=255) cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggg > 1:53881/1‑143 (MQ=255) cggcgcCTCCCGTTATTCTTCCCCCCAAACACCCGCGCTTTTACCAAACGCCATGTGCAAAATTCCCTT‑C‑GCGGCGGCTTTCTGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGc < 2:779921/148‑1 (MQ=255) aGCGCCTCCGGTGTTTCTTTCCCCCCAACACCCGCCCTTTCACCCAACCCCATTTAAAAAAATCCCTC‑C‑GGGGCGGCTGTCTCCCGTTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga < 2:800615/143‑1 (MQ=255) cGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaaaattt > 2:1242825/1‑143 (MQ=255) ggCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggacgggaaaagtc > 2:591182/1‑133 (MQ=255) cTACCGTGTTTCTTTCCCCCCACCACCCGCGATTTTACCCAACACCTTGTGAAAAATTCCCTC‑G‑GGGGCGGCTGCTTGCAAGTTTTTTTTCCCCCAAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg < 2:176929/150‑1 (MQ=255) cGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGaaaaaaattttaaacc > 2:277174/1‑136 (MQ=25) ttccccccAACACCCGCGCTTTTACCCACCCCCTTGTGAAAAAATCCCTC‑C‑GGGGCGGTTTTCCGCCATTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttgt < 2:168462/144‑1 (MQ=255) cccacaaacaccCGCGCTTCTACCAAACCCCAGGTAAAAAAAGCCCTC‑C‑GCGGCGGCTGTCTCCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgtt < 2:141035/140‑1 (MQ=255) cacccccgcttttccccACCGCCTTGT‑CAAAAATGCCCT‑CGGGGGCGGTGGTCGGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCg < 1:987038/135‑1 (MQ=255) acCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAACCCCTCCTGTTTTGTggggaactc > 2:197298/1‑142 (MQ=25) ccGCGCTTTTCCCCAACCCTTTGTAAAAAAATCCTTG‑C‑GGGGCGCCTGTCCGCAGGTTTTTTCCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc < 2:50995/146‑1 (MQ=255) tttACCCACCCCCATGGGAAAAAATGCTTT‑C‑GCGGCGGCTGTCTGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt < 2:646651/148‑1 (MQ=255) aaCGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCGGGGAAAAAAATTTTAAACCCCCCTGTTGTGTGGCGATCgctgcttgggacaaaat > 1:142527/1‑138 (MQ=21) aaaaaaTGCCTT‑G‑GGGGCGGTTGTCTCCGGGTTTTTTCCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt < 1:685720/149‑1 (MQ=255) aaaaaaTCCCTC‑C‑GGGGCGGCTGTCTCCCATTTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg < 2:553014/150‑1 (MQ=255) aaaaaTGCTTG‑C‑GGGGGGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa < 2:1061174/146‑1 (MQ=255) aaaaaTCCCGC‑C‑GCGGGGGCTGGCTCCCACTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa < 1:1189248/147‑1 (MQ=255) gcggggGCTGTCTGCCGGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc < 1:541482/147‑1 (MQ=255) caacattttttcCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc < 2:430608/137‑1 (MQ=255) cagtttttttcCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg < 1:1178385/139‑1 (MQ=255) ccgccaAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt < 1:1101045/146‑1 (MQ=255) cgcAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAAAAAATTTAAACCCCCCCTCTTTTGGGGGGATCGCTTCTGGGGACAAAGTGGGAAGGGGAGGGGGGGGGGTCAGCCCCCGAAACGCCCCAATTTAGTTTGcccc > 1:387227/1‑147 (MQ=255) | ATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTT‑C‑GCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC > NC_000913/357743‑358033 |
| 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: ‑ |