Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,578,433 | 0 | G | T | 29.5% | 97.7 / 23.2 | 61 | intergenic (+110/‑214) | PP_3158/PP_5589 | hypothetical protein/rhs repeat‑associated core domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base G (17/26); new base T (18/0); total (35/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.95e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGT > NC_002947/3578298‑3578577 | tgatgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTTCAAATGGCGTTTTTTCCAGTTGGAGCgc > 7:150224/1‑148 (MQ=255) tgatgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTTCCAGTTGCAGCGa > 3:233826/1‑149 (MQ=255) gatgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTTAAAAATGCATTTTTTCCAGTTGCAGCGAg > 4:357867/1‑149 (MQ=255) atgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGGCAGTTGGAGCGAgt > 3:282901/1‑149 (MQ=255) aaCTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAgtgt < 4:62823/148‑1 (MQ=255) ccGGTGCATTAAGAGGCCTCTGCCGCACGCCGTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGcc < 7:152866/148‑1 (MQ=255) gTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCATCAACAGGCAATTTTTCCAGTTGCAGCGAATGTAAGACACCCTAAAAATc > 4:135331/1‑148 (MQ=255) aaGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa < 5:66740/148‑1 (MQ=255) agagGCCTCTGTCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCAGATTTTCAACTTTCAGCGAGTGTATGAAAGCCGAGAAATCAAgg > 2:388416/1‑149 (MQ=255) agagGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAgg < 3:357867/149‑1 (MQ=255) gagGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGc < 4:286875/111‑1 (MQ=255) gagGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGc > 3:286875/1‑111 (MQ=255) gCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGGATTTTTGGCAGTTGCAGCGAGTGGAt > 4:320194/1‑126 (MQ=255) tctGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAATGGCATTTTTGCCAGTTGCAGCGAGTGTATTAAAGCCTAGAATTCAAGGGAAATgg > 3:185606/1‑149 (MQ=255) cGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTGCTCTAAGGTCAtt < 3:135331/148‑1 (MQ=255) gCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCAGCAAAAAGCAATTTTTACAGTTTCAGCGAGTTTATGAACACCAAGAATTAAAGGGAAATGGTGTTAAAAAACGTCATTAAtatt > 5:409069/1‑145 (MQ=255) ccATACATGAATGAGGATTGAGTTGATGAGTGTTGATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTTCCAGTTGCAGCGAGTGGATGACAGGCTAGAATTAAAGGGAAATGGTGTTGTCTAAGGTCAGTAATAGTGa > 4:267704/1‑148 (MQ=255) atgaatgaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCAGCAACATGCAATTTTTACAGTTGCAGCGAGTGTATGACAGCCAAGAATTAAAAGGAAAAGTTGTTGTCTAAGGTCATTAATAATTAAACGTTc > 4:84455/1‑149 (MQ=255) tgaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTTCAAAATGCAATTTTTACACATTCACCGAGTGTATGACAGCCAAAAAATAAAAGGAAAAAGTGTTATAAAAAGTCAATAATAGGGAAACCTTCTCAAg > 2:233878/1‑149 (MQ=255) agtaTTTAGTTGATGGGTGTTCATGAGCTAATCAACGTTTCGTCGGTCTTCAACTTGCATTTTTGCCATTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGc < 1:189685/145‑1 (MQ=255) aGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTTGCAGTTGCAGCGAGTGTATGACAGCCTAGAATTAAAGGGAAATTGTGTTCTGTAAGGTAATTAATAGTGaaaa > 8:95981/1‑137 (MQ=255) aGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAAc < 7:95981/138‑1 (MQ=255) aGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGc < 7:389004/148‑1 (MQ=255) gATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTTCCAGTTGCAGCGAGTGTATTACAGCCTAGAATTAAAgg > 6:85511/1‑101 (MQ=255) gATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAgg < 5:85511/101‑1 (MQ=255) ttGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTTGCAATTTTGCCAGTTGCAGCGAGTGTATTACAGCCTAGAATTAAATGTAAATTGTGTTATATAAGGTAATTAATAGTGAAACGTTCGCAAGCGTaca > 6:419550/1‑149 (MQ=255) gAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAAAGGCATTTTTTCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTGTCTAAGGTCATTAATCGTGAAAAGTTCGCAAGCGTACACg > 8:319481/1‑149 (MQ=255) gTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGc < 4:233826/148‑1 (MQ=255) gtTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTa < 1:19772/142‑1 (MQ=255) gtTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTTGCAGTTGCAGCGAGTGTTTTTCAGCCAAGAATTAAAGGGTAAATGTGTTGTCTAAGGTGAATAATAGTGAAACGTTAGCAAGCGTACACCCCCTaaagaaag > 2:19772/1‑142 (MQ=255) ttCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTTCCAGTTGCAGCGAGTGGATGACAGCCGAGAATTAAAGGGAAATGGTGTTGTCTAAGGTCATTAATAGTGAAACGGTGGCAAGCCTACACGCCCTACGgcgcga > 8:185129/1‑147 (MQ=255) aTGAGCGAATCAACGTTTCGTCGGGCTGCTACTGGCATTTTTGCCAGTTGCAGCGATTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTAtt < 7:248092/148‑1 (MQ=255) gAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTTCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTAAAGGGAAATGGTGTTGTATAAGGTCATTAATAGTTAAAAGTTAGCAAGCGTACACCCCCTAACGCGCTAttttg > 4:86997/1‑147 (MQ=255) aGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTTCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATTGTGTTCTCTAAGGTCATTAATCCTGAAACGTTAGCAAGCGTACACGCCCTACGGCGCTATTTggg > 1:304286/1‑149 (MQ=255) aaTCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACCGCCTAGAATTCAAGGGAAATGGTGTTGTGTAAGGGCATTAATAGGGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGGGa > 8:81323/1‑149 (MQ=255) tCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATg < 8:150224/149‑1 (MQ=255) cAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCaa < 7:16193/114‑1 (MQ=255) cAACGTTTCGTCGGGCTGCAACATGCATTTTTTGCAGGTGGAGCGAGTGTTTTACAGCCTATAAATAAAGGGAAATTGTGTTATGTAATGTCATTTATAGTGAAAAGTTAGCaa > 8:16193/1‑114 (MQ=255) aaCGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATgg < 5:419550/148‑1 (MQ=255) aaCGTTTCGTCGGGCTGCAACAGGCATTTTTTGCAGTTGCAGCGAGTGTATTACAGCCTAGAAATAAAGGGTAAATGTGTTGTCTAAGGGGATTAATAGTGAAAAGTTAGCAATCGTAAACCCCCAACGGCGCGATTTTGGGTGATggg > 8:46738/1‑148 (MQ=255) aaCGTTTCGTCGGGCTGCAACAGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAAAGGTGTTCTGTAAGGTGATTAATCGTGAAACGTTAGCAAGCGTACACGCCCTACGGCGCTATTTTGGGTGATGGc > 6:318777/1‑149 (MQ=255) cGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCa < 3:84455/149‑1 (MQ=255) tttCGTCGGGCTGGAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGt < 8:82725/96‑1 (MQ=255) tttCGTCGGGCTGGAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGt > 7:82725/1‑96 (MQ=255) gtcgGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCt < 7:185129/148‑1 (MQ=255) gggCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCg < 3:267704/148‑1 (MQ=255) gCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCgtgt > 4:55238/1‑149 (MQ=255) cTGCAACTGGCAAAAATGCCAGTTGCAGCGAGTGTTTGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGtt > 6:269617/1‑149 (MQ=255) tGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGtt > 2:67037/1‑148 (MQ=255) gCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTg > 3:413603/1‑148 (MQ=255) cAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGaa < 4:282901/149‑1 (MQ=255) aaCTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGATATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAg > 2:378870/1‑149 (MQ=255) aaCTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAg < 7:46738/149‑1 (MQ=255) aCTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGt < 3:55238/149‑1 (MQ=255) tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTc > 2:220381/1‑70 (MQ=255) tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTc < 1:220381/70‑1 (MQ=255) tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTa > 5:236078/1‑148 (MQ=255) tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTa > 5:230134/1‑148 (MQ=255) tttGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGt < 3:86997/148‑1 (MQ=255) atGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGtt < 3:267324/133‑1 (MQ=255) atGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGtt > 4:267324/2‑134 (MQ=255) | TGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGT > NC_002947/3578298‑3578577 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |