| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 2,085,589 | C→T | 12.2% | intergenic (+298/‑297) | PP_1864 → / → PP_1865 | hypothetical protein/transposase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 2,085,589 | 0 | C | T | 12.2% | 101.2 / 5.6 | 41 | intergenic (+298/‑297) | PP_1864/PP_1865 | hypothetical protein/transposase |
| Reads supporting (aligned to +/- strand): ref base C (21/15); new base T (3/2); total (24/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01 | |||||||||||
AGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCCAAAAG‑CTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTT > NC_002947/2085442‑2085737 | aggTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTACGCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGc < 1:195587/149‑1 (MQ=21) gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGAcc > 6:272110/1‑148 (MQ=34) gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCa < 5:239485/149‑1 (MQ=34) ggATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGCCCAAAGTaa > 8:88720/1‑149 (MQ=21) ggATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GAGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTaa > 6:54473/1‑149 (MQ=21) gATGCCCGTCAAGCGGTGGGCCCCTGGAGGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTaa < 2:159237/148‑1 (MQ=255) gATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGGGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAAc < 2:216530/149‑1 (MQ=21) gATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGTCGTGGCCAAAGTaa < 5:98845/148‑1 (MQ=255) gATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAAc > 8:141138/1‑149 (MQ=34) cccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACc‑‑cgc < 1:56044/149‑1 (MQ=34) ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACc‑‑cgcc > 6:132445/1‑149 (MQ=34) ccccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAgg > 2:140390/1‑149 (MQ=34) aGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGaa > 8:223550/1‑149 (MQ=34) tGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTg > 6:80471/1‑148 (MQ=35) gcgcGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCg < 5:272110/149‑1 (MQ=37) tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGgcgc > 3:212351/1‑132 (MQ=34) tccCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGgcgc < 4:212351/132‑1 (MQ=34) aaGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑GGCCGGAAGGGCGGAAAGGTGAATAGGCGGCTGCCTCGTGAATGAAGGTCTACATGa > 6:224594/1‑149 (MQ=11) aaGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCAGGATGAGGGGAGCCGTGCGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACtt > 1:208986/1‑148 (MQ=255) aGGGGACCCCGGAGCGTA‑GCTTAGGGGCCGGATGTGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAgcgc < 7:141138/148‑1 (MQ=255) gACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAAc > 4:225755/1‑149 (MQ=39) cGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCaaa < 6:77113/149‑1 (MQ=39) cGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCa > 6:112156/1‑147 (MQ=39) ggAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAg < 5:132445/149‑1 (MQ=40) gAGCGTA‑GCGTTGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAgg < 8:240081/149‑1 (MQ=255) gCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAGCCCAAAGGCCcaaaag‑ctcaa > 5:44351/1‑147 (MQ=255) gAGGGGACCCG‑ACGTGTTTGGTTACTTTGCCCACGACCAAATTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCCAAAAG‑CTCAAAAGCTTTa < 7:88720/148‑1 (MQ=14) ccG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCCAAAAG‑CTCAAAAGCTTTAGCGCCGGtt > 5:150776/1‑149 (MQ=255) cG‑ACGGCTTTGGTTACTTTGGCCACGACCAAATTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCCAAAAG‑CTCAAAAGCTTTAGCGCCGGTTc < 2:208986/149‑1 (MQ=255) g‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCCAAAAG‑CTCAAAAGCTTTAGCGCCGGTTCt > 4:139788/1‑149 (MQ=255) ggCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCCAAAAGCCTCAAAAGCTTTAGCGCCGGTTCTGa < 1:140390/149‑1 (MQ=255) cTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCcaaaag‑ctcaaaagct < 5:142768/130‑1 (MQ=255) cTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCcaaaag‑ctcaaaagct > 6:142768/1‑130 (MQ=255) gcgggTTACTTTGGTCGTGGCCAAAGTAACCAAAGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCCAAAAG‑CTCAAAAGCTTTAGCGCCGGTTCTGATc > 4:141510/4‑149 (MQ=17) gggTTACTTTGGTCGTGGCCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACtt > 1:150933/2‑84 (MQ=11) ttACTTTGGTCGTGGCCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTCGAAAGCGCACAACCCAAAGGCCCAAAAG‑CTCAAAAGCTTTAgcgc < 1:191641/132‑1 (MQ=14) ttACTTTGGTCGTGGCCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCCAAAAG‑CTCAAAAGCTTTAgcgc > 2:191641/1‑131 (MQ=255) ttACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGGCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCCAAAAG‑CTCAAAAGCTTTAGCCCCGGTTCTGATCGGGCGTa > 2:27257/1‑149 (MQ=255) cTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCCAAAAG‑CTCAAAAGCTTTAGCGCCGGTTCTGATCgg < 3:152065/141‑1 (MQ=255) cTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCCAAAAG‑CTCAAAAGCTTTAGCGCCGGTTCTGATCgg > 4:152065/1‑141 (MQ=255) tGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCCAAAAG‑CTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTcc > 3:101552/1‑149 (MQ=255) gcccACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCCAAAAG‑CTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTcc < 7:223550/146‑1 (MQ=255) gCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCCAAAAG‑CTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGt > 2:4575/1‑149 (MQ=255) ccACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCCAAAAG‑CTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGtt > 3:77705/1‑149 (MQ=255) | AGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTA‑GCGTAGGGGCCGGATGAGGGGAGCCG‑ACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGT‑GAAAGCGCACAACCCAAAGGCCCAAAAG‑CTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTT > NC_002947/2085442‑2085737 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |