Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 899,657 G→A 10.0% intergenic (+12/+187) rlmC → / ← artJ 23S rRNA m(5)U747 methyltransferase, SAM‑dependent/arginine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913899,6570GA10.0% 147.0 / 6.1 50intergenic (+12/+187)rlmC/artJ23S rRNA m(5)U747 methyltransferase, SAM‑dependent/arginine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (20/25);  new base A (3/2);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TCCAGCTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTAC  >  NC_000913/899511‑899801
                                                                                                                                                  |                                                                                                                                                
tCCAGCTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCAtgtg                                                                                                                                                >  2:50819/1‑149 (MQ=255)
tCCAGCTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCAtgtg                                                                                                                                                >  1:465017/1‑149 (MQ=255)
 ccAGCTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGc                                                                                                                                               <  2:586498/149‑1 (MQ=255)
  cAGCTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCg                                                                                                                                              >  3:133897/1‑149 (MQ=255)
  cAGCTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCg                                                                                                                                              >  1:501654/1‑149 (MQ=255)
     cTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCt                                                                                                                                           >  1:203389/1‑149 (MQ=255)
       gTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCACATGCGGCttt                                                                                                                                         <  1:442698/149‑2 (MQ=255)
        tAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCa                                                                                                                                        <  2:30736/149‑1 (MQ=255)
        tAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCa                                                                                                                                        <  1:571037/149‑1 (MQ=255)
              ccAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGc                                                                                                                                  >  3:59646/1‑149 (MQ=255)
                 aaCTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGCTTGCTGa                                                                                                                               <  2:89142/149‑1 (MQ=255)
                 aaCTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGa                                                                                                                               <  1:221999/149‑1 (MQ=255)
                  aCTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGAc                                                                                                                              <  2:13194/149‑1 (MQ=255)
                       ggCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGt                                                                                                                         <  4:576891/149‑1 (MQ=255)
                       ggCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGt                                                                                                                         <  1:144738/149‑1 (MQ=255)
                                tatCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCgttgt                                                                                                                <  1:514457/149‑1 (MQ=255)
                                     gcgAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATg                                                                                                           >  1:438288/1‑149 (MQ=255)
                                      cgAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGc                                                                                                          <  2:203389/149‑1 (MQ=255)
                                          cTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGa                                                                                                      >  4:337213/1‑149 (MQ=255)
                                           tGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGAt                                                                                                     <  2:50713/149‑1 (MQ=255)
                                              cTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCgttgtt                                                                                                               <  1:259610/136‑1 (MQ=255)
                                              cTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCgttgtt                                                                                                               >  2:259610/1‑136 (MQ=255)
                                                gggTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATgcggc                                                                                                <  4:22383/149‑1 (MQ=255)
                                                   tttCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTa                                                                                             >  2:364782/1‑149 (MQ=255)
                                                       gTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACg                                                                                         >  1:267321/1‑149 (MQ=255)
                                                       gTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACg                                                                                         >  2:447982/1‑149 (MQ=255)
                                                        tATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGc                                                                                        <  1:23543/149‑1 (MQ=255)
                                                          ttGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCt                                                                                      >  4:378045/1‑149 (MQ=255)
                                                          ttGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCt                                                                                      <  1:447982/149‑1 (MQ=255)
                                                                          tttCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                                      <  2:267321/149‑1 (MQ=255)
                                                                              gATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGt                                                                  >  1:460433/1‑149 (MQ=255)
                                                                                   gtgcccGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCaaa                                                             >  3:388768/4‑149 (MQ=255)
                                                                                    ttCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAAt                                                            >  2:117033/1‑149 (MQ=255)
                                                                                        cGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCaa                                                        <  4:459889/149‑1 (MQ=255)
                                                                                        cGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCaa                                                        <  2:465017/149‑1 (MQ=255)
                                                                                               cgcgCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTg                                                 >  2:143157/1‑149 (MQ=255)
                                                                                                    aCTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATgg                                            <  1:364782/149‑1 (MQ=255)
                                                                                                          aaGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATg                                      >  3:368587/1‑149 (MQ=255)
                                                                                                            gTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTa                                    <  3:337213/149‑1 (MQ=255)
                                                                                                              gCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAgg                                  <  4:388768/149‑1 (MQ=255)
                                                                                                                 gACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCt                               <  3:324266/149‑1 (MQ=255)
                                                                                                                   cGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGa                             <  2:501654/149‑1 (MQ=255)
                                                                                                                        ctgGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGc                        <  2:460433/149‑1 (MQ=255)
                                                                                                                           gTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTa                     >  2:269841/1‑149 (MQ=255)
                                                                                                                              aaGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAgcg                  >  2:543620/1‑149 (MQ=255)
                                                                                                                                   aTAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCa             >  2:296176/1‑149 (MQ=255)
                                                                                                                                     tctgaaGCCGCACATGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAgg           >  4:185052/5‑149 (MQ=255)
                                                                                                                                      aaaaaGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGc          >  2:158373/1‑149 (MQ=255)
                                                                                                                                        gaaGCCGCACATGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa        >  4:15567/2‑149 (MQ=255)
                                                                                                                                           gCCGCACATGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAg                    <  2:409641/134‑1 (MQ=255)
                                                                                                                                           gCCGCACATGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAg                    >  1:409641/1‑134 (MQ=255)
                                                                                                                                              gCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTAc  <  1:296176/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TCCAGCTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTAC  >  NC_000913/899511‑899801

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: