Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 529,930 | G→A | 11.4% | intergenic (+24/‑160) | PP_5SC → / → birA | 5S ribosomal RNA/bifunctional biotin‑[acetyl‑CoA‑carboxylase] ligase/biotin operon repressor BirA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 529,930 | 0 | G | A | 11.4% | 87.8 / 3.6 | 35 | intergenic (+24/‑160) | PP_5SC/birA | 5S ribosomal RNA/bifunctional biotin‑[acetyl‑CoA‑carboxylase] ligase/biotin operon repressor BirA |
Reads supporting (aligned to +/- strand): ref base G (15/16); new base A (2/2); total (17/18) | |||||||||||
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.71e-01 |
GACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGCCCATGG > NC_002947/529792‑530054 | taaggcaggcatgatagccatgggctgttggcttgcataggaaaagggcctgcaggagaggccgcgtatcgcaatcacagggaggagcattggcaaggccaggatttttaagcccacaaaggCAAAACCCCTATCTGCGCATGCAGGTAg > 1:332409‑M1/123‑150 (MQ=255) aggcaggcatgatagccatgggctgttggcttgcataggaaaagggcctgcaggagaggccgcgtatcgcaatcacagggaggagcattggcaaggccaggatttttaagcccacaaagaCAAAACCCCTATCTGCGCATGCAGGTAggg > 1:159458‑M1/121‑150 (MQ=255) aaCCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGc > 1:379608/1‑150 (MQ=32) ccTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTaaa > 1:482877/1‑150 (MQ=37) cTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCACAATCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTaaaa > 1:346447/1‑150 (MQ=255) tGATCCCATCCCGATCTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTaaaaa < 2:253563/150‑1 (MQ=38) tcccatcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATcc < 1:254883/150‑1 (MQ=40) ccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCtt < 2:307969/150‑1 (MQ=255) tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGctc < 2:308535/150‑1 (MQ=255) cGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCtgt < 2:379608/150‑1 (MQ=255) aTGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGa < 2:482877/150‑1 (MQ=255) tGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGAt > 2:105145/1‑150 (MQ=255) ccGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATa > 1:225651/1‑150 (MQ=255) ggTAGTGTGGGGTTTCCCCATGTGAGAGTAAGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCgg < 2:408823/150‑1 (MQ=255) gggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAAt < 2:239675/114‑1 (MQ=255) gggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAAt > 1:239675/1‑114 (MQ=255) tttCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAg > 2:148855/1‑150 (MQ=255) tgtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCtttt < 2:346447/150‑1 (MQ=255) gtgAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTc < 1:105145/150‑1 (MQ=255) gagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTAt < 2:341253/150‑1 (MQ=255) gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTg > 1:117684/1‑109 (MQ=255) gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTg < 2:117684/109‑1 (MQ=255) gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTcc < 2:372756/141‑1 (MQ=255) gTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTcc > 1:372756/1‑141 (MQ=255) aggatttttaagcccacaaagaCAAAACCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTc > 1:512191‑M1/23‑140 (MQ=255) aggatttttaagcccacaaagaCAAAACCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTc < 2:512191‑M1/118‑1 (MQ=255) tCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGc > 1:577710/1‑150 (MQ=255) tCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCACACAGc > 2:95925/1‑150 (MQ=255) aGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGCCCa > 2:502787/1‑150 (MQ=255) gATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGCCCAt < 2:577710/150‑1 (MQ=255) ttCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGCCCATgg < 1:148855/150‑1 (MQ=255) ccacaaagaCAAAACCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCt > 2:162320‑M1/10‑114 (MQ=255) ccacaaagaCAAAACCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCt < 1:162320‑M1/105‑1 (MQ=255) tttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTAt < 1:59986/127‑1 (MQ=255) tttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTAt > 2:59986/1‑127 (MQ=255) | GACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGCCCATGG > NC_002947/529792‑530054 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |