| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 357,887 | 0 | C | A | 42.9% | 27.1 / 10.4 | 28 | L268F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
| Reads supporting (aligned to +/- strand): ref base C (11/5); new base A (1/11); total (12/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.05e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACC‑GTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC > NC_000913/357738‑358033 | gCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACC‑GTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc > 2:488710/1‑150 (MQ=255) gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACC‑GTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca < 2:1123027/150‑1 (MQ=255) gATAAGCGTAGCGAATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACC‑GTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa < 1:1233611/151‑1 (MQ=255) gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACC‑GTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtct > 1:1301032/1‑151 (MQ=255) caacagcaaaTCCCACCCGAAGCCCTGTGCCAGCGGCTCCG‑GTGTTTCTTTCCCGCCAACAA‑CCGCCCTTTTACCCAACCCCATGTCAAAAAATGCCTTCGCGGCGGCTGTCTCCCAGCTTTTTTTCCGCCGCAAAAAAAccgccg < 2:1449744/138‑1 (MQ=255) atcggtcattcccACCCCCAGACCTGTGTCGGCGGCTCCC‑GTGATTCATTCCCGCAAACAC‑CCGCGCATTTACCCAACGCAATGTGAAAAAATGCTTTCGCGGCGGCTGCCTGCCGGTTTTTTTCCCGCCGCAACAAAAccgccg < 2:1329850/134‑1 (MQ=255) aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACC‑GTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCCGCTATTTTTCCGCCGAAACAAAACCGCCGTtctc > 2:530462/1‑150 (MQ=255) ccGCACTCCTGTTGCAGCGG‑TACCGGTGTTTCTTTCCCACCACCAC‑CCGCGCATTTACCAACCCCCTTGTCCAAAAATCCCTGCGGGGCGTCTGTCACCCGATTATTTTCCCGCCAAAAAAAATCCGCCGTTCTCTCCAGTAGTg < 1:1457655/145‑1 (MQ=255) gCAGACCTGTGTCAGCGGCTACC‑GTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 1:522872/1‑149 (MQ=255) gACCTGTGTCAGCGGCTACC‑GTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 1:1172263/1‑146 (MQ=255) aCCTGTGTCAGCGGCTACC‑GTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggg > 2:91509/1‑145 (MQ=255) tCAGCGGCTACC‑GTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggaaggaaaaa > 1:1218774/1‑138 (MQ=255) cAGCGGCCCCC‑GTATTTCTTTCCCGCCAACAC‑CCGCGCATTTACCCAACCCCATGTGAAAAAATCCCTCCGGGGCGGCTGTCTGCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaa < 1:22386/145‑1 (MQ=255) ttCCCCCAACAACCCCCGATTTTATCCACCCCCTTTTCAAAAAATCCCTCGGGGGGGGCTGTCTCCCATTTATTTTCCCCCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg < 2:313059/149‑1 (MQ=255) cccccAACAC‑CCGCGCTTTTCCCCAACGCCATGTAAAAAAATGCTTGGGGGGCGGCTGTCTGCCGTTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgtt < 2:566333/142‑1 (MQ=255) acaa‑cCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAATTTAAAACCGCCAtactttgggggact > 2:1073182/1‑136 (MQ=25) tCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAATTTTAACCCCTCCTGTTTTGTGGGGATCCCTGCGGGTGACa > 1:806138/1‑148 (MQ=25) ccaaccgcctggtgAAAAAATCCCTGGGGGGGGGCTGTCGCCCAGTTTTTTCCCCGCGGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACa < 2:688208/137‑1 (MQ=255) aaCGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTTAAACCCCCATGTTTTTTGGCAAACGCCGCTGGGAaaaa > 1:490211/1‑146 (MQ=25) cccccttgccaaaaaaTCCTTCCCCGGCGGCTGTCTCCCATTTTTTTTCCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg < 2:830527/139‑1 (MQ=255) cgcattgtaaaaaatTGCCTTCGCGGGGGCTGTCGCCCAGTTTTTTTTCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa < 1:1425396/140‑1 (MQ=255) cATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAACCCCGCCTGTTTTGTGGGGATTGCTGCGGGGAaaaaaaggggaa > 1:79831/1‑141 (MQ=19) ggggcggtTGTCGGCAAGTTTTTTTCCCCCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt < 2:116566/144‑1 (MQ=255) ggttGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg < 1:1450322/147‑1 (MQ=255) tgtcgccaagtttttttcCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGa < 1:91454/123‑1 (MQ=255) tcccATCTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 2:1300927/145‑1 (MQ=255) gcaagtttttttcCCCCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 2:360993/138‑1 (MQ=255) ccAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 1:374837/145‑1 (MQ=255) ccAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAAATTAAATCCGCCCCGTTGTGTGGGGAATTCCGCCGGGGAAAAAGTGGGAAGGGAAGAGGGGGGGAACAACTCCAGaaacc > 2:819000/1‑143 (MQ=255) ccgtttttttcCCCCCACAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc < 1:965859/138‑1 (MQ=255) cTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACTTTTAACCCCTCCTTTTGTGTGGGAATCCCTGCTGGGTAAAAAGTGGGAAAGGAAGAGGGGGGGATCAGACCCACAAACgcgcca > 1:419721/1‑145 (MQ=255) ccgaAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTc < 2:932063/150‑1 (MQ=255) | GCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACC‑GTGATTCATTTCCGCCAACAA‑CCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC > NC_000913/357738‑358033 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |