Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,806,842 C→T R226C (CGC→TGC)  kdtA → 3‑deoxy‑D‑manno‑octulosonic‑acid transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,806,8420CT97.0% 577.0 / 1.4 166R226C (CGC→TGC) kdtA3‑deoxy‑D‑manno‑octulosonic‑acid transferase
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base T (67/94);  total (71/95)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

Only 100 of 175 total aligned reads displayed.

GCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTG  >  NC_000913/3806752‑3806931
                                                                                          |                                                                                         
gcAAAAAATAATCAGGTGACCGTTACCGATAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACgctgcgctgc                                                                                         <  1:4631239/93‑1 (MQ=255)
  aaaaaaTAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCa                                                                                       <  2:3284139/93‑1 (MQ=255)
     aaaTAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTg                                                                                    <  1:120916/93‑1 (MQ=255)
     aaaTAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTg                                                                                    <  1:682988/93‑1 (MQ=255)
     aaaTAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTg                                                                                    <  2:3939421/93‑1 (MQ=255)
     aaaTAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTg                                                                                    <  2:1741125/93‑1 (MQ=255)
      aataatCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTgg                                                                                   <  1:446639/93‑1 (MQ=255)
      aataatCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTgg                                                                                   <  2:2622518/93‑1 (MQ=255)
       ataatCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTggg                                                                                  >  1:987368/1‑93 (MQ=255)
        taatCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGc                                                                                 >  2:1695179/1‑93 (MQ=255)
          atCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGcac                                                                               <  1:3892275/93‑1 (MQ=255)
           tCAGATGACCGTTGCCGGTAGCCTGAAATTCGATATGTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGcacc                                                                              <  2:2191079/93‑1 (MQ=255)
            cAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGcacca                                                                             >  1:3757549/1‑93 (MQ=255)
            cAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGcacca                                                                             <  1:3993140/93‑1 (MQ=255)
              ggTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACcaca                                                                           <  2:2658934/93‑1 (MQ=255)
               gTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCacac                                                                          <  2:605043/93‑1 (MQ=255)
                tGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACAcc                                                                         <  1:3571816/93‑1 (MQ=255)
                  aCCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGc                                                                       >  2:3810008/1‑93 (MQ=255)
                      ttACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCgg                                                                   >  1:1794564/1‑93 (MQ=255)
                      ttACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCgg                                                                   >  1:1586268/1‑93 (MQ=255)
                      ttACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCgg                                                                   <  2:3971350/93‑1 (MQ=255)
                        aCCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTGACGCCCGGACGCTGCGGGGCCAGTAGGCACGCCCCCGCCCggat                                                                 >  1:2802004/1‑91 (MQ=255)
                          cGGTAGCCTGAAATTTGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTCCGCTGCCAGTGGGCACCACACCGCCCGGTATg                                                               <  2:4084868/93‑1 (MQ=255)
                          cGGTAGCCTGAAATTCGATATTTCTGTACCGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGCGGGCACCACACCGCCCGGTATg                                                               <  2:2992593/93‑1 (MQ=255)
                          cGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATg                                                               >  2:252752/1‑93 (MQ=255)
                          cGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATg                                                               <  1:128464/93‑1 (MQ=255)
                           ggTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATgg                                                              >  2:4574770/1‑93 (MQ=255)
                           ggTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATgg                                                              >  1:4244560/1‑93 (MQ=255)
                            gTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGa                                                             >  1:259981/1‑93 (MQ=255)
                               gCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTg                                                          >  1:4678969/1‑93 (MQ=255)
                                 cTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGcc                                                        >  2:1941958/1‑93 (MQ=255)
                                 cTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGcc                                                        >  1:1119401/1‑93 (MQ=255)
                                  tGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGcca                                                       >  2:3660077/1‑93 (MQ=255)
                                   gAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGccac                                                      >  1:3635338/1‑93 (MQ=255)
                                   gAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGccac                                                      >  2:1075308/1‑93 (MQ=255)
                                     aaTTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGccacca                                                    <  1:1783992/93‑1 (MQ=255)
                                      aTTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAg                                                   <  1:1025767/93‑1 (MQ=255)
                                       ttCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGc                                                  >  2:432657/1‑93 (MQ=255)
                                       ttCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCAGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGc                                                  >  1:1917876/1‑93 (MQ=255)
                                         cGATATTTCTGTAACGCCGCAGTTTGCTGCTAAAGCCGTCACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCAc                                                <  2:4038195/93‑1 (MQ=255)
                                         cGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCAc                                                <  2:3763357/93‑1 (MQ=255)
                                          gATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACt                                               >  2:3496579/1‑93 (MQ=255)
                                          gATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCAACAGCACt                                               >  1:2663612/1‑93 (MQ=255)
                                           atatTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTc                                              <  2:3885/93‑1 (MQ=255)
                                            tatTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCa                                             <  1:4385334/93‑1 (MQ=255)
                                            tatTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCa                                             >  1:169281/1‑93 (MQ=255)
                                             atTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCAc                                            <  2:220501/93‑1 (MQ=255)
                                             atTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCAc                                            <  1:719670/93‑1 (MQ=255)
                                             atTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCAc                                            <  2:4137133/93‑1 (MQ=255)
                                             atTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCAc                                            <  1:415938/93‑1 (MQ=255)
                                             atTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCAc                                            >  2:736803/1‑93 (MQ=255)
                                              tttCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACg                                           <  2:150411/93‑1 (MQ=255)
                                                tCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGaa                                         >  1:3639145/1‑93 (MQ=255)
                                                 cTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAg                                        <  1:4536362/93‑1 (MQ=255)
                                                  tGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAgg                                       <  2:4178320/93‑1 (MQ=255)
                                                     aaCGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGa                                    >  2:2877557/1‑93 (MQ=255)
                                                     aaCGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGa                                    >  2:961686/1‑93 (MQ=255)
                                                     aaCGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGa                                    >  1:2893222/1‑93 (MQ=255)
                                                      aCGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGaa                                   >  2:3659568/1‑93 (MQ=255)
                                                       cgccgcAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGaaa                                  >  2:3577726/1‑92 (MQ=255)
                                                        gccgcAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGcgacga                                 >  1:3723183/1‑93 (MQ=255)
                                                         ccgcAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAagag                                <  1:2768397/93‑1 (MQ=255)
                                                         ccgcAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAagag                                <  2:3096364/93‑1 (MQ=255)
                                                         ccgcAGTTGACTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACGGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAagag                                <  2:2816472/93‑1 (MQ=255)
                                                          cgcAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAagaga                               >  2:3044961/1‑93 (MQ=255)
                                                           gcAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAagagag                              <  2:2808317/93‑1 (MQ=255)
                                                             aGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAgtg                            <  2:2564631/93‑1 (MQ=255)
                                                             aGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAgtg                            >  1:798332/1‑93 (MQ=255)
                                                              gTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAgagt                           >  1:2372711/1‑93 (MQ=255)
                                                               ttGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGtgtg                          >  1:1449181/1‑93 (MQ=255)
                                                               ttGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGtgtg                          <  2:1632894/93‑1 (MQ=255)
                                                               ttGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGtgtg                          <  1:3793033/93‑1 (MQ=255)
                                                               ttGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGtgtg                          >  1:583925/1‑93 (MQ=255)
                                                                tGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTgtgg                         <  1:2088625/93‑1 (MQ=255)
                                                                 ggCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTgtggt                        >  1:2149684/1‑93 (MQ=255)
                                                                 ggCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTgtggt                        <  1:2958626/93‑1 (MQ=255)
                                                                  gctgctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTgtggtg                       >  1:1373290/1‑93 (MQ=255)
                                                                   ctgctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGa                      <  2:2686529/93‑1 (MQ=255)
                                                                    tgctAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGCAGGGGAAGCGAGTgtggtgct                     >  1:1634968/1‑91 (MQ=255)
                                                                     gctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATc                    <  2:686772/93‑1 (MQ=255)
                                                                     gctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATc                    <  2:1279634/93‑1 (MQ=255)
                                                                      ctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGaccc                   >  1:2170168/1‑90 (MQ=255)
                                                                      ctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcg                   <  2:3126222/93‑1 (MQ=255)
                                                                      ctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcg                   >  1:1709476/1‑93 (MQ=255)
                                                                       tAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTTCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcgc                  >  2:2594885/1‑93 (MQ=255)
                                                                        aaaGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCAACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcgcc                 <  1:4091019/93‑1 (MQ=255)
                                                                          aGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAGGAGAGGGGGGTGATggccgc               >  1:1562486/1‑93 (MQ=255)
                                                                          aGCCGTGACGCAGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcgccgc               <  1:2653872/93‑1 (MQ=255)
                                                                           gCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAACAGAGGGTGGTGATCGCCGca              >  1:4450823/1‑93 (MQ=255)
                                                                            ccGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGcac             <  2:4272134/93‑1 (MQ=255)
                                                                             cGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGcaca            <  2:3882250/93‑1 (MQ=255)
                                                                             cGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGcaca            <  2:2389107/93‑1 (MQ=255)
                                                                              gcgACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCCCCAGCACTCCAGAAGGACACGAGAGCGTGGTGCTCGTCGCACAt           >  1:4569111/3‑93 (MQ=255)
                                                                              gTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACAt           >  1:3302473/1‑93 (MQ=255)
                                                                               tGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATc          <  1:3376617/93‑1 (MQ=255)
                                                                               tGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATc          <  2:1357597/93‑1 (MQ=255)
                                                                                gACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCa         <  1:1360273/93‑1 (MQ=255)
                                                                                 aCGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAg        <  2:170026/93‑1 (MQ=255)
                                                                                     tGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCAt    <  2:573064/93‑1 (MQ=255)
                                                                                      gcgcTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCAtt   >  2:3199480/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTG  >  NC_000913/3806752‑3806931

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