Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 357,890 | 0 | C | A | 54.3% | 19.0 / 47.2 | 35 | L267F (TTG→TTT) | cynR | transcriptional activator of cyn operon; autorepressor |
Reads supporting (aligned to +/- strand): ref base C (10/6); new base A (0/19); total (10/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.36e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-05 |
AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC > NC_000913/357745‑358033 | aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc > 1:497992/1‑151 (MQ=255) aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc < 1:451938/151‑1 (MQ=255) cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc < 2:594183/151‑1 (MQ=255) cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc > 1:536046/1‑151 (MQ=255) gTAGCGCATCAGGCAATCCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc < 2:142299/151‑1 (MQ=255) cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc > 1:130580/1‑151 (MQ=255) gCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggg > 2:472613/1‑151 (MQ=255) cctcagggttgtgCCAGCGGCTACCGTTATTCATTCCCGCCAACACCCCCTCATTTATCCACCCCCATGACAAAAAATGCCTTCGGGGGGGCTCTCGGCCGTCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg < 1:372606/142‑1 (MQ=255) ccGGAGACTTGTGTCGGCGCCTCCCGTGATTCTTTCCCGCCAACAACCGCGCTTTTACCCACCGCCATGGGCAAAAATGCCTCCGGGGGGGCTGTCTCCCGGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg < 2:727116/149‑1 (MQ=255) ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg > 1:238996/1‑151 (MQ=255) caccacccctttttCAGCGCATCCCGTTATTATTTCCCCCCAACACCCGCCCTTTTCCCCACCCCCTTGTGAAAAAACCCCCCCGCGGCGTCTCTCGGCCATTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgg < 2:54896/136‑1 (MQ=255) tgtCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTATTTTATCCCCAACAAGACCGCCGCTCCCTCCACTGGcggcggagggcac > 2:372544/1‑138 (MQ=11) gcgcCTACCGTGTTTCTTTCCCCCAAACACCCGCGCTTTTACCAAACCCCTTGTGAAAAAATCCCTCCGGGGCGGCTGTCGCCCAGCTTTTTTTCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga‑g < 2:872968/141‑1 (MQ=255) ggCTCCCGTGATTCATTCCCGCCAACACCCGCGCATTTACCCAACGCCATGTGAAAAAATCCCTCCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTa < 1:668988/150‑1 (MQ=25) gcctcccGTGATTCATTCCCCCCAACACCCGCGCATTTACCCACCGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTa < 1:269100/145‑1 (MQ=255) aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGAGCAAGA‑GAAATAGCTTTAAAccc > 2:758188/1‑151 (MQ=255) ccGGGATTCTTTCCCCCCACCACCCGCGCATTTACCCACCGCCATGTCAAAAAATCCCTCCGGGGCGGCTGTCTCCCGGCTATTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCg < 2:719175/151‑1 (MQ=255) ttccccccAACACCCACGATTTTACCAACCGCCTTGTACAAAAATGCCTGCGGGGGGTCTTCCGCCCAGCTTTTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCAtgttg < 1:364364/143‑1 (MQ=255) ttcccGCAAACACCCGCGCTTTTACCCACCGCCTTGTGAAAAAAGGCTTTGGCGGCGGCTGTCTGCAAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCAtgtt < 1:140/145‑1 (MQ=255) ccccgtttttACCCAACCCAATGTAAAAAAGGCCCTGCGGGGCGGTTGTCGCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAAt < 2:821007/138‑1 (MQ=255) tttcccaaaCGCCATGTGAAAAAATCCCTTCGCGGGGCTTTTCTGCCAGTTTTTTTCCCCCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa < 2:445559/143‑1 (MQ=255) ccAACCCCAGGTGAAAAAATGCCTTCGGGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg < 1:838357/150‑1 (MQ=255) aCGCAATGTGAAAAATTGCCTCCGGGGCGGCTGTCTGCGAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGa < 1:372723/151‑1 (MQ=255) aaaaaTGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGGAAAAAATTTTAACCCCTCCAT‑TTTTGTGGGCATCGCTGCCGGGGACAATTTGGGAAGGGGAgggg > 2:52731/1‑148 (MQ=255) aacacgCCTGCGGGGGGGGTGTCTGCCTCTTTTTTGCCCTCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc < 2:218481/145‑1 (MQ=255) aaaaTGCCTTCGGGGGGCCTGTCTGCCAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg < 2:527261/151‑1 (MQ=255) gcttGTCTCCCACTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCa < 1:891660/146‑1 (MQ=255) tcccattttttttcCCCCCGCAAAAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg < 2:687437/139‑1 (MQ=255) ccAGCTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGa < 1:949809/142‑1 (MQ=255) ttttttcCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa < 2:407952/145‑1 (MQ=255) tttttcCCGCCGAAACAAAGCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc < 2:605589/140‑1 (MQ=255) tttttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GAAAAAAATTTTAAACCCCCCTGTTTTGTTGGGATCGCTGCTCGGAAAAAAGAGGGAAAGGGAGGGGGGCGAGTCAGCCCCCCAACCGCCCCAATTGAg > 1:100075/1‑150 (MQ=255) ttCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGA‑GAAAAAAATTTTAAACCCCCCTGTTTTGTGGGGAACGCTGCTGGGGAAAAAATGGGAAGGGGAGAGGGGGGGATCAACCCCCGAACCGCGCTAAttttgttt > 1:293637/1‑144 (MQ=255) cccgccgAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTc < 1:766048/148‑1 (MQ=255) ccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTc < 1:745670/150‑1 (MQ=255) | AAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTC > NC_000913/357745‑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: ‑ |