Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,032,731 | 0 | T | G | 28.3% | 115.1 / 18.5 | 60 | G80G (GGT→GGG) | yigZ | IMPACT family member YigZ |
Reads supporting (aligned to +/- strand): ref base T (14/29); new base G (17/0); total (31/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.19e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTG > NC_000913/4032596‑4032866 | gAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAAc < 2:336564/140‑1 (MQ=255) ttgttgAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTc > 1:862252/1‑140 (MQ=255) gttgAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGc < 2:1193870/140‑1 (MQ=255) tCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATg < 2:924184/129‑1 (MQ=255) tCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATg > 1:924184/1‑129 (MQ=255) gCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCAGCTAATgg > 1:168045/1‑140 (MQ=255) ggCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGcc < 1:73199/127‑1 (MQ=255) ggCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGcc > 2:73199/1‑127 (MQ=255) gCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAg > 1:77786/1‑140 (MQ=255) aaCACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAgcggc > 1:1032683/1‑140 (MQ=255) acacCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAgcggcg > 2:296891/1‑140 (MQ=255) cccGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCCCGCCCCGCTAATGGGGAGCGGGGGCgg > 2:1254978/1‑140 (MQ=255) ccGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCggg > 1:568849/1‑140 (MQ=255) cGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTAc < 2:568849/140‑1 (MQ=255) cGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTAc < 1:1172496/140‑1 (MQ=255) ggCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCggggg > 1:1134097/1‑132 (MQ=255) ggCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGa < 2:1134097/133‑1 (MQ=255) ggCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCCGCTAATGGGCAGCGGCGGCGGGGGAAATAcc > 2:214184/1‑140 (MQ=255) gCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCg < 2:1031324/140‑1 (MQ=255) aCCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGAGGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCt < 1:1465298/140‑1 (MQ=255) aCCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCt < 2:862252/140‑1 (MQ=255) ttGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGAGGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCCGCTAATGGGCAGCGGGGGCGGGGGAATTACCGCTGTgg > 2:1465298/1‑140 (MQ=255) ttGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCCGCTAATTGGCAGCGGCGGCGGGGGAAATAACCCTGTgg > 2:614516/1‑140 (MQ=255) ttGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTgg > 1:507430/1‑140 (MQ=255) tGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGt > 1:1220945/1‑140 (MQ=255) gcgtggcgtgGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTa < 1:476349/140‑1 (MQ=255) cgtggcgtggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAg < 2:127604/140‑1 (MQ=255) cgtggcgtggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGGAATTACCGCTGTGGGAg > 1:442476/1‑140 (MQ=255) tggcgtggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCg < 2:1054416/129‑1 (MQ=255) tggcgtggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCg > 1:1054416/1‑129 (MQ=255) tggGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGcta > 1:124817/1‑140 (MQ=255) ggTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGctact > 1:1123426/1‑140 (MQ=255) gTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGctacta < 2:1037886/140‑1 (MQ=255) gcgGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTAcgg < 1:463925/140‑1 (MQ=255) cgGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTAcggc < 2:1123426/140‑1 (MQ=255) gggTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTAcggcg < 1:296891/140‑1 (MQ=255) gggTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGGAAAACGATGCTCGCCCAGCTAATGGGGAGCGGGGGCGGGGGAAATACCGCTGTGGGAGTGCGCTACTAcggcg > 2:1567355/1‑140 (MQ=255) tGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACTGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCa < 2:168045/140‑1 (MQ=255) cgcCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCata < 2:507430/140‑1 (MQ=255) tgatTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTg < 2:124817/140‑1 (MQ=255) tctcAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCAc < 1:1175876/140‑1 (MQ=255) tcAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCACCg < 2:77786/140‑1 (MQ=255) tcAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCACCg < 2:567215/140‑1 (MQ=255) aaCAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCGCCGGt < 1:198635/140‑1 (MQ=255) aaCAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGt < 2:882515/140‑1 (MQ=255) cTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCCGCTAATGGGCAGCGGCGTCGGGGAAAATACCGCTGTGGGAGTGCGCTACTACGGCGGCATATTGCTTGGCACCgggggggt > 2:1133989/1‑138 (MQ=255) cTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGtt > 1:1278156/1‑140 (MQ=255) gTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTg < 1:208043/140‑1 (MQ=255) tctcTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGaaa < 2:1422163/140‑1 (MQ=255) tctcTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGaaa < 1:1254978/140‑1 (MQ=255) ctctGACGACGGGGAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGGAATTACCCCTGTGGGAGTGCGCTACTACGGCGGCATATTTCTTGGGACCGGGGGGGTTGTGgaaa > 2:1011408/1‑139 (MQ=255) ctGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCg > 1:1512570/1‑140 (MQ=255) acgacgGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGGTAGTGAAAGCGTAt > 1:1018284/1‑140 (MQ=255) cgggAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATggcgg < 1:446320/139‑1 (MQ=255) gggAGCCGGCGGGAACGGCAGGTAAACCCATGCTCGCCCAGCTAATGGGGAGCGGCGTCGGGGGAATTACCGCTGTGGGAGGGCGCTACCACGGCGGCATATTGCTTGGCACCGGGGGGGTAGGGGAAACGGATgggggg > 2:631894/1‑139 (MQ=255) gggAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCCGCTAATGGGCAGCGGGGTCGGGGGAATTACCGCTGTGGGAGTGCGCTACTACGGCGGCATATTGCTTGGGACCGGGGGGGTTGTGAAAGCGGATggcggg > 2:884188/1‑139 (MQ=255) gggAGCCGGCGGGAACGGCAGGGAAACCGATGCTCGCCCAGCTAATGGGGAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTTGCACCGGGGGGGTAGGGGAAGCGTATggcggg > 2:1057/1‑139 (MQ=255) gggAGCCGGCGGGAACGGCAGGGAAAACGATGCTCCCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGG‑GGGGGAGTGAAAGCGTATgggggcg > 2:215644/1‑140 (MQ=255) ggAACGGCAGGGAAAACCATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGGAAATACCGCCGTGGGAGTGCGCTACTACGGCGGCATATTTCTTTGCACCGGGGGGGTAGTGAAAGCGGATTgcggcggcggggaaca > 2:1283237/1‑133 (MQ=255) cAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTg > 1:1565939/1‑140 (MQ=255) | GAAAGCGTTTGTTGAATCGGTGCGGGCAGAACACCCCGATGCCCGGCACCATTGCGTGGCGTGGGTCGCGGGTGCGCCGGATGATTCTCAACAGCTGGGTTTCTCTGACGACGGGGAGCCGGCGGGAACGGCAGGTAAACCGATGCTCGCCCAGCTAATGGGCAGCGGCGTCGGGGAAATTACCGCTGTGGTAGTGCGCTACTACGGCGGCATATTGCTTGGCACCGGTGGGTTAGTGAAAGCGTATGGCGGCGGCGTGAATCAGGCGCTG > NC_000913/4032596‑4032866 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |