| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 993,913 | 0 | G | C | 44.4% | 88.0 / 25.6 | 54 | G44G (GGC→GGG) | ssuB | aliphatic sulfonate ABC transporter ATPase |
| Reads supporting (aligned to +/- strand): ref base G (27/3); new base C (0/24); total (27/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.17e-12 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCT > NC_000913/993763‑994057 | ccACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGAGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCgg > 2:1149075/1‑151 (MQ=255) aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCccacca < 2:1106329/151‑1 (MQ=255) gcgcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCcac > 1:1035808/1‑151 (MQ=255) aaaCATCATCCGGGTGTCTACCTAAATTTCCCCCGGCGTTGT‑GTGCCCCGCAACCCCACCCCCGGGGGTGGCCTTC‑ACACGCCCTCCCACCGGGCCCAGGAGGGTAC‑TTTTC‑CCACCC‑CGGCTGGGCCCCACCACCGCCACAAACTGa < 1:214408/148‑1 (MQ=255) atcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCg > 2:1218896/1‑151 (MQ=255) catGCGCGTGCCTTCTATATTTTCACCCGGGGGTGTGGGCCCCGTAAACTCACCCCCGGGGGTGGGTTTTC‑CCA‑CCCCGGCCAGCAGGCGGGGAGGGGTAC‑TTTCG‑CCACCA‑CCGCCGCGCCCCCCCACCGCCACAAACTGACCTGCCGGa < 2:790926/151‑1 (MQ=255) tCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAat > 1:553688/1‑150 (MQ=255) gcgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATg > 2:1037599/1‑151 (MQ=255) tgtCTTCCTAAATTCCAGCCGGGGGTGGGGTCCCCGCTACCCCACCGCCTGGGGTTGGCTTTC‑CCA‑GCCCGCCCAGCAGCCGCAGCGGGGTTT‑TTTGC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGaa < 2:347527/144‑1 (MQ=255) ttttCCTGAATTTCAGCCGGGGGTGGGGTCCCCGCAAACCCACCGCCTGGGGTGGGCTTTC‑CCA‑GCCCCCCCGGCGGGCGCAGCAGGGTTT‑TTTTC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAat < 2:76845/141‑1 (MQ=255) cTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATc > 2:1090913/1‑151 (MQ=255) ttCCTAATTTCCACCCGGGGGGGTGGCCCCCGCTAACCCCCCCCCTGGGGTGGGCTTTC‑CCA‑GCCCCGCCAGCGGGCGCAGCAGGTTTC‑TTTCC‑CCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAatat < 1:1114408/143‑1 (MQ=255) cctgatttccccccaggGGTGGGGCCCCCTCAAACACACCGCCGGGGGTTGGCGTTC‑CCA‑GCCCCCCCGGAGGGCGCGGCGGGTTAC‑TTTCC‑CCACCC‑CCGCTGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAg < 1:254826/139‑1 (MQ=255) gAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg > 2:675246/1‑151 (MQ=255) attttCAGCCGGGGGTGTGGTCCCCGCAAACACTCCGCCGGGGTTGGGCTTTC‑CCA‑GCCCGGCCAGCAGGCGCAGCGGGGTTT‑TTTTG‑CCCCCC‑CCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt < 2:58845/147‑1 (MQ=255) aaTTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg > 1:1249844/1‑151 (MQ=255) tttCACCCGGGGGTGTGGTGCCCGCAAACACATCGCCTGGGGTTGGTTTTC‑CCA‑GCCCGGCCAGCGGGCCCAGCGGGTTTT‑TTTCG‑CCCCCA‑CCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt < 2:496481/150‑1 (MQ=255) cagccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc > 2:286474/1‑149 (MQ=255) cagccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc > 2:970385/1‑149 (MQ=255) agcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg > 1:18025/1‑150 (MQ=255) cGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat > 1:840224/1‑150 (MQ=255) gtggtgCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt > 2:482111/1‑151 (MQ=255) gggggccccGCAACCCCACCGCCTGCGGTGGGCTTTC‑CCA‑CCCCTCCCAGCGGGCCCAGCAGGTTTT‑TTTCC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt < 1:119830/143‑1 (MQ=255) ggtgCCCGCTAACACACCGCCTGCGGTGGGCTTTC‑CCA‑GCCCTGCCAGCAGGCGCGGCAGGTTTT‑TTTTC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt < 2:546661/148‑1 (MQ=255) cccGCTACCACATCCCCGGGGGTTGGCTTTC‑CCA‑GCCCCCCCAGCAGGCGCAGCGGGTTTT‑TTTTC‑CCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt < 2:146794/151‑1 (MQ=255) cTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg > 1:658677/1‑151 (MQ=255) accacaCCGCCGGGGGTGGGCTTTC‑CCA‑GCCCTCCCAGCGGGCGCAGAGGGTTTT‑TTTTC‑CCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAAt < 2:1055750/146‑1 (MQ=255) aaCACACCGCCGGGGGTTGGTTTTC‑CAA‑GCCCGCCCAGCGGGCGCAGCGGGGTAT‑TTTTC‑CCCCCC‑CCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt < 1:915289/150‑1 (MQ=255) gggggtGGCGTTTCCCA‑GCCCTCCCAGCGGGCGCA‑CGGGGTTC‑TTTTGCCCCCCC‑CCGCGGCGCCCCCCGACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACt < 1:1129579/143‑1 (MQ=255) gCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt > 2:481346/1‑151 (MQ=255) ggTTGGCCTTTCCCA‑GCCCTGCCGGCAGGCGCGTCGGGCTTC‑TTTGC‑CCACCC‑CCGGTGCCCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt < 1:451683/143‑1 (MQ=255) ggTGGGCTTTC‑CCA‑GCCCTGCCAGCGGGCGCAGAGGGGTAC‑TTTTC‑CCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt < 2:622488/138‑1 (MQ=255) ttgaggttcccaccccTCCCAGCGGCCGCAGCGGGTTTT‑TTTTC‑CCCCCG‑CGGCTGCGCCCCCCTGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc < 1:7732/131‑1 (MQ=255) tgggctttccaagcccTGCCAGCAGGCGCAGCGGGTTTT‑TTTTC‑CCCCCC‑CCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa < 2:129246/136‑1 (MQ=255) ggctttcccaccccTGCCAGCAGGCGCAGCGGGTTTC‑TTTGC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAg < 2:530056/137‑1 (MQ=255) ttt‑CCA‑GCCCGGCCGGCGGGCGCAG‑AGGGTTCTTTTTG‑CCCCCC‑CCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac < 2:612424/149‑1 (MQ=255) ccggCCGGCGGGCGCAGCGGGGTTT‑TTTTC‑CCCCCC‑CCGCCGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATg < 2:95898/142‑1 (MQ=255) catcctttccggCGCGGAAGGCTAT‑TTTTC‑CCACCC‑CTGTTTGCGCCCCCCGCCGCCGCAAACGGGCCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAtt < 2:101394/124‑1 (MQ=255) cTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGcc > 1:1186805/1‑151 (MQ=255) tGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGccc > 1:960066/1‑151 (MQ=255) ccAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTg > 2:689655/1‑151 (MQ=255) cagGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc > 1:1237611/1‑151 (MQ=255) ggcgcagggggtttttttG‑CCCCCCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt < 2:519050/133‑1 (MQ=255) ggttcttgag‑CCACCC‑CGGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAg < 2:792697/137‑1 (MQ=255) aC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc > 2:375034/1‑151 (MQ=255) tttG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt > 1:1019560/1‑150 (MQ=255) ccaccc‑CCGGTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg < 2:598795/150‑1 (MQ=255) ccaccc‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATa < 1:982685/147‑1 (MQ=255) tacca‑ccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg > 2:1115957/2‑151 (MQ=255) cacca‑cCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg > 2:574914/1‑151 (MQ=255) acca‑ccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg > 2:610264/1‑151 (MQ=255) acca‑ccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg > 1:18567/1‑151 (MQ=255) acca‑ccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg > 1:750028/1‑151 (MQ=255) ca‑ccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc > 2:820476/1‑151 (MQ=255) gCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGGGGcctcct > 1:1122769/1‑151 (MQ=255) | CCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCA‑GACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCT > NC_000913/993763‑994057 |
| 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: ‑ |