Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 741,626 T→C S52G (AGC→GGC)  ybgH ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913741,6260TC99.5% 692.5 / ‑4.8 216S52G (AGC→GGC) ybgHputative transport protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (120/95);  total (120/96)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

Only 100 of 219 total aligned reads displayed.

CCAGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAAT  >  NC_000913/741536‑741715
                                                                                          |                                                                                         
ccAGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCaa                                                                                         >  2:65200/1‑93 (MQ=255)
 cAGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCaaa                                                                                        <  1:4495841/93‑1 (MQ=255)
 cAGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCaaa                                                                                        >  1:6156275/1‑93 (MQ=255)
  aGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAAt                                                                                       <  2:4781335/93‑1 (MQ=255)
  aGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAAt                                                                                       >  2:2942227/1‑93 (MQ=255)
   gCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATa                                                                                      <  1:5729058/93‑1 (MQ=255)
     aTCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAAc                                                                                    <  1:6044619/93‑1 (MQ=255)
      tCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACt                                                                                   >  2:602593/1‑93 (MQ=255)
       cACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTc                                                                                  <  1:5879106/93‑1 (MQ=255)
        aCCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCg                                                                                 <  1:4656439/93‑1 (MQ=255)
        aCCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCg                                                                                 >  2:2148398/1‑93 (MQ=255)
         ccGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGt                                                                                >  1:2724568/1‑93 (MQ=255)
          cGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTa                                                                               <  1:89083/93‑1 (MQ=255)
            ccATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAgg                                                                             >  1:3215122/1‑93 (MQ=255)
             cATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGc                                                                            <  2:876543/93‑1 (MQ=255)
             cATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGc                                                                            >  2:1540774/1‑93 (MQ=255)
              aTGCGATTGCCGAGAACTTTATCCGCCAAAAGGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCg                                                                           <  2:3465079/93‑1 (MQ=255)
              aTGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCg                                                                           <  2:1742130/93‑1 (MQ=255)
              aTGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCg                                                                           <  1:4813158/93‑1 (MQ=255)
              aTGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCg                                                                           <  1:3571763/93‑1 (MQ=255)
               tGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGt                                                                          <  1:2263653/93‑1 (MQ=255)
                 cGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGa                                                                        <  1:5845737/93‑1 (MQ=255)
                  gATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAg                                                                       >  2:1135954/1‑93 (MQ=255)
                  gATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAg                                                                       >  2:742685/1‑93 (MQ=255)
                  gATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAg                                                                       >  1:1667548/1‑93 (MQ=255)
                    ttGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTa                                                                     >  1:4046106/1‑93 (MQ=255)
                    ttGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTa                                                                     >  2:233599/1‑93 (MQ=255)
                    ttGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTa                                                                     >  1:2223591/1‑93 (MQ=255)
                     tGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGCAt                                                                    >  1:1944812/1‑93 (MQ=255)
                      gCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATc                                                                   >  2:1402538/1‑93 (MQ=255)
                        cGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGt                                                                 >  2:619840/1‑93 (MQ=255)
                        cGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGt                                                                 >  1:1090256/1‑93 (MQ=255)
                        cGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGt                                                                 <  1:4107048/93‑1 (MQ=255)
                        cGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGt                                                                 <  1:32664/93‑1 (MQ=255)
                         gagaACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGtt                                                                <  1:4154078/93‑1 (MQ=255)
                         gagaACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGtt                                                                >  2:5579607/1‑93 (MQ=255)
                         gagaACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGtt                                                                <  1:1888451/93‑1 (MQ=255)
                          agaACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTCGGCGTGAGTATCGTTg                                                               >  1:408781/1‑93 (MQ=255)
                           gaACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGt                                                              >  1:2250942/1‑93 (MQ=255)
                           gaACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGt                                                              <  1:1754602/93‑1 (MQ=255)
                            aaCTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTa                                                             <  1:3872317/93‑1 (MQ=255)
                               tttATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTAttt                                                          >  1:159545/1‑93 (MQ=255)
                               tttATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTAttt                                                          >  1:5034031/1‑93 (MQ=255)
                               tttATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTAttt                                                          >  2:5557807/1‑93 (MQ=255)
                                ttATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTAtttt                                                         >  2:3606194/1‑93 (MQ=255)
                                ttATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTAtttt                                                         >  2:2898760/1‑93 (MQ=255)
                                 tATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTa                                                        <  2:445520/93‑1 (MQ=255)
                                 tATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTa                                                        <  2:2687826/93‑1 (MQ=255)
                                   tCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGt                                                      >  1:370385/1‑93 (MQ=255)
                                   tCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGt                                                      >  1:4686556/1‑93 (MQ=255)
                                     cGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTg                                                    <  1:2216331/93‑1 (MQ=255)
                                       ccAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGAt                                                  >  1:303800/1‑93 (MQ=255)
                                        cAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGGAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGAtt                                                 <  1:3680609/93‑1 (MQ=255)
                                          aaaaGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTgg                                               >  2:4924893/1‑93 (MQ=255)
                                           aaaGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTggg                                              >  2:5907626/1‑92 (MQ=255)
                                           aaaGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGt                                              <  2:5467758/93‑1 (MQ=255)
                                            aaGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTg                                             <  1:2453281/93‑1 (MQ=255)
                                            aaGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTg                                             >  2:1968496/1‑93 (MQ=255)
                                              gCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTgag                                           >  2:61413/1‑93 (MQ=255)
                                               ccACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTgaga                                          <  1:5896461/93‑1 (MQ=255)
                                                 aCCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATa                                        >  1:5403244/1‑93 (MQ=255)
                                                 aCCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATa                                        <  1:2024558/93‑1 (MQ=255)
                                                  ccGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATag                                       >  1:1635482/1‑93 (MQ=255)
                                                  ccGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATag                                       >  1:345534/1‑93 (MQ=255)
                                                   cGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATaga                                      <  2:908086/93‑1 (MQ=255)
                                                   cGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATaga                                      >  2:3379294/1‑93 (MQ=255)
                                                   cGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATaga                                      >  2:2056065/1‑93 (MQ=255)
                                                    gAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATagag                                     >  2:2662913/1‑93 (MQ=255)
                                                     aGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAgaga                                    >  2:5750919/1‑93 (MQ=255)
                                                      ggATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGaa                                   <  2:1496137/93‑1 (MQ=255)
                                                       gATTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAAt                                  <  1:4804126/93‑1 (MQ=255)
                                                        aTTGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATc                                 <  2:328114/93‑1 (MQ=255)
                                                         ttGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATca                                <  1:689500/93‑1 (MQ=255)
                                                         ttGGCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATca                                >  2:1317414/1‑93 (MQ=255)
                                                           ggCGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGCGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagc                              <  1:3391515/93‑1 (MQ=255)
                                                             cGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagcag                            <  1:750213/93‑1 (MQ=255)
                                                             cGTGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATcagcag                            <  1:2555204/93‑1 (MQ=255)
                                                               tGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAggg                          >  1:290119/1‑93 (MQ=255)
                                                               tGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAggg                          <  2:2446368/93‑1 (MQ=255)
                                                               tGACATACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAggg                          >  2:3630070/1‑93 (MQ=255)
                                                                   aTACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACg                      <  2:820765/93‑1 (MQ=255)
                                                                    tACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGc                     >  2:4853699/1‑93 (MQ=255)
                                                                    tACACCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGc                     >  1:2666229/1‑93 (MQ=255)
                                                                      cacCAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCAt                   >  1:1214107/1‑93 (MQ=255)
                                                                         cAGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGcc                >  1:580613/1‑93 (MQ=255)
                                                                          aGCGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCg               >  2:5994006/1‑93 (MQ=255)
                                                                            cGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTa             <  1:952455/93‑1 (MQ=255)
                                                                            cGAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTa             >  1:5271863/1‑93 (MQ=255)
                                                                             gAACAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaa            >  1:4024856/1‑93 (MQ=255)
                                                                              aaCAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaaa           <  2:511466/93‑1 (MQ=255)
                                                                              aaCAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaaa           <  1:2075885/93‑1 (MQ=255)
                                                                               aCAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaaaa          <  2:3052667/93‑1 (MQ=255)
                                                                               aCAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaaaa          <  1:1017128/93‑1 (MQ=255)
                                                                               aCAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaaaa          >  1:2786378/1‑93 (MQ=255)
                                                                                cAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAg         <  1:3768499/93‑1 (MQ=255)
                                                                                cAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAg         >  1:5302441/1‑93 (MQ=255)
                                                                                cAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAg         <  1:2380625/93‑1 (MQ=255)
                                                                                cAGTAGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAg         >  1:5483454/1‑93 (MQ=255)
                                                                                    aGGCGCCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaa     <  2:5612235/93‑1 (MQ=255)
                                                                                      gcgcCAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTaaaa   <  2:3758649/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCAGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAAT  >  NC_000913/741536‑741715

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