Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,043,403 A→G H260H (CAT→CAC yihG ← putative endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,043,4030AG100.0% 846.1 / NA 262H260H (CAT→CACyihGputative endonuclease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (137/125);  total (137/125)

Only 100 of 278 total aligned reads displayed.

TGAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAG  >  NC_000913/4043313‑4043492
                                                                                          |                                                                                         
tGAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCACCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTg                                                                                         <  2:1219028/93‑1 (MQ=255)
 gAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAGCTCATCGGCAATAGGCTGTAAATCCACgtgt                                                                                        <  2:2291747/93‑1 (MQ=255)
   gTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTAc                                                                                      >  1:3076099/1‑93 (MQ=255)
   gTTCAACCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGCAc                                                                                      >  2:1138673/1‑93 (MQ=255)
    ttCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACg                                                                                     >  1:4464895/1‑93 (MQ=255)
     tCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGa                                                                                    <  2:3226950/93‑1 (MQ=255)
     tCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGa                                                                                    <  1:3532528/93‑1 (MQ=255)
      cAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGAc                                                                                   <  2:3275817/93‑1 (MQ=255)
       aGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACa                                                                                  >  2:3615503/1‑93 (MQ=255)
        gCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACaa                                                                                 <  2:220070/93‑1 (MQ=255)
        gCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACaa                                                                                 >  1:263969/1‑93 (MQ=255)
         ccATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAAt                                                                                >  1:652521/1‑93 (MQ=255)
          cATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATc                                                                               <  2:1976336/93‑1 (MQ=255)
           aTTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATcg                                                                              >  1:1695298/1‑92 (MQ=255)
           aTTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATcc                                                                              >  2:476075/1‑93 (MQ=255)
           aTTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATcc                                                                              >  2:762764/1‑93 (MQ=255)
            ttGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCcg                                                                             >  1:1866926/1‑93 (MQ=255)
             tGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCcgc                                                                            <  2:1750987/93‑1 (MQ=255)
             tGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCcgc                                                                            <  1:2285118/93‑1 (MQ=255)
             tGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCcgc                                                                            <  2:1051050/93‑1 (MQ=255)
               cTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGt                                                                          <  2:4203813/93‑1 (MQ=255)
               cTGAAAATGGCGCTTAAAACTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGt                                                                          <  2:3628686/93‑1 (MQ=255)
                 gAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCa                                                                        >  2:500376/1‑93 (MQ=255)
                   aaaTGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGt                                                                      >  1:1393968/1‑93 (MQ=255)
                     aTGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGttt                                                                    <  1:4376000/93‑1 (MQ=255)
                      tGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTa                                                                   >  1:2952455/1‑93 (MQ=255)
                           cTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCt                                                              >  1:2208618/1‑93 (MQ=255)
                           cTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCt                                                              >  1:3051998/1‑93 (MQ=255)
                           cTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCt                                                              <  2:3803432/93‑1 (MQ=255)
                            ttAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCtt                                                             >  1:2039476/1‑93 (MQ=255)
                            ttAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCtt                                                             >  1:1150770/1‑93 (MQ=255)
                              aaaGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTaa                                                           <  1:3112771/93‑1 (MQ=255)
                              aaaGCTTTTGTCATTGATGTAGCCGCCATGTAACTCATCTGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTaa                                                           <  2:4391700/93‑1 (MQ=255)
                               aaGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAAc                                                          >  1:1745190/1‑93 (MQ=255)
                               aaGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAAc                                                          >  1:4642979/1‑93 (MQ=255)
                                 gCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACat                                                        <  1:2552960/93‑1 (MQ=255)
                                  cTTTTGTCATTGATGTAGTGGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTCTACCGCTTAACata                                                       >  1:1642752/1‑93 (MQ=255)
                                  cTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACata                                                       >  2:4600600/1‑93 (MQ=255)
                                   ttttGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACatat                                                      <  1:710106/93‑1 (MQ=255)
                                    tttGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATc                                                     <  1:3905420/93‑1 (MQ=255)
                                     ttGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCg                                                    <  2:927202/93‑1 (MQ=255)
                                      tGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCga                                                   >  1:27853/1‑93 (MQ=255)
                                         cATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCgaaga                                                <  2:3954583/93‑1 (MQ=255)
                                         cATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCgaaga                                                <  1:608217/93‑1 (MQ=255)
                                         cATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCgaaca                                                <  2:2446724/93‑3 (MQ=255)
                                          aTTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCgaagaa                                               <  1:4152789/93‑1 (MQ=255)
                                          aTTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCgaagaa                                               >  2:2629581/1‑93 (MQ=255)
                                          aTTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCgaagaa                                               >  2:4433443/1‑93 (MQ=255)
                                          aTTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCgaagaa                                               >  1:243378/1‑93 (MQ=255)
                                           ttGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAAc                                              <  1:1788305/93‑1 (MQ=255)
                                           ttGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAAc                                              >  1:2695730/1‑93 (MQ=255)
                                            tGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACg                                             <  2:3696813/93‑1 (MQ=255)
                                            tGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACg                                             <  2:3449241/93‑1 (MQ=255)
                                             gATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACgg                                            >  1:2973186/1‑93 (MQ=255)
                                             gATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACgg                                            <  2:692608/93‑1 (MQ=255)
                                             gATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAAACCGCGTCAGTTTACCGCTTAACATATCGAAGAACgg                                            >  2:3894086/1‑93 (MQ=255)
                                              aTGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGc                                           <  2:3864890/93‑1 (MQ=255)
                                              aTGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGc                                           <  1:2174418/93‑1 (MQ=255)
                                                 tAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGa                                        >  2:2122802/1‑93 (MQ=255)
                                                 tAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGa                                        >  1:4453403/1‑93 (MQ=255)
                                                  aGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGAc                                       <  1:1539717/93‑1 (MQ=255)
                                                  aGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGAc                                       >  2:1449204/1‑93 (MQ=255)
                                                   gTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACg                                      >  2:658131/1‑93 (MQ=255)
                                                   gTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACg                                      <  1:1078242/93‑1 (MQ=255)
                                                   gTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACg                                      <  2:4120461/93‑1 (MQ=255)
                                                   gTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACg                                      >  1:2852093/1‑93 (MQ=255)
                                                     cGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGt                                    <  1:1004492/93‑1 (MQ=255)
                                                      gCCATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACggct                                   >  1:3562337/1‑91 (MQ=255)
                                                       ccATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTa                                  >  2:1107879/1‑93 (MQ=255)
                                                       ccATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTa                                  >  1:3355001/1‑93 (MQ=255)
                                                        cATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTAt                                 <  2:2047249/93‑1 (MQ=255)
                                                        cATGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTAt                                 >  1:1999255/1‑93 (MQ=255)
                                                         aTGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTAtt                                >  1:125562/1‑93 (MQ=255)
                                                          tGTAACTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTAttt                               >  2:3020090/1‑92 (MQ=255)
                                                             aaCTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTcc                            <  1:4059312/93‑1 (MQ=255)
                                                              aCTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCg                           >  2:682725/1‑93 (MQ=255)
                                                              aCTCATCGGCAAGAGGCTGTAAAACCACGTGTACGACAATCCGCGGCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTACg                           >  2:3782100/1‑93 (MQ=255)
                                                               cTCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGGTATTGTCCgg                          >  2:206325/1‑93 (MQ=255)
                                                                tCATCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGa                         >  1:1442233/1‑93 (MQ=255)
                                                                  aTCGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATa                       >  1:1380056/1‑93 (MQ=255)
                                                                    cGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGTCGGTTATTGTTCGGATAac                     <  2:4583846/93‑1 (MQ=255)
                                                                    cGGCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAac                     <  2:61105/93‑1 (MQ=255)
                                                                     ggCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAaca                    >  1:786797/1‑93 (MQ=255)
                                                                     ggCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAaca                    >  1:370984/1‑93 (MQ=255)
                                                                     gacAATAGGCTGTAAATCCACGTGTCCGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTACTGTCCGGATAaca                    <  2:4521413/91‑1 (MQ=255)
                                                                      gCAATAGGCTGTAAGTCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAacac                   <  1:1833090/93‑1 (MQ=255)
                                                                      gCAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAacac                   <  1:3901044/93‑1 (MQ=255)
                                                                       cAATAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAAcaca                  <  1:1890273/93‑1 (MQ=255)
                                                                        aaTAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAg                 <  2:1531529/93‑1 (MQ=255)
                                                                         aTAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGc                <  2:3198442/93‑1 (MQ=255)
                                                                         aTAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGc                <  2:1831537/93‑1 (MQ=255)
                                                                         aTAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGc                <  2:4265133/93‑1 (MQ=255)
                                                                          tAGGCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCg               >  1:4023630/1‑93 (MQ=255)
                                                                            ggCTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTa             <  1:3611763/93‑1 (MQ=255)
                                                                              cTGTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAAc           <  1:1629069/93‑1 (MQ=255)
                                                                                gTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACAt         <  1:2594097/93‑1 (MQ=255)
                                                                                gTAAATCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACAt         <  1:3237170/93‑1 (MQ=255)
                                                                                   aaTCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTaa      <  1:3784505/93‑1 (MQ=255)
                                                                                     tCCACGTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACAGTAAGc    >  2:1888068/1‑93 (MQ=255)
                                                                                       cacgTGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAg  >  2:2814981/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGAGTTCAGCCATTGCTGAAAATGGCGCTTAAAGCTTTTGTCATTGATGTAGTCGCCATGTAACTCATCGGCAATAGGCTGTAAATCCACATGTACGACAATCCGCGTCAGTTTACCGCTTAACATATCGAAGAACGGCTGACGGTTATTGTCCGGATAACACAGCGTAACATTAAGCAG  >  NC_000913/4043313‑4043492

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: