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_0009133,234,5070AT18.2% 77.5 / 13.9 33intergenic (‑52/+232)higB/rlmGmRNA 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: