Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,039,960 Δ1 bp intergenic (+42/‑36) yihD → / → yihE hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,039,9590A.99.2% 539.8 / ‑3.1 121intergenic (+41/‑37)yihD/yihEhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (61/58);  minor base C (1/0);  total (63/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 127 total aligned reads displayed.

CGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTAC  >  NC_000913/4039884‑4040029
                                                                           |                                                                      
cgcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATa                                                                        >  1:4690182/1‑76 (MQ=255)
cgcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATa                                                                        >  1:226612/1‑76 (MQ=255)
 gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATat                                                                       >  2:5706668/1‑75 (MQ=255)
   tGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatcc                                                                     <  1:1607056/76‑4 (MQ=255)
     ttACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCtt                                                                  <  2:5807381/76‑1 (MQ=255)
     ttACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCtt                                                                  <  2:488557/76‑1 (MQ=255)
     ttACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCtt                                                                  <  1:965471/76‑1 (MQ=255)
      tACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTc                                                                 <  1:188648/76‑1 (MQ=255)
       aCGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCg                                                                <  1:3880066/76‑1 (MQ=255)
         gcgcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTc                                                              >  2:1274896/1‑76 (MQ=255)
            cTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCAtt                                                           >  1:3186978/1‑76 (MQ=255)
            cTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCAtt                                                           >  1:4220892/1‑76 (MQ=255)
            cTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCAtt                                                           <  1:2320630/76‑1 (MQ=255)
            cTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCAtt                                                           <  2:3267190/76‑1 (MQ=255)
             tCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTc                                                          <  2:3496606/76‑1 (MQ=255)
             tCGCGGCCTACTTAAAGAGCAGAAACTAGGATGCGGCGACACCAAAATCATCGTGAAATGATATCCCTCGTCagta                                                          >  2:1950134/1‑73 (MQ=255)
               gcggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGt                                                        <  2:3355494/76‑1 (MQ=255)
               gcggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGt                                                        >  1:5343063/1‑76 (MQ=255)
               gcggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGt                                                        >  2:5639703/1‑76 (MQ=255)
                cggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTa                                                       <  2:5195516/76‑1 (MQ=255)
                 ggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTaa                                                      >  2:1819760/1‑76 (MQ=255)
                  gcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAAt                                                     <  1:5816875/76‑1 (MQ=255)
                  gcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAAt                                                     >  1:1208710/1‑76 (MQ=255)
                   cgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATg                                                    >  1:4519411/1‑76 (MQ=255)
                    gTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGt                                                   <  2:3205321/76‑1 (MQ=255)
                    gTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGt                                                   >  1:710574/1‑76 (MQ=255)
                    gTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGt                                                   >  1:4760904/1‑76 (MQ=255)
                    gTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGt                                                   <  1:4559532/76‑1 (MQ=255)
                     tAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtt                                                  <  1:3613557/76‑1 (MQ=255)
                     tAATTAAAGAGTAAAAGCTGGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtt                                                  <  1:5760679/76‑1 (MQ=255)
                      aaTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGttt                                                 >  1:2975439/1‑76 (MQ=255)
                      aaTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGttt                                                 >  1:1430458/1‑76 (MQ=255)
                      aaTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGttt                                                 <  1:320853/76‑1 (MQ=255)
                       aTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtttt                                                >  2:921096/1‑76 (MQ=255)
                       aTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGtttt                                                <  1:3773735/76‑1 (MQ=255)
                        tgaaaGAGTAAAAGCGTGGAAGCGGAGACACCAAAAGCGTAGTGGAATGAGATCAGTCGGCATTCGTAATGTTGTc                                               >  1:1869028/3‑76 (MQ=255)
                         tAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTcc                                              <  2:862243/76‑1 (MQ=255)
                           aaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCgg                                            >  1:1694245/1‑76 (MQ=255)
                           aaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCgg                                            >  2:2300926/1‑76 (MQ=255)
                           aaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCgg                                            <  1:522719/76‑1 (MQ=255)
                            agagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGga                                           <  2:3626358/76‑1 (MQ=255)
                             gagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgat                                          >  1:2793961/1‑76 (MQ=255)
                             gagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgat                                          <  1:4779884/76‑1 (MQ=255)
                              agTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatg                                         <  2:317727/76‑1 (MQ=255)
                              agTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatg                                         >  1:5137349/1‑76 (MQ=255)
                              agTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatg                                         <  2:5282721/76‑1 (MQ=255)
                               gTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatga                                        <  2:4246766/76‑1 (MQ=255)
                               gTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGgatga                                        >  1:946730/1‑76 (MQ=255)
                                 aaaaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGatgatg                                      >  2:4553097/1‑76 (MQ=255)
                                 aaaaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGatgatg                                      >  2:3543664/1‑76 (MQ=255)
                                   aaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATggg                                    >  1:1213502/1‑76 (MQ=255)
                                    aGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGa                                   <  1:921096/76‑1 (MQ=255)
                                     gCTTGTGAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGAt                                  <  1:3436540/76‑1 (MQ=255)
                                      cTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATg                                 <  1:1476775/76‑1 (MQ=255)
                                       ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGa                                >  1:1875224/1‑76 (MQ=255)
                                       ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCGTCGTCATTCGTAATGTTTTCCGGATGATGGGATGa                                >  1:3063305/1‑76 (MQ=255)
                                         gTAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAAt                              <  2:2108360/76‑1 (MQ=255)
                                          taggcggggCCACCAAAATCAGCGGGAAATGATATCCTTCGTCATTCGTAATGTGTTCCGGATGAGGGGATGAATa                             <  2:3861407/73‑1 (MQ=255)
                                          tAAGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATa                             <  2:3849176/76‑1 (MQ=255)
                                           aaGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATaa                            >  1:3219409/1‑76 (MQ=255)
                                           aaGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATaa                            >  1:5477773/1‑76 (MQ=255)
                                           aaGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATaa                            >  2:454115/1‑76 (MQ=255)
                                           aaGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATaa                            <  2:68258/76‑1 (MQ=255)
                                            aGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAAc                           >  2:5611618/1‑76 (MQ=255)
                                            aGCGGCGCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAAc                           >  2:598052/1‑76 (MQ=255)
                                             gcggcgCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACa                          <  2:5023402/76‑1 (MQ=255)
                                               ggcgCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgc                        <  2:5357158/76‑1 (MQ=255)
                                               ggcgCCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgc                        >  1:4387055/1‑76 (MQ=255)
                                                gcgcGACCAAAAGCATCTTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcg                       <  2:2819882/76‑1 (MQ=255)
                                                gcgcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcg                       >  2:4662580/1‑76 (MQ=255)
                                                 cgcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcgc                      <  2:1868969/76‑1 (MQ=255)
                                                 cgcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcgc                      >  2:3126868/1‑76 (MQ=255)
                                                  gcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCt                     >  2:4663925/1‑76 (MQ=255)
                                                  gcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCt                     >  1:1013157/1‑76 (MQ=255)
                                                  gcCACCAAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCt                     >  2:5404530/1‑76 (MQ=255)
                                                   ccaccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCtt                    >  1:4752309/1‑76 (MQ=255)
                                                   ccaccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCtt                    >  2:385438/1‑76 (MQ=255)
                                                    caccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCttt                   <  1:1740959/76‑1 (MQ=255)
                                                     accaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCtttt                  <  2:3529710/76‑1 (MQ=255)
                                                     accaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCtttt                  >  2:4255905/1‑76 (MQ=255)
                                                     accaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCtttt                  <  1:2633360/76‑1 (MQ=255)
                                                      ccaAAATCATCGTGAAATTATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGCATGAGTAACATCGCTTTTa                 <  1:2531181/76‑1 (MQ=255)
                                                      ccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTa                 >  1:4230896/1‑76 (MQ=255)
                                                      ccaAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTa                 <  2:1065984/76‑1 (MQ=255)
                                                       caAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTAc                >  1:1591820/1‑76 (MQ=255)
                                                       caAAATCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTAc                >  2:289782/1‑76 (MQ=255)
                                                        aaaaTCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACt               >  1:2689453/1‑76 (MQ=255)
                                                         aaaTCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACtt              >  2:1782798/1‑76 (MQ=255)
                                                          aaTCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACttt             >  1:4309599/1‑76 (MQ=255)
                                                          aaTCATCGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACttt             <  1:2310130/76‑1 (MQ=255)
                                                           atcatcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTc            >  1:482073/1‑76 (MQ=255)
                                                            tcatcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTcc           >  1:1335471/1‑76 (MQ=255)
                                                             catcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCa          <  1:1465725/76‑1 (MQ=255)
                                                             catcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCa          >  1:3933789/1‑76 (MQ=255)
                                                              atcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAg         >  1:2802366/1‑76 (MQ=255)
                                                              atcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAg         >  2:3383152/1‑76 (MQ=255)
                                                               tcGTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGa        >  1:1222366/1‑76 (MQ=255)
                                                                 gTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGaca      >  2:398849/1‑76 (MQ=255)
                                                                 gTGAAATGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGaca      >  2:4296857/1‑76 (MQ=255)
                                                                     aaTGATATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTac  <  1:2506224/76‑1 (MQ=255)
                                                                           |                                                                      
CGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTAC  >  NC_000913/4039884‑4040029

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