Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,234,507 | A→T | 18.2% | intergenic (‑52/+232) | higB ← / ← rlmG | mRNA interferase toxin of the HigB‑HigA toxin‑antitoxin system/23S rRNA m(2)G1835 methyltransferase, SAM‑dependent |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,234,507 | 0 | A | T | 18.2% | 77.5 / 13.9 | 33 | intergenic (‑52/+232) | higB/rlmG | mRNA interferase toxin of the HigB‑HigA toxin‑antitoxin system/23S rRNA m(2)G1835 methyltransferase, SAM‑dependent |
Reads supporting (aligned to +/- strand): ref base A (13/14); new base T (3/3); total (16/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 = 1.03e-01 |
TAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTT > NC_000913/3234362‑3234602 | tAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAtt < 2:938415/148‑1 (MQ=255) tCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAt < 1:1014914/141‑1 (MQ=255) tCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAt > 2:1014914/1‑141 (MQ=255) cAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATCTGGCaa > 4:310890/1‑137 (MQ=255) cAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATCTGGCaa < 2:163215/137‑1 (MQ=255) accaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCtgtg > 2:996540/1‑148 (MQ=255) ccGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGga > 4:14901/1‑148 (MQ=255) tttATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGgaagaa < 3:14901/148‑1 (MQ=255) gTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATCTGGCAACgg > 1:186454/1‑117 (MQ=255) gTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATCTGGCAACgg < 2:186454/117‑1 (MQ=255) ttttCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTg > 3:682687/1‑99 (MQ=255) ttttCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGaa < 4:645117/111‑1 (MQ=255) ttttCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGaa > 3:645117/1‑111 (MQ=255) ttttCCGCAGCATCTATCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTg < 4:682687/99‑1 (MQ=255) tCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAgaga < 1:996540/148‑1 (MQ=255) aGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAtt > 3:700239/1‑69 (MQ=255) aGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAACTGTCAATATGATAACATCCGTTGCCAGAtt < 4:700239/69‑1 (MQ=255) ttATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTa > 2:376930/1‑84 (MQ=255) ttATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTa < 1:376930/84‑1 (MQ=255) ggTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGa < 1:780154/59‑1 (MQ=255) ggTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGa > 2:780154/1‑59 (MQ=255) gCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAt > 4:688185/1‑57 (MQ=255) gCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAt < 3:688185/57‑1 (MQ=255) aTTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGaa < 2:208542/95‑1 (MQ=255) aTTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGaa > 1:208542/1‑95 (MQ=255) aaCATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGAtt < 4:464480/104‑1 (MQ=255) aaCATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGAtt > 3:464480/1‑104 (MQ=255) ccGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGAt > 3:612050/1‑86 (MQ=255) ccGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGAt < 4:612050/86‑1 (MQ=255) ccGTTGCCAATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGa < 4:282548/82‑1 (MQ=39) ccGTTGCCAATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGa > 3:282548/1‑82 (MQ=38) aGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGAtttttt > 2:622361/1‑95 (MQ=255) aGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGAtttttt < 1:622361/95‑1 (MQ=255) gATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGtt > 4:101345/1‑97 (MQ=255) gATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGtt < 3:101345/97‑1 (MQ=255) | TAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGATTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTT > NC_000913/3234362‑3234602 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |