Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,111,747 | 0 | G | T | 35.9% | 68.6 / 25.4 | 39 | N156K (AAC→AAA) | yiiQ | DUF1454 domain‑containing protein YiiQ |
Reads supporting (aligned to +/- strand): ref base G (23/2); new base T (0/14); total (23/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.96e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG > NC_000913/4111621‑4111885 | aTTAAACCCTTCAAGCGATTCAGATAGCCCCAGTTAAATCGGTCCAACAGCGAAGGCCACCCCCTTTC‑CGCCGTTGCCCGCGACAACATACCGCAGTCCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTCCCccgccg < 2:1142641/140‑1 (MQ=255) cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc > 1:418559/1‑140 (MQ=255) cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc > 1:278107/1‑140 (MQ=255) cAAGCAATTCAGATACCGCCAGTTAAATCGTTTAAACAGCGAAGGTCCCCCCTTTTT‑CGCCGTTGCCCGCAACAAAATAACGCGGTGCCCTTTCTGTCTGGGTGAAATATTGTTTTTTTTCCCCCGCCGTTAGTAGCGAc < 2:73743/140‑1 (MQ=255) ttCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtt > 2:1157491/1‑140 (MQ=255) cAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGGAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGtttt > 2:587208/1‑140 (MQ=255) gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt > 1:84027/1‑140 (MQ=255) gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt > 2:85618/1‑140 (MQ=255) ggtaaaaCAGCGAAGGCCACCCCTTTTC‑CGCTGTTGCCCGCAACAACAAAACGCAGTGCCCTTTCTGTCTCGGTTTAATATTTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt < 2:1014083/135‑1 (MQ=255) ttCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACGTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa > 2:840517/1‑140 (MQ=255) ttCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa > 1:1134713/1‑140 (MQ=255) ttCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa > 1:750722/1‑140 (MQ=255) ttCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa > 1:650954/1‑140 (MQ=255) ttCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa > 1:721391/1‑140 (MQ=255) tacAACAGCGAAGGTCATCCCCTTTC‑CGCCGTTGCCCGCGACAACATAACGCAGTGCCCCTTCTGTCTCGGTGTAATACTTTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa < 2:1150560/138‑1 (MQ=255) taaaCAGCGAAGGTCACCCCTTTTC‑CGCCGTTGTCCGCAACAACATACCGCAGTGCCCCTTCTGTCTGGGTGTAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAt < 2:750722/138‑1 (MQ=255) tCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAt > 1:302750/1‑140 (MQ=255) tCAACAGCGAAGCTCAGCCCCTTCT‑CGCCGTTGTCCGCGACAACATAGCGCAGTGCACCTTCTGCGTCGGTGCAATACCGTTTGTTTTTCGCCGCTGATAGGCGCGGCTGCCGTTTTTTCTGGCTTTGTGTTTTGGCCAt > 2:1134509/1‑140 (MQ=255) ccAACGGCGAAGGTCCCCCCTTTTT‑CGCCGTTCTCCGCAACAACATAAGGCATTGCCCTTTCGTTCTGGTTTTAATAGTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAt < 2:1134713/139‑1 (MQ=255) caaaCAGGGAGGGCCACCCCTTTTC‑CGCCTTTGCCCGCGACAAAATAACGCAGTCCCCTTTCTGTCTGGGTTTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAt < 2:721391/138‑1 (MQ=255) aCAGCGAAGGCCACCCCTTTTT‑CGCCGTTGCCCGCGACAACATAACGCAGTCCCCCTTCTGTCTCGGTTTAATATCTTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTaa < 2:84027/140‑1 (MQ=255) gAAGGTCCCCCCTTTTC‑CGCCGTTGCCCGCGACAACATAACGCAGTGCCCCTTTTGTCTCGGTTTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATgggg < 2:834764/140‑1 (MQ=255) aaGGTCCCCCCTTTTT‑CGCCGTTGCCCGCGACAACATAACGCAGTGCCCCTTCTGTCGGGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt < 2:27707/140‑1 (MQ=255) aaGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt > 2:912488/1‑140 (MQ=255) ccaccccTTTTCCGCGGTTGCCCGCAAAAAAATACCGCAGTGCCCCTTCTGTCTGGGTATAATAGCGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAgt < 2:433519/137‑1 (MQ=255) ccctTTTC‑CGCCGTTGCCCGCGACAACATAACGCGGTCCCCCTTCTGTCTGGGTGTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGc < 2:650954/140‑1 (MQ=255) tttc‑CGCCGTTGCCCGCAAAAACATAACGCAGTCCCCCTTCTTTCTGGGTGTAATAATTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGAt < 2:418559/140‑1 (MQ=255) ttt‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc > 1:13242/1‑140 (MQ=255) ttt‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGGGAGTGTGCGGTTc > 1:91344/1‑140 (MQ=255) tt‑CCCCTTTGTCCGCGAAAAAAAAACGCAGTCCCCCTTTTGTCTGGGTTTAATATGTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCa < 1:887453/140‑1 (MQ=255) ttGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACtgctgc > 1:4856/1‑140 (MQ=255) tCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCAt > 1:931156/1‑140 (MQ=255) ccGCGAAAACATACCGCGGTCCCCTTTCTGTCTGGGTTTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATg < 2:278107/140‑1 (MQ=255) caTAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGa > 1:283567/1‑140 (MQ=255) tAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgc > 1:1156668/1‑140 (MQ=255) cGCAGTCCCCCTTCTGTCTGGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgcttt < 1:85618/140‑1 (MQ=255) gCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgcttt > 1:741536/1‑140 (MQ=255) tttCTGTCTGGGTTAAATGACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg < 1:1134509/139‑1 (MQ=255) gTCTGGGTTTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt < 2:125956/140‑1 (MQ=255) gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTTGCTTTCGCGGCttt > 2:809252/1‑140 (MQ=255) gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt > 1:141352/1‑140 (MQ=255) ggTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct > 1:125956/1‑140 (MQ=255) gtgtAATAACGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTTGCTTTCGCGGCTTTTTGctc < 1:809252/140‑1 (MQ=255) gtgtAATAACGTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctc < 1:1010869/140‑1 (MQ=255) atttttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCt < 1:46187/137‑1 (MQ=255) ttttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg < 2:302750/138‑1 (MQ=255) tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa > 2:352231/1‑140 (MQ=255) tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa > 2:972736/1‑140 (MQ=255) tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTg > 2:626214/1‑124 (MQ=255) tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg > 2:688177/1‑140 (MQ=255) tgtttTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTg < 1:626214/124‑1 (MQ=255) | ATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTT‑CGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATG > NC_000913/4111621‑4111885 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |