| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 1,241,973 | Δ1 bp | 10.4% | intergenic (+32/+26) | grxD → / ← bfr‑II | monothiol glutaredoxin/bacterioferritin |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 1,241,973 | 0 | A | . | 10.4% | 144.3 / 11.5 | 48 | intergenic (+32/+26) | grxD/bfr‑II | monothiol glutaredoxin/bacterioferritin |
| Reads supporting (aligned to +/- strand): ref base A (19/24); new base . (3/2); total (22/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.49e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01 | |||||||||||
CCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGG > NC_002947/1241832‑1242104 | ccGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGttt < 1:1028606/150‑1 (MQ=255) aGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTAt < 1:233826/150‑1 (MQ=255) tGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATggg < 2:405141/150‑1 (MQ=255) gcggcAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAAt > 1:118638/1‑150 (MQ=255) gcggcAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAAt > 1:118668/1‑150 (MQ=255) cggcAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATa < 2:118638/150‑1 (MQ=255) cggcAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATa < 2:118668/150‑1 (MQ=255) cGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAAGAGTAAct < 1:607569/150‑1 (MQ=255) catGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGc > 1:243190/1‑150 (MQ=255) aaTGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATgggg < 1:958302/123‑1 (MQ=255) aaTGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATgggg > 2:958302/1‑123 (MQ=255) gTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACt < 1:303064/150‑1 (MQ=255) ttCGAGAAGGGTGAGCTGTAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTg > 2:970979/1‑150 (MQ=255) aaGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTa < 1:896120/150‑1 (MQ=255) tGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTtct < 2:243190/150‑1 (MQ=255) gCTGTAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGa < 1:970979/150‑1 (MQ=255) tGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCg > 2:279717/1‑105 (MQ=255) tGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCg < 1:279717/105‑1 (MQ=255) tGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATa > 1:855252/1‑150 (MQ=255) gCAAACCCTGATCACGGCGGCGGCCGCGAAGACCAAGTCTGCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATAc < 2:131959/150‑1 (MQ=255) cAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGATTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGCTTCGGATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATAcc > 1:979591/1‑150 (MQ=255) tGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCTGATC‑GCGGGTTTTTTTa < 2:414150/89‑1 (MQ=37) tGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGCTGATC‑GCGGGTTTTTTTa > 1:414150/1‑89 (MQ=37) gATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGt < 1:482261/150‑1 (MQ=255) tCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCg > 2:840448/1‑150 (MQ=255) caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTggg > 1:1004390/1‑150 (MQ=255) gccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGAc > 1:595821/1‑100 (MQ=255) gccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGAc < 2:595821/100‑1 (MQ=255) ccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTcc > 2:1063906/1‑150 (MQ=255) caaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGc < 2:190542/53‑1 (MQ=255) caaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGc > 1:190542/1‑53 (MQ=255) caaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTcca > 1:883138/1‑150 (MQ=255) aaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCTAATAGTTACTCTTCGCCCATCTGCGACTGCAg > 2:923363/1‑102 (MQ=255) aaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCTAATAGTTACTCTTCGCCCATCTGCGACTGCAg < 1:923363/102‑1 (MQ=255) ggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTccagcc < 2:1004390/150‑1 (MQ=255) gcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTccagcca < 1:440971/150‑1 (MQ=255) ttCTGAAGCCTGATCTTCCGGAAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGt < 2:855252/150‑1 (MQ=255) tGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTccag < 2:1118226/142‑1 (MQ=255) tGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTccag > 1:1118226/1‑142 (MQ=255) gCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGAtgtg > 2:908399/1‑150 (MQ=255) tCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTc < 2:883138/150‑1 (MQ=255) ggCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGa > 1:943263/1‑150 (MQ=255) cccATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGt > 1:99351/3‑127 (MQ=255) cccATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGt < 2:99351/125‑1 (MQ=255) gCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGa < 1:840448/150‑1 (MQ=255) aTAAAAAAACCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAgg > 2:526921/1‑150 (MQ=255) aaaCCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTActc < 2:175248/44‑1 (MQ=14) aaaCCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTActc > 1:175248/1‑44 (MQ=14) aGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAg < 2:900203/131‑1 (MQ=255) aGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAg > 1:900203/1‑131 (MQ=255) | CCGGTGAGCTGGTTGGCGGCAGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGG > NC_002947/1241832‑1242104 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |