Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,008,473 +G intergenic (+40/‑193) 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,4721.G96.5% 223.5 / 2.6 87intergenic (+39/‑194)b4555/yigMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (50/33);  total (52/35)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAA  >  NC_000913/4008403‑4008541
                                                                              |                                                                       
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                           <  2:4309193/76‑1 (MQ=255)
                                             gacgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATc                           >  1:1664774/1‑76 (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                          >  1:1802523/1‑76 (MQ=255)
                                              acgCCGAAATGCTGTCGATGGCGGGG‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATcc                          >  2:3083447/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)
                                                                       gg‑‑AAAGGCTGCATTATGGGCAGTGCGCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGGAAGTGATTGGTACTaa  <  2:3582901/76‑1 (MQ=255)
                                                                              |                                                                       
GGCTACAGCTGAAGGATTGATTGCGTT‑‑GTGACGGGATGAACGACGCCGAAATGCTGTCGATGGCGGG‑‑AAAGCTGCATTATGGGCAGTGGCCACCAGCGTCTGAAAGACCTTCATCCCGAGCTGAAGTGATTGGTACTAA  >  NC_000913/4008403‑4008541

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: