| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,604,344 | 0 | C | T | 52.6% | 2.5 / 21.2 | 19 | S62F (TCT→TTT) | bglJ | bgl operon transcriptional activator |
| Reads supporting (aligned to +/- strand): ref base C (4/5); new base T (10/0); total (14/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.08e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-02 | |||||||||||
AATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAC > NC_000913/4604206‑4604492 | aaTGCGTCATGAGTAGTATCGGCATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGACAGGAGTAATTTCAGCAGGCCATGTCGCGGATCTCGTTTGCGCCGGGCAATttttttttttctttcct > 1:14845/1‑146 (MQ=255) cGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGACATTTTTTTTTTTTCTGCCATga > 1:170093/1‑149 (MQ=255) gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCATGAGAGGGGTGGGCAGAGAggg > 2:174001/1‑151 (MQ=255) ttGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGGAGTGAGCGCCGAGAGGGAtt > 2:78292/1‑151 (MQ=255) cGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGGAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCg > 2:315277/1‑149 (MQ=255) gggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGAtt > 2:148823/1‑151 (MQ=255) cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAg < 1:100984/148‑1 (MQ=255) aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCg < 2:288996/151‑1 (MQ=255) aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGGGAGCGCCGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCggt > 2:277596/1‑149 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCATTTTTTTTTTTTTCCCCTGAGAAGAGGGCGCCGAGAGGGGGTTTCTTTCCCGACTGAAATGGCGGTTAAAGTTCCCCGGAACCGGCGgtttg > 2:371818/1‑148 (MQ=255) aGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGgatga < 2:107801/151‑1 (MQ=255) ggATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGGAGTGAGCGCCGAGGGGGGTTATTTTTGCTTACTGAACTGGCGGTTAAATTTCCGCGGACCCGGCCTTTTGTTATTGCGGGTGATGATAtt > 1:281029/1‑148 (MQ=255) aTGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGc < 1:289455/151‑1 (MQ=255) cgGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATTAGGAGTGAGCGCAGGGGGGGGTTATCTTTCCTTACTGAACTGGCGATTAAGGTTCCGCGGACCCGGCGGTTTGTTATTGCGGGTGATGATATTTAAACTCGGCTGGtt > 1:334360/1‑150 (MQ=255) aTCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTTTCTTTCCTTACTTAAATGGCGGTTAAGTTTTCGCGTACCCGGCGGTTAGTTTTTGCGGATGATGATATTTAAACTCGGCTGATTGGtt > 1:281422/1‑151 (MQ=255) cGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGGTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCAtt > 2:76082/1‑151 (MQ=255) gcggcgGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCg < 2:296223/150‑1 (MQ=255) gTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGACGCTCGGCTGATTGGTTCATTGTCGCCATCAc > 1:258454/1‑151 (MQ=255) gTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAc > 1:258430/1‑151 (MQ=255) tCATTTTTTCTTTTTCTGCCATGAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAc > 2:213893/1‑150 (MQ=255) tttttctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCt < 2:34275/150‑1 (MQ=255) tttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTgg < 2:79265/150‑1 (MQ=255) tctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAc < 1:78267/150‑1 (MQ=255) | AATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGAC > NC_000913/4604206‑4604492 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |