| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 993,913 | 0 | G | C | 53.2% | 54.1 / 76.0 | 47 | G44G (GGC→GGG) | ssuB | aliphatic sulfonate ABC transporter ATPase |
| Reads supporting (aligned to +/- strand): ref base G (17/5); new base C (0/25); total (17/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.61e-09 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.57e-07 | |||||||||||
AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCT > NC_000913/993770‑994057 | aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGCGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑Cccacc > 1:685847/1‑150 (MQ=255) cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGcc > 2:626072/1‑150 (MQ=255) cgcgCATCCTGAAACATCATTCGCGTGTCTGCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGcc < 2:1001022/150‑1 (MQ=255) catTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCgg > 2:26816/1‑150 (MQ=255) ttCGCTTTTCTTCCTAATTTCCCGCCGGGGGTGTGGTGCCCGCTAACACCTCGCCGGCGTTTGGCTTTTCC‑AGCCCGGCCAGCGGGCGCGGCGGGTTTTTTTTCCCCCCCCCGCGGCGC‑CCCCCCACCGCCACAAACTGACCTGcc < 1:666772/146‑1 (MQ=255) tgcgctTGTCTTCCTAATTTTCAGCCGGCGGTGGGGTGCCCGAAACCACCCCGCCGGGGGTGGGCTTTCCC‑AGCCCTGCCAGCGGGCGCAGCGGGTTTTTTTTCCCCCCCCCGCGGCGC‑CCCCCCACCGCCACAAACTGACCTGCCg < 1:1044361/145‑1 (MQ=255) cgggtttctcCCTAAATTTCACCCAGCGGTGTGGTCCCCGCTAACACATCGCCGGGGGTTGGCTTTTCC‑AGCCCTGCCAGCAGGCGCAGCGGGTTCTTTTTCCCCCCCCCGCTGCGC‑CCCCCCCCCGCCACAAACTGACCTGCCGGAat < 1:977922/143‑1 (MQ=255) gtgtCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGt > 1:566476/1‑150 (MQ=255) ttccTTCCTGAGTTCCAGCGAGCGGTGTGGTGCCGGCTACCGCCCCGCGTGCGGTTGGCTTTCCC‑ACCCCAGCCAGCAGGCGCAGCGGGGTTCTTTGGCCCCCACCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTa < 1:745707/147‑1 (MQ=255) tCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc > 2:657683/1‑150 (MQ=255) gAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg > 1:239183/1‑151 (MQ=255) ttccacccgggggTGTGGTCCCCGCAACCACACCCCCTGGGTTTGGCTTTCCC‑AGCCCGGCCAGCAGGCGCAGCGGGGTTCTTTTCCCCCCCCCGCGGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg < 1:587303/137‑1 (MQ=255) ttCACCCGGGGGTGGGGTCCCCGCTAACACCCCGCCTGGGGTTGGGTTTCCC‑AGCCCTGCCAGCGGGCGCAGCGGGTTTTTTTTGCCCCCACCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt < 1:46944/150‑1 (MQ=255) ccaggGGTGTGGTCCCCGCAAACCCACCGCCTGGGGTTGGCGTTCCC‑AGCCCTGCCAGCGGGCGCAGCAGGGTACTTTTCCCCCCCCCGCGGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGa < 1:31263/150‑1 (MQ=255) ggggTGTGGTCCCCGCAAACCCACCGCCTGCGGTTGGCTTTTCC‑AGCCCCGCCAGCAGGCGCAGCGGGTTTCTTTTGCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg < 1:665156/144‑1 (MQ=255) ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat > 2:487145/1‑151 (MQ=255) ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat > 2:674478/1‑151 (MQ=255) cccGCAAACGCACCGCCTGGGGTGGGCGTTCCC‑ACCCCCGCCAGCGGGCGCAGGAGGTTCTTTTTGCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgc < 2:373663/149‑1 (MQ=255) ccGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt > 2:770741/1‑150 (MQ=255) aaCACCCCGCCGGGGGTTGGCTTTCCC‑ACCCCTGCCAGCGGGCGCAGCGGGGTTTTTTTCCCCCCCCCCCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg < 1:398239/149‑1 (MQ=255) acaTCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc > 1:753668/1‑151 (MQ=255) cccgggttgggCTTTCCC‑AGCCCCCCCAGCAGGCGCGGCAGTTTTTTTTCCCCCCCACCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAAt < 2:613608/130‑1 (MQ=255) ccTGGGGTGGGCTTTCCC‑AGCCCTCCCGGCAGGCGCAGCAGGTTTTTTTTGCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt < 1:1035566/147‑1 (MQ=255) tGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc > 1:1042088/1‑150 (MQ=255) ggCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc > 2:122686/1‑150 (MQ=255) gCGTTTCCAAAACCTGCCGGCAGGCGCAGAGGGGTCCTTTTCCCCCCCCCGGTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc < 1:822985/150‑1 (MQ=255) ttccc‑AGCCCTGCCGGCGGGCGCGGAGGGGTTTTTTTCCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac < 2:1044382/146‑1 (MQ=255) ttccc‑AGCCCCGCCAGCGGGCGCAGCGGGTTTTTTTCCCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAg < 2:1084108/142‑1 (MQ=255) ttcccagcccgcccAGCAGGCGCAGCAGGTTCTTTTTCCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaa < 1:506732/136‑1 (MQ=255) cc‑AGCCCTGCCGGCAGGCGCAGCAGGTTACTTTTCCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt < 1:176426/149‑1 (MQ=255) agcccTGCCAGCGGGCGCAGCGGGTTTTTTTCCCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc < 2:701341/147‑1 (MQ=255) aCCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg > 1:845576/1‑151 (MQ=255) cccgCCAGCAGGCGCAGCAGGGTACTTTTCCCCCCCCCGCTGCGG‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg < 1:352485/147‑1 (MQ=255) ccccccAGCGGGCGCAGCGGGTTCTTTTTCCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg < 2:71551/144‑1 (MQ=255) gcagGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGtt > 1:368001/1‑151 (MQ=255) ggggcagggggttttttttCCCCCCCCCGCTGGGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa < 1:678960/135‑1 (MQ=255) cgcAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc > 2:542270/1‑151 (MQ=255) actggaggTACTTTGTCCCTCCCCGCCGCGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa < 2:1116178/143‑1 (MQ=255) ttatttttCCCCCCCCCGCTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTa < 1:1064645/146‑1 (MQ=255) ttcTTTTGCCCCCCCCGCTGCGC‑CCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTa < 1:2900/147‑1 (MQ=255) tttcGCCCCCCCCGCTGCGG‑CCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt < 1:41770/147‑1 (MQ=255) ttccccccccccGCGGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt < 2:228714/140‑1 (MQ=255) ccccccccGCTGCGG‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt < 2:678830/145‑1 (MQ=255) caccacCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg > 2:723479/1‑151 (MQ=255) ccaccGCCGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc < 2:669202/146‑1 (MQ=255) gCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct > 1:634938/1‑151 (MQ=255) cTGCGC‑CCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcc < 1:312617/149‑1 (MQ=255) | AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCC‑AGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGG‑CCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCT > NC_000913/993770‑994057 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |