Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,224,100 | C→T | 7.0% | W149* (TGG→TGA) | PP_1962 ← | site‑specific recombinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,224,100 | 0 | C | T | 7.0% | 164.0 / 2.3 | 57 | W149* (TGG→TGA) | PP_1962 | site‑specific recombinase |
Reads supporting (aligned to +/- strand): ref base C (24/29); new base T (2/2); total (26/31) | |||||||||||
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.78e-01 |
GAAAGAGGAACCCACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGG > NC_002947/2223963‑2224244 | gAAAGAGGAACCCACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCaa < 7:147289/149‑1 (MQ=255) aaaGAGGAACCCACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGTCAAAATCGCaa > 4:165085/1‑148 (MQ=255) aaaGAGGAACCCACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGTCAAAATCGCaa > 4:159601/1‑148 (MQ=255) gagGAACCCACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑At > 5:12903/1‑148 (MQ=255) gAACCCACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGaa > 7:21708/1‑148 (MQ=255) cccACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGaaaaa > 4:273998/1‑148 (MQ=255) cACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAg < 8:179519/149‑1 (MQ=255) aCGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAgg < 6:12903/149‑1 (MQ=255) aCGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAgg < 4:175501/149‑1 (MQ=255) ttGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCTAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGTGCGATAAGGATAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTg < 7:102042/149‑1 (MQ=255) aaTCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGTCAAAATCGCAAC‑ATGAAAAACAGGCCTGaa < 3:159601/148‑1 (MQ=255) aaTCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGTCAAAATCGCAAC‑ATGAAAAACAGGCCTGaa < 3:165085/148‑1 (MQ=255) gAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGaataat < 8:21708/148‑1 (MQ=255) gAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGaataat > 7:19232/1‑148 (MQ=255) aGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGa < 3:273998/149‑1 (MQ=255) cccTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACt > 4:336450/1‑149 (MQ=255) cTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCg > 6:99110/1‑149 (MQ=255) aaaTTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGc > 6:36386/1‑130 (MQ=255) aaaTTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGc < 5:36386/130‑1 (MQ=255) gTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTg > 5:104160/1‑149 (MQ=255) ccTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGTTAAGGAAAGCCAAAATCGCAACTTTGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCtttt < 1:255050/147‑1 (MQ=255) ttGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGtt > 2:8198/1‑130 (MQ=255) ttGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGtt < 1:8198/130‑1 (MQ=255) gtgGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGTAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGACAACTTGCTTTTGcgtcg > 7:18732/1‑148 (MQ=255) tgGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCgtcgtc > 4:138375/1‑149 (MQ=255) gcgcTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGa < 3:138375/149‑1 (MQ=255) aTCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGAtt > 8:302800/1‑147 (MQ=255) tCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTaa < 5:99110/148‑1 (MQ=255) tgatgTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCg > 3:223645/1‑149 (MQ=255) tgatgTGATATGTTCGGGTAGACTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAATAATGCCTGAATAATGAATCGGTTGTACACTTGCATGTGCTTATTCGAGTGAGCTATTAATCATCg > 8:335796/1‑149 (MQ=255) tgATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTc > 8:127287/1‑105 (MQ=255) tgATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTc < 7:127287/105‑1 (MQ=255) tgATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGAtt < 7:250249/148‑1 (MQ=255) gATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTg > 3:232199/1‑149 (MQ=255) atttGTTCGGGTAGCCTTAAACTTTCAGGAA‑TGGCTTGTTGAGCGATAAGGATAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGt < 4:94969/145‑1 (MQ=255) atGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGtt < 8:19232/148‑1 (MQ=255) tGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTa < 4:223645/148‑1 (MQ=255) tGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTa < 7:150511/148‑1 (MQ=255) ttCGGGTAGCCTTAAACTTTCAGGAATTAGGTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATga < 7:191696/149‑1 (MQ=255) gggTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGAt < 7:302800/149‑1 (MQ=255) cTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGt < 3:336450/148‑1 (MQ=255) ttAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGt < 3:315482/133‑1 (MQ=255) ttAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGt > 4:315482/1‑133 (MQ=255) ttAAACGTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGGTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGc > 4:315561/1‑149 (MQ=255) aCTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAg < 4:232199/149‑1 (MQ=255) tttCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAgg > 1:173754/1‑148 (MQ=255) aGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCg < 6:104160/149‑1 (MQ=255) gAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCt > 4:176873/1‑149 (MQ=255) gCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGTGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCa < 5:263862/134‑1 (MQ=255) gttgAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCAc < 1:111511/68‑1 (MQ=255) gttgAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCAc < 1:115461/68‑1 (MQ=255) gttgAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCAc > 2:111511/1‑68 (MQ=255) gttgAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCAc > 2:115461/1‑68 (MQ=255) ggAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCaa > 5:158191/1‑147 (MQ=255) gAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAg > 3:278812/1‑149 (MQ=255) aaaGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAgg < 3:176873/149‑1 (MQ=255) aaaGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGCCCCAATAgg < 3:315561/149‑1 (MQ=255) ccAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTggg < 8:117124/135‑1 (MQ=255) ccAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTggg > 7:117124/1‑135 (MQ=255) | GAAAGAGGAACCCACGATTGAATCGAGCCCTTGAAATTATCTAGGACGTCCTTGTGGCGCTCATGGCTATCGAATGATGTGATATGTTCGGGTAGCCTTAAACTTTCAGGAATTAGCTTGTTGAGCGATAAGGAAAGCCAAAATCGCAAC‑ATGAAAAACAGGCCTGAATAATGACTCGGTTGTCCACTTGCTTTTGCGTCGTCGAGGGCCCGATTAAGCATCGGATTTGTTATGAGATCCAAGTGCTCAGGAATCGCTGAAAGCCTATTGGGGCCCAATAGG > NC_002947/2223963‑2224244 |
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: ‑ |