Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,085,594 | A→G | 14.6% | intergenic (+303/‑292) | 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,594 | 0 | A | G | 14.6% | 73.8 / 5.0 | 34 | intergenic (+303/‑292) | PP_1864/PP_1865 | hypothetical protein/transposase |
Reads supporting (aligned to +/- strand): ref base A (18/11); new base G (3/2); total (21/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01 |
ACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGAC‑CCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTT > NC_002947/2085447‑2085737 | aCGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGAC‑CCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGa < 8:34843/148‑1 (MQ=34) gATGCCGGTCAAGCGCTGGGCGCATGGTGGGGGCGTGCTGGGGGGTCCTCCGGGGAGCAAGCCGGTAGTGTGGGGAC‑CCCGTAGCGTTGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACCTTGGCCACGACCAAAGTAAc < 4:25022/149‑1 (MQ=11) gcTGGGCCCCTGGATGGGGCCGGCAGCGCGGGCCTCCCGGGAGCAAGCCCGGAGTGAGGGGAC‑CCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCg < 3:50352/149‑1 (MQ=255) cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGAC‑CCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGt > 4:329430/1‑149 (MQ=34) gATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGAC‑CCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGaa > 3:186370/1‑147 (MQ=34) gggCCTGCAGCGCGGTCCTCACGGGAGCAAGCCCGGAGTGAGGGGAC‑CCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGTCCACGACCAAAGTAACC‑‑CGCCGGAAGGGAGGAAAGGTGAATAGgct > 2:56627/1‑148 (MQ=255) ccTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGAC‑CCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCt > 7:297244/1‑149 (MQ=34) tGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGAC‑CCCGGAGCGTAGCGTAGGGGCCGGCTGAGGGGAGCCGCCGGCATTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGATTAGGCGCCTg > 7:330436/1‑148 (MQ=255) gcgGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGAC‑CCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGt > 8:76587/1‑148 (MQ=37) tccCGGGAGCAAGCCCGGAGTGGGGGGAC‑CCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGt > 7:55996/1‑141 (MQ=37) ccGGGAGCAAGCCCGGTGTTTGGGGAC‑CCCGGAGCGTAGCGTAGGGGCCGGAGGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTgaatgaat < 8:334756/147‑1 (MQ=255) ccGGGAGCAAGCCCGGAGTGAGGGGACGCCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGaatgaatg > 1:17172/1‑149 (MQ=255) aaGCCCGTTGTTATGGTAC‑CCCGGATCGTAGCGTAGGGGCCGGATGAGGGGAGCCGAAGGCTTTGTTTACTTTGTCCACTACCAAAGTAACA‑‑CGCCGGAATGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGa < 1:56627/149‑1 (MQ=11) aGCCCGGAGTGAGGGGAC‑CCCGGCGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGAt < 2:17172/149‑1 (MQ=37) ccGGAGTGAGGGGAC‑CCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGAtgt > 2:264315/1‑148 (MQ=38) tGAGGGGAC‑CCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAgc < 1:264315/148‑1 (MQ=39) aGGGGAC‑CCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGca > 2:154986/1‑149 (MQ=39) aGGGGAC‑CCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGca > 4:21682/1‑149 (MQ=39) ggggAC‑CCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGGCCAAAGTACCC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGGGCCGGCATCGTGAATGAAGGTCTACTTGAGGTGATCGCGca > 4:77169/1‑148 (MQ=14) cccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATa > 2:332472/1‑98 (MQ=37) cccGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATa < 1:332472/98‑1 (MQ=37) gcgtggGGCCCGGAGGTGGGCCCCCGACGGGTTTGGTAACTTTCCCCACCACAAAATTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCcaaaa < 4:186370/149‑1 (MQ=11) gggAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAgcgc > 3:344667/1‑149 (MQ=255) gAGCCGACGGCTTTTGTTACTTTGGCCACGACCAAAGTAACC‑‑CCCCGGAAGGGCGGAAATGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTTAAAGCGCAACACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCg > 2:186017/1‑149 (MQ=17) ccGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGtt < 3:21682/149‑1 (MQ=255) gACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACACCCCAAAGGCCCAACAGGTCAAAAGCTTTAGCGCCGGTTCt > 3:189640/1‑149 (MQ=255) tttGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCgg < 8:322165/149‑1 (MQ=255) tttGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGAGCCGGTTCTGATCgg > 7:159564/1‑149 (MQ=255) ttACTTTGGTCGTGGCCAAAGTAACCAAAGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCg > 3:16820/1‑149 (MQ=18) ttACTTTGGTCGTGGCCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACaa < 4:362711/98‑1 (MQ=11) ttACTTTGGTCGTGGCCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACaa > 3:362711/1‑98 (MQ=11) aCTTTGGTCGTGGCCAAAGTAACCAAAGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGc < 4:16820/146‑1 (MQ=18) cTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACa < 4:344667/148‑1 (MQ=255) cTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAg > 2:44089/1‑149 (MQ=255) ccACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGtt > 8:120336/1‑149 (MQ=255) | ACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGAC‑CCCGGAGCGTAGCGTAGGGGCCGGATGAGGGGAGCCGACGGCTTTGGTTACTTTGGCCACGACCAAAGTAACC‑‑CGCCGGAAGGGCGGAAAGGTGAATAGGCGCCTGCATCGTGAATGAATGTCTACTTGATGTGAAAGCGCACAACCCAAAGGCCCAAAAGCTCAAAAGCTTTAGCGCCGGTTCTGATCGGGCGTACAGTCCGTT > NC_002947/2085447‑2085737 |
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: ‑ |