Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,469 +G intergenic (+36/‑197) b4555 → / → yigM hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,4661.G96.5% 229.0 / 2.5 87intergenic (+33/‑200)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base G (51/33);  total (53/34)
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.00e+00

GGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATT  >  NC_000913/4008403‑4008533
                                                                     |                                                                        
ggCTACAGCCTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑gaaag                                                                >  2:915022/1‑74 (MQ=255)
        ccTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCAtt                                                        <  1:1872596/75‑1 (MQ=255)
         cTGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTa                                                       >  1:3886085/1‑76 (MQ=255)
          tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTAt                                                      >  1:4499619/1‑76 (MQ=255)
          tGAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTAt                                                      <  1:2935877/76‑1 (MQ=255)
           gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATg                                                     >  1:1245432/1‑76 (MQ=255)
           gAAGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATg                                                     >  1:1065087/1‑76 (MQ=255)
            aaGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATgg                                                    <  1:3820759/76‑1 (MQ=255)
            aaGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATgg                                                    <  1:4232602/76‑1 (MQ=255)
             aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATggg                                                   <  1:2548493/76‑1 (MQ=255)
             aGGATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATggg                                                   >  2:3253005/1‑76 (MQ=255)
             aGGATTGTATTGCGTTTGGCGACGGGCTGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTCTggg                                                   >  1:5543071/1‑76 (MQ=255)
              ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGc                                                  >  1:4593855/1‑76 (MQ=255)
              ggATTGTATTGCGTTTGGTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGc                                                  >  1:3112/1‑76 (MQ=255)
                            ttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCg                                    <  2:1217182/73‑1 (MQ=255)
                            ttggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCg                                    >  1:3697610/4‑76 (MQ=255)
                             tggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGt                                   >  1:2569580/3‑76 (MQ=255)
                              ggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                  >  1:5050634/2‑76 (MQ=255)
                              ggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                  >  2:3992761/2‑76 (MQ=255)
                              ggtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                  <  1:3622319/75‑1 (MQ=255)
                              ggggACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTc                                  <  2:5529668/73‑1 (MQ=255)
                               gtgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCt                                 <  1:4076250/76‑1 (MQ=255)
                                tgACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTg                                >  2:4333679/1‑76 (MQ=255)
                                 gaggggAGGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGa                               <  2:5394620/73‑1 (MQ=255)
                                 gaggggAGGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCAGTATGGGCAGTGCGCACCAGCGTCTGa                               <  2:5449293/73‑1 (MQ=255)
                                  aCGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACTAGCGTCTGaa                              >  1:2869631/1‑76 (MQ=255)
                                  aCGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGaa                              <  1:3406146/76‑1 (MQ=255)
                                   cGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGaaa                             >  1:1231634/1‑76 (MQ=255)
                                   cGGGATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGaaa                             <  1:4096738/76‑1 (MQ=255)
                                    gggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCgcacca                                       >  2:5655796/1‑65 (MQ=255)
                                    gggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAg                            >  1:2080745/1‑76 (MQ=255)
                                    gggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAg                            >  2:2069747/1‑76 (MQ=255)
                                    gggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAg                            >  2:1092253/1‑76 (MQ=255)
                                    gggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATAATGGGCAGTGCGCACCAGCGACTGAAAg                            >  2:3901876/1‑76 (MQ=255)
                                     ggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGa                           >  2:1874076/1‑76 (MQ=255)
                                     ggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGa                           <  2:5842576/76‑1 (MQ=255)
                                     ggATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGa                           >  2:3855996/1‑76 (MQ=255)
                                      gATGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAACGCTGCATTATGGGCAGCGCGCACCAGCGTCTGAAAGAc                          >  2:4254916/1‑76 (MQ=255)
                                       aTGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGAcc                         >  1:1331522/1‑76 (MQ=255)
                                       aTGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGAcc                         >  2:1230256/1‑76 (MQ=255)
                                        tGAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCt                        >  1:5779737/1‑76 (MQ=255)
                                        gggacgacgCCGAAATGCTGGCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCt                        <  2:5428302/73‑1 (MQ=255)
                                         gAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCtt                       >  1:4304521/1‑76 (MQ=255)
                                         gAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCtt                       >  1:3437078/1‑76 (MQ=255)
                                         gAACGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCtt                       >  2:3564124/1‑76 (MQ=255)
                                         gAACGACGCCGAAATGCTGTCGATGGCGGGGCAAGGCTGCATTATGGGCCGTGCGCCCCAGCGTCTGCAAGAcata                       >  2:3597522/1‑73 (MQ=255)
                                          aaCGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTc                      >  1:4959963/1‑76 (MQ=255)
                                          aaCGACGCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTc                      >  1:3855544/1‑76 (MQ=255)
                                           acgactCCTAAATGCGGTCGAGGGCGGGG‑‑AAAGGCTTCATTATGGCCAGTGCGCACCAGCGTCTGAAAGACCTTCa                     <  2:1062961/76‑1 (MQ=255)
                                           acgacgCCGAAATGCTGTCGATGNGGGGG‑‑NAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCa                     <  2:1046323/76‑1 (MQ=255)
                                            cgacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAt                    <  1:972753/76‑1 (MQ=255)
                                            cgacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAt                    >  2:2730353/1‑76 (MQ=255)
                                            cgacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAt                    >  2:1530832/1‑76 (MQ=255)
                                            cgacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAt                    >  1:5670577/1‑76 (MQ=255)
                                             gacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATc                   >  1:1664774/1‑76 (MQ=255)
                                             gacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATc                   <  2:4309193/76‑1 (MQ=255)
                                             gacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCCTCATc                   >  2:4835137/1‑76 (MQ=255)
                                             gacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGGCTGAAAGACCTTCATc                   >  1:3541422/1‑76 (MQ=255)
                                              acgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCAtac                  >  2:3453286/1‑74 (MQ=255)
                                              acgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                  >  2:3083447/1‑76 (MQ=255)
                                              acgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                  >  1:1802523/1‑76 (MQ=255)
                                               cgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATccc                 <  1:1735251/76‑1 (MQ=255)
                                               cgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATccc                 <  1:5288584/76‑1 (MQ=255)
                                               cgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATccc                 >  1:3237918/1‑76 (MQ=255)
                                               cgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGGGCGCACCAGCGTCTGAAAGACCTTCATccc                 >  1:3252514/1‑76 (MQ=255)
                                                gCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCg                <  2:1684409/76‑1 (MQ=255)
                                                gCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCg                >  1:4952937/1‑76 (MQ=255)
                                                gCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCg                >  1:300196/1‑76 (MQ=255)
                                                 ccGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa               >  1:579879/1‑76 (MQ=255)
                                                 ccGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGa               <  2:4959962/76‑1 (MQ=255)
                                                  cGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAg              >  1:993413/1‑76 (MQ=255)
                                                  cGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAg              <  2:3867434/76‑1 (MQ=255)
                                                    agaTGCTGTCGGTGGCGGGG‑‑AGAGTCTGCAGTAGGGTCAGTGCGCACCAGCGTGTGAAAGACCTTCAGCCCGAGcg            <  1:236779/74‑2 (MQ=255)
                                                    aaaTGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCt            <  1:3796540/76‑1 (MQ=255)
                                                    aaaTGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCt            <  2:4472166/76‑1 (MQ=255)
                                                    aaaTGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCt            <  1:198609/76‑1 (MQ=255)
                                                    aaaTGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCt            <  2:857560/76‑1 (MQ=255)
                                                     aaTGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTg          <  2:2210680/76‑1 (MQ=255)
                                                      aTGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgg         <  1:5734751/76‑2 (MQ=255)
                                                      aTGCTGGCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgg         <  1:5739028/76‑2 (MQ=255)
                                                       tGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTgga        >  1:1095076/1‑74 (MQ=255)
                                                        gCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTggaa       >  1:5112817/1‑73 (MQ=255)
                                                           gTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTg     >  2:4588871/1‑76 (MQ=255)
                                                           gTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTg     <  1:586117/76‑1 (MQ=255)
                                                              gATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt  >  2:3665213/1‑76 (MQ=255)
                                                              gATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGGGAtt  >  2:613244/1‑76 (MQ=255)
                                                              gATGGCGGGG‑‑AAAGTCTGCATTATGGGCAGTGCGCACGAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGAtt  <  1:4254917/76‑1 (MQ=255)
                                                                     |                                                                        
GGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATT  >  NC_000913/4008403‑4008533

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