Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,315,936 | 0 | T | A | 26.7% | 93.3 / 11.9 | 45 | K34* (AAA→TAA) | yciG | protein YciG |
Reads supporting (aligned to +/- strand): ref base T (26/7); new base A (0/12); total (26/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.75e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAATATG > NC_000913/1315799‑1316074 | cATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaa > 1:809978/1‑140 (MQ=255) gCCCGCAGATGAAATTATTTTTATCAGTTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTCCCCCCGCTTTGTTGACCGCTTTTTTTCCCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAc < 2:809978/140‑1 (MQ=255) cccGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAcc > 1:779661/1‑140 (MQ=255) cGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccgc > 2:511163/1‑140 (MQ=255) gCAGATGCATTTTTCTTTATCAGCTGATTTTAATAATTACCGGTGAATCAGGTTTTGCTTGTTTTATTCCCCCCGTTTTGTGGCCCGCTTTTTTTCCCCGCTTCAGATGCGCGTTGCGGTTCTTTTTAAAAATTAccgcc < 2:739727/140‑1 (MQ=255) aTCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCttt > 1:77245/1‑140 (MQ=255) ttGATATTAATAATAAGCGGTGAATCAGTTTTGGCCTGTTTTATCCCCCCCGTTTTGTTGCCCGCTTTTTTTCCCCGCTTCAGATGCGCGTGGCGGATCTTTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGc < 2:271072/139‑1 (MQ=255) gATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCt > 2:106930/1‑140 (MQ=255) aatatTAATAATTACGGGTAAATCAGGATTTGCCTGATTTATTCCCCCCTTTTTGTTGCCCGCTTTTTTTCCCCGCTTCAGATGCGCGTTGCGGATCTTTTTAAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCt < 2:595328/139‑1 (MQ=255) aataatTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCg > 1:890475/1‑140 (MQ=255) ataatTAATTGTAAATAAGGTTTTCCTTGTTTTATTCCCCCCTTTTTGTTGCCCGCTTTTTTTCCCCCCTTAAGAGGCGCGTTGGGGATTTTTTTAAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCgg < 2:995108/140‑1 (MQ=255) aaatTACGGGTAAATGAGTATTTGCCTTTTTTATTCCCCCCGTTTTGTTCCCCCCTTTTTTCCCCCGCTTCAGATGCGCTTTGCGAATCTTTTTAAAAATTCCCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGa < 2:455944/139‑1 (MQ=255) tAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGcc > 1:332028/1‑140 (MQ=255) gaaaTCAGGTTTTGCTTGTTTTATTCCCCCCGTTTTGTTGACCCCTTTTTTCCCCCGCTTCAGATGCGCGTTGCGGATTTTTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCa < 2:332028/138‑1 (MQ=255) aaTCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACg > 2:460579/1‑140 (MQ=255) ccTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGaa > 1:524888/1‑140 (MQ=255) aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAAttt > 2:988576/1‑140 (MQ=255) aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAAttt > 1:17776/1‑140 (MQ=255) ccGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAaccacc > 2:623354/1‑140 (MQ=255) ttttGTTGCCCGCTTTTTTCCCCCTCTTCAGATGCGCGTTGCGGATCTTTTTAAAAATTCCCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACg < 2:894129/139‑1 (MQ=255) tttGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGa > 1:772176/1‑140 (MQ=255) cccGCTTTTTTTCCCCGCTTCAGATGCGCGTTGCGGATTTTTTTAAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCgg < 2:779661/139‑1 (MQ=255) ccccTTTTTTCCCCGCTTCAGATGCGCGTTGGGGATCTTTTTAAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGTGAAATTTCCTGAACCACCACGATGTTCGGcc < 2:67327/138‑1 (MQ=255) gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCa > 1:347028/1‑140 (MQ=255) ccTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCAt > 1:111321/1‑140 (MQ=255) tttCCCCGCTTCAGAGGCGCGTTGCGGTTCTTTTTAAAATTTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTAt < 2:564721/140‑1 (MQ=255) ttcccccGCTTCAGATGCGGGTTGCGAATCTTTTTAAAAATTCCCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTAt < 2:212840/136‑1 (MQ=255) gCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTccc > 1:564721/1‑140 (MQ=255) tttCAGTTGCGCTTGGCGGTTCTTTTTAAAATTTCCCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcg < 2:111321/139‑1 (MQ=255) tttCAGATGCGTTTTGCGGATTTTTTTAAAATTTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcg < 2:966502/139‑1 (MQ=255) cTTCAGATGCGCTTTGGGGATTTTTTTAAAATTTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcg < 1:511163/140‑1 (MQ=255) cTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcg > 2:366113/1‑140 (MQ=255) cTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcg > 2:71859/1‑140 (MQ=255) ttCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcgt < 1:988576/140‑1 (MQ=255) ttCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcgt > 2:561833/1‑140 (MQ=255) ttCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcgt > 2:637183/1‑140 (MQ=255) cgcgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCgttgcgttgc > 1:670750/1‑140 (MQ=255) cgcgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCgttgcgttgc < 2:890475/140‑1 (MQ=255) cgcgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCgttgcgttgc > 2:304872/1‑140 (MQ=255) cgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCAt > 1:77016/1‑140 (MQ=255) cgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCAt > 2:871845/1‑140 (MQ=255) cgTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCAt > 1:862685/1‑140 (MQ=255) tGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGt > 1:273064/1‑140 (MQ=255) ggATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTc < 1:106930/140‑1 (MQ=255) aTCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCAt < 1:366113/140‑1 (MQ=255) ttttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAata < 1:623354/140‑1 (MQ=255) ttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAATATg < 1:851739/140‑1 (MQ=255) | CATACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGGATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAATATG > NC_000913/1315799‑1316074 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |