| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 993,913 | 0 | G | C | 63.0% | 9.0 / 28.1 | 27 | G44G (GGC→GGG) | ssuB | aliphatic sulfonate ABC transporter ATPase |
| Reads supporting (aligned to +/- strand): ref base G (6/4); new base C (0/17); total (6/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.09e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.16e-04 | |||||||||||
CCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCC > NC_000913/993763‑994053 | 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) ttcccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc < 2:218383/142‑1 (MQ=255) ttGCCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc < 2:833716/151‑1 (MQ=255) ccgCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTCTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc < 1:249506/149‑1 (MQ=255) cccccccccGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg < 2:17501/144‑1 (MQ=255) ccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc < 1:251665/149‑1 (MQ=255) | CCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTG‑CCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCC > NC_000913/993763‑994053 |
| 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: ‑ |