Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 899,656 T→C 10.0% intergenic (+11/+188) 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,6560TC10.0% 143.9 / 6.8 50intergenic (+11/+188)rlmC/artJ23S rRNA m(5)U747 methyltransferase, SAM‑dependent/arginine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (20/25);  new base C (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.61e-01

TACTCCAGCTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTAC  >  NC_000913/899508‑899801
                                                                                                                                                    |                                                                                                                                                 
tACTCCAGCTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCAt                                                                                                                                                   <  1:437087/149‑1 (MQ=255)
   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                    >  1:409641/1‑134 (MQ=255)
                                                                                                                                              gCCGCACATGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAg                    <  2:409641/134‑1 (MQ=255)
                                                                                                                                                 gCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTAc  <  1:296176/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TACTCCAGCTGTAACGCCCAAACTATGGCGAAAGATATCCGCGAACTGCCTGGGTTTCGTATTGAACGGGTGCAGCTTTTCGATATGTTCCCGCATACCGCGCACTATGAAGTGCTGACGCTGCTGGTGAAGCAATAAAAAAGCCGCATGTGCGGCTTCAGATTGCTGACAAAGTGCGCGTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTAC  >  NC_000913/899508‑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: