| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,223,634 | 0 | T | G | 56.2% | 10.7 / 19.2 | 16 | intergenic (‑6/‑194) | iclR/metH | transcriptional repressor/homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent |
| Reads supporting (aligned to +/- strand): ref base T (1/6); new base G (9/0); total (10/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.74e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-03 | |||||||||||
ACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTT > NC_000913/4223498‑4223783 | acacTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGGTTCGCGGGGATGGGGGGGGCCATGACAGTCTC‑ttttttttttt > 2:182314/1‑146 (MQ=255)acacTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTGTTTTGCCGCGTTTCGCGGGGGTGGGGGGGACCATGACAGGCTC‑TTTTTTTTGtt > 1:79030/1‑150 (MQ=255) ggATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGCTTTGCGGGGAAGGGGGGGGGTCTGGCAAGCTCTTTTTTTTTGTATTGTGGaaa > 2:174303/1‑150 (MQ=255) cGGCAATCCCCTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGGTTCGCGGGGATGGGGGGGGCCCTGCCAGGCTC‑TTTTTTCTGTTTCTTGGAAAtcatttt > 2:73246/1‑151 (MQ=255) tttCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGATGGGGGGGGCCATGACAGGCTC‑TTTTTTTTGTTTCGTGGGAATCATTTTTATttttttttttttctta > 2:276654/1‑145 (MQ=255) gTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTCGCGGGGATGGGGGCGGCCATGACAGGCTC‑TTTTTTTTGTATCGTGGGAATCATTTTTATTTTTTTTTTTTGCTAATGCAATAGttt > 1:144042/1‑148 (MQ=255) gTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCCCGGTTTGCGGGGATGGGGGCGGCCCTGACCGGCTC‑TTTTTTTTGTTTCGTGGGAATCCTTTTTCTTTTTTTTTTTTGCTTATTCAATAGTTa > 1:112669/1‑149 (MQ=255) aaGACTGAACCTGTCCAGTCGCGGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaac < 2:25558/151‑1 (MQ=255) cTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATg < 2:130241/151‑1 (MQ=255) gTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGGGTTTTCCGGGGAGGGGGGGGGCCCTGACAGGCTC‑TTTTTTTTTTTTTTTGGGAATCATTTTTTTTTTTTTTTTTTTCTTATTCAATAGTTTCTTAACTGATCCggttgg > 1:185978/1‑146 (MQ=255) tCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGGCGCGTTTCGCGGGAATGGGGGCGACCCTGGGCGGCTC‑TTTTTTTTGTATCGTGGGAAACATTTTTATTTTTTTTTTTTGCTAATGCAATTGGTATTGTACTGATTCGATggggttatt > 2:314374/1‑143 (MQ=255) cTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTGCGGGGATGGGGGCGGCCCTGACAGGCTC‑TTTTTTTTGTATCGTGGGAATCATTTTTATTTTTTTTTTTTGCTAATTCAATAGTTACTTAACTTATTCGATGAGTTTATGttt > 1:10943/1‑150 (MQ=255) tttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtggtg < 2:400360/151‑1 (MQ=255) ttCGCGGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCTTTTTTATTGTTAGCTAATGCAATAGTTTCTGAGCTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTAcc < 2:101259/151‑1 (MQ=255) cgcgGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACcac < 2:144327/151‑1 (MQ=255) aCCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGcc < 2:261263/150‑1 (MQ=255) tctc‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATTAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCCCAGATGCGTTTATGCCAGTATGGttt < 2:392293/150‑1 (MQ=255) | ACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTT > NC_000913/4223498‑4223783 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |