Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,868,920 | G→C | 10.5% | intergenic (+416/‑515) | PP_2520 → / → PP_2522 | hypothetical protein/transposase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,868,920 | 0 | G | C | 10.5% | 106.9 / 4.5 | 38 | intergenic (+416/‑515) | PP_2520/PP_2522 | hypothetical protein/transposase |
Reads supporting (aligned to +/- strand): ref base G (16/18); new base C (2/2); total (18/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01 |
CAGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTT > NC_002947/2868780‑2869056 | cAGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCACCAATGCCCCGGCATTgat < 2:1073784/150‑3 (MQ=255) ccTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGTGCGTTGGTGAt > 2:63174/1‑150 (MQ=255) ttCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTg > 1:1155041/1‑140 (MQ=255) ttCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTg < 2:1155041/140‑1 (MQ=255) attggccattCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATatcaat < 2:959655/150‑1 (MQ=255) gAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTcgcg < 1:264720/150‑1 (MQ=255) aaGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCgg > 1:276984/1‑150 (MQ=255) gCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGcc < 2:304755/150‑1 (MQ=255) cATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGcccc > 2:1100603/1‑150 (MQ=255) tttGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGccccc < 2:989409/149‑1 (MQ=255) tttGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGccccc > 1:989409/1‑149 (MQ=255) gggAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTg < 2:1096892/107‑1 (MQ=255) gggAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTg > 1:1096892/1‑107 (MQ=255) tGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATgc < 2:276984/150‑1 (MQ=255) ccTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATatcaat < 1:309477/105‑1 (MQ=255) ccTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATatcaat > 2:309477/1‑105 (MQ=255) cTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCACCCACGTTGTCGATGCGCGATTTTGTcc > 1:171969/1‑150 (MQ=255) gggCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACgg < 1:1100603/150‑1 (MQ=255) gggCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACgg > 1:395999/1‑150 (MQ=255) aGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCttt < 2:395999/150‑1 (MQ=255) gCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGAtt < 1:836312/134‑1 (MQ=255) gCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGAtt > 2:836312/1‑134 (MQ=255) ggAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAAcc > 2:427727/1‑134 (MQ=255) ggAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAAcc < 1:427727/134‑1 (MQ=255) gAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACAtct > 1:531361/1‑136 (MQ=255) gAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACAtct < 2:531361/136‑1 (MQ=255) ggcCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGAta < 1:362112/150‑2 (MQ=255) ccGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCg < 1:1039790/150‑1 (MQ=255) aCCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATc > 1:218670/1‑150 (MQ=255) gTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTc < 1:875183/85‑1 (MQ=255) gTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTc > 2:875183/1‑85 (MQ=255) atcaccaATGCCCCGGCATTGATGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGttt > 2:123629/6‑150 (MQ=255) atcaATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACAtc < 2:699519/108‑1 (MQ=255) atcaATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACAtc > 1:699519/1‑108 (MQ=255) ccaATGCCCCGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGaa < 2:500491/136‑1 (MQ=255) ccaATGCCCCGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGaa > 1:500491/2‑137 (MQ=255) caATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCa > 2:252027/1‑91 (MQ=255) caATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCa < 1:252027/91‑1 (MQ=255) | CAGCCTAGTTCAGGCCATTGGCCATTCGGGTGAAGCCATTTGGGAGCGGCTGGGGCATGAGCCTTCGGGGCAGCCGTTCAACTCGTTGGTGCAGCTGGAGTTCGAGCGCGGGCGGCCGCGTAACCCGTTCATCAATGCCGGGGCATTGGTGATCTGCGATATCAATCAGTCGCGCTTCGCGGCCCCCACGTTGTCGATGCGCGATTTTGTCCGACGGCTTTCAGGCAACCCTCACATCTCGATCAATACCCGGGTCGCTGATTCGGAATATCAGTTT > NC_002947/2868780‑2869056 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |