Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 993,904 | 0 | A | C | 58.3% | 24.4 / 30.5 | 25 | G47G (GGT→GGG) | ssuB | aliphatic sulfonate ABC transporter ATPase |
Reads supporting (aligned to +/- strand): ref base A (6/4); new base C (0/14); total (6/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.56e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-05 |
CCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAC > NC_000913/993763‑994046 | ccACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCgg < 2:585378/151‑1 (MQ=255) aCGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGGCCTG‑CCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGc < 1:827208/150‑1 (MQ=255) aCGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGCCCTG‑CCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGc < 1:325272/150‑1 (MQ=255) aaaccatcttTCGCGTGTCTTCCTAAATTCCACCCAGCGGTGTGGGGCCCCCAAACACACCGCCGGCGGTTGGCTTTTCCAGCCCTC‑CCAGCGGGCGCAGCAGGTTACTTTTGCCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGAcc < 2:339667/147‑1 (MQ=255) catTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGa > 2:659732/1‑151 (MQ=255) tCGCGTTTCTTCCTGAATTTCACCCAGGGGTGTGGTCCCCGCTAACACATCGCCTGCGTTTGGCGTTCCCAGCCCTG‑CCAGCGGGCGCAGCGGGGTACTTTTGCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAata < 2:140219/151‑1 (MQ=255) gAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg > 1:823967/1‑151 (MQ=255) ccgggggggtggcccccGCAAACACACCGCCTGGGTTGGGCGTTCCCAGCCCCG‑CCAGCGGGCGCAGCAGGGTATTTTTGCCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc < 1:868479/137‑1 (MQ=255) gcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat > 1:521772/1‑151 (MQ=255) gtggtgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt > 1:120415/1‑151 (MQ=255) ccccgctggCGGTGGGCGTTTCCAGCCCTG‑CCAGCGGCCCCAAGGGGGTATTTTTGCCACCCCCGCTGCGCCCCACCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCt < 1:648392/143‑1 (MQ=255) gCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTa > 1:218383/1‑149 (MQ=255) cggCGTTTGGCTTTCCCAGCCCTC‑CCAGCGGGCGCAGCGGGTTTCTTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc < 1:386018/149‑1 (MQ=255) ggcggtgggggttCCAAACCCGCCCAGCCGGCCCAGAGGGTTTCTTTCGCCCCCCCCGCGGCGGCCCACAACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc < 1:100245/139‑1 (MQ=255) tGGCTTTTCCACCCCGC‑CCAGCAGGCGCAGCGGGGTTTTTTTCCCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc < 2:607596/151‑1 (MQ=255) ggCTTTCCCAGCCCTG‑CCAGCAGGCGCAGCGGGGTATTTTTCCCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc < 2:689999/150‑1 (MQ=255) ttcccagcccccccAGCAGGCGCAGCGGGGTCCTTGTGCCCCCACCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa < 1:73340/139‑1 (MQ=255) accccgcccAGCAGGCGCAGCGGGGTTCTTTGCCCCCCCCCGTGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg < 1:461012/144‑1 (MQ=255) gcgcAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa > 2:306910/1‑151 (MQ=255) cgcAGCGGGGTACTTTCCCCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc < 1:834723/151‑1 (MQ=255) cgggttttttttCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc < 2:22329/142‑1 (MQ=255) ggttACTTTGCCCCCCCCCCCTCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAg < 1:837422/147‑1 (MQ=255) ttccTTTGGCCCCCGCGGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt < 1:91931/148‑1 (MQ=255) ttGCCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc < 2:833716/151‑1 (MQ=255) ccgCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTCTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc < 1:249506/149‑1 (MQ=255) | CCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAC > NC_000913/993763‑994046 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |