Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 557,339:1 | +A | 14.3% | 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 | 14.3% | 107.2 / 16.7 | 35 | intergenic (+79/‑130) | rpsN/rpsH | 30S ribosomal protein S14/30S ribosomal protein S8 |
Reads supporting (aligned to +/- strand): ref base . (15/15); new base A (3/2); total (18/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01 |
TCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAG > NC_002947/557209‑557480 | tCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGCCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGTCCCTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTg > 1:250750/1‑150 (MQ=255) gTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGc > 1:270100/1‑150 (MQ=255) ccAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCtt < 2:70814/129‑1 (MQ=255) ccAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCtt > 1:70814/1‑129 (MQ=255) cAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCAc < 2:262863/150‑1 (MQ=255) aGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCAcc > 2:76857/1‑150 (MQ=255) gTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTg > 2:221042/1‑150 (MQ=255) cAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTgg < 1:221042/150‑1 (MQ=255) cAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTgg < 1:399080/150‑1 (MQ=255) cAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGAAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTgg < 1:231760/150‑1 (MQ=255) aaGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTggg > 1:401786/1‑150 (MQ=255) ggCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCt < 2:663946/110‑1 (MQ=255) ggCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCt > 1:663946/1‑110 (MQ=255) ctcctGGTAATCGTTGCAGGCTAGGGGAGGGTGGGAGGACAGAGATCTTCCGTCGTGCGGTGACTTTTCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGttttt < 2:250750/150‑1 (MQ=255) ctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGa > 2:225422/1‑135 (MQ=255) ctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGa < 1:225422/135‑1 (MQ=255) gTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCaa > 1:523598/1‑150 (MQ=255) gCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGc > 1:109775/1‑133 (MQ=255) gCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGc < 2:109775/133‑1 (MQ=255) cAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGggtggt < 2:270100/150‑1 (MQ=255) aaGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGtt > 2:280280/1‑150 (MQ=255) gAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCa > 2:598068/1‑150 (MQ=255) aGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCtaa < 1:410780/150‑1 (MQ=255) ggTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAg < 1:280280/150‑1 (MQ=255) ttGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCt > 1:798143/1‑132 (MQ=255) ttGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCt < 2:798143/132‑1 (MQ=255) aaacgaaacaaGCCCC‑ATATAGGGGCTTGTTT‑GTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGc > 1:628000/6‑110 (MQ=37) aaacgaaacaaGCCCC‑ATATAGGGGCTTGTTT‑GTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGc < 2:628000/105‑1 (MQ=37) caacaaACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCa < 1:124426/150‑1 (MQ=255) acaaacaaGCCCCTATAT‑GGGGCTTGTATCGTGTCTGCCTTATGTCTAGAATGACAGTCTCACCTTAGCCTTGGTTTTGGCCCAAACGGATCGGGTGGTTGTCATAGCCGCACGGCTAATAATTCTTGTATCAGGAGCATCTATCCCATg > 2:13991/1‑150 (MQ=255) caaacaaGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGa < 1:76857/150‑1 (MQ=255) aaacaaGCCCC‑ATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGa > 2:422081/1‑149 (MQ=255) aaacaaGCCCC‑ATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGa < 1:422081/149‑1 (MQ=255) acaaGCCCC‑ATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTa > 2:594396/1‑150 (MQ=255) ccccTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAg > 1:651891/1‑150 (MQ=255) | TCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATAT‑GGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAG > NC_002947/557209‑557480 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |