Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 557,339:1 | +A | 15.4% | intergenic (+79/‑130) | rpsN → / → rpsH | 30S ribosomal protein S14/30S ribosomal protein S8 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 557,339 | 1 | . | A | 15.4% | 121.2 / 21.3 | 39 | intergenic (+79/‑130) | rpsN/rpsH | 30S ribosomal protein S14/30S ribosomal protein S8 |
Reads supporting (aligned to +/- strand): ref base . (15/18); new base A (3/3); total (18/21) | |||||||||||
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.91e-01 |
TGGGCCGTAACATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAGGACCCG > NC_002947/557190‑557486 | tGGGCCGTAACATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTatat < 2:588288/150‑1 (MQ=255) gctgcGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTgtttcg < 2:584033/150‑1 (MQ=255) tGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTa < 2:964041/150‑1 (MQ=255) ggggcgGTGGCGTACCAGGTATGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAAGCTTCCGACCGCCGGTGATCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTa < 2:186942/147‑1 (MQ=255) gcgGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGaa > 2:1059664/1‑150 (MQ=255) cctcctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCC‑ATATAGGGGCTTgtttgttg < 2:87114/105‑5 (MQ=38) cctcctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCC‑ATATAGGGGCTTgtttgttg > 1:87114/1‑101 (MQ=38) ctcctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTtgt < 1:1073354/114‑1 (MQ=255) ctcctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTtgt > 2:1073354/1‑114 (MQ=255) ctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGcc > 2:1116504/1‑150 (MQ=255) tGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTgg < 2:642605/130‑1 (MQ=255) tGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTgg > 1:642605/1‑130 (MQ=255) ggCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGggtg < 1:223663/150‑1 (MQ=255) aaGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGtt > 2:581443/1‑150 (MQ=255) ggCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCAGGTGAGccgg < 2:684965/120‑3 (MQ=255) ggCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCAGGTGAGccgg > 1:684965/1‑118 (MQ=255) cgcgGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGTGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCa < 1:1059664/150‑1 (MQ=255) ggTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGc < 1:1116504/150‑1 (MQ=255) gTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGcc > 2:1182283/1‑150 (MQ=255) gggAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTgtt < 1:520849/62‑1 (MQ=255) gggAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTgtt > 2:520849/1‑62 (MQ=255) gggAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGc < 1:1049708/150‑1 (MQ=255) gAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTACCCAAACAAATCGGGTGGTTGTCATAGCCGCAAGGCta > 2:561587/1‑150 (MQ=255) cGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTTGTTTCCGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGt < 1:1113641/150‑1 (MQ=255) cGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTg < 2:6066/149‑1 (MQ=255) cGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTg > 1:6066/1‑149 (MQ=255) ccTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCt < 1:839866/150‑1 (MQ=255) cTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCAACTGAAGCTGGGTTTTTGCCCAAACGGTTAGAAGAGTGTGCAGAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTa > 1:570850/1‑150 (MQ=255) aacaaacaAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAg > 2:35421/1‑115 (MQ=255) aacaaacaAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAg < 1:35421/115‑1 (MQ=255) gaaacaaGCCCC‑ATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCtt > 2:419813/2‑126 (MQ=255) gaaacaaGCCCC‑ATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCtt < 1:419813/125‑1 (MQ=255) aaacaaGCCCC‑ATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGttttt > 2:773343/1‑77 (MQ=37) aaacaaGCCCC‑ATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGttttt < 1:773343/77‑1 (MQ=37) gCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGtt < 1:884911/83‑1 (MQ=255) gCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGtt > 2:884911/1‑83 (MQ=255) gCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGt > 1:848806/1‑82 (MQ=255) gCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGt < 2:848806/82‑1 (MQ=255) cTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAGGAc < 1:581443/150‑1 (MQ=255) cTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAGGAc > 2:598277/1‑150 (MQ=255) tat‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAGGACCCg > 2:1009034/1‑150 (MQ=255) | TGGGCCGTAACATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAGGACCCG > NC_002947/557190‑557486 |
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: ‑ |