Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,980,954 A→G L128L (TTA→TTG aslB → putative arylsulfatase regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,980,9540AG100.0% 808.4 / NA 254L128L (TTA→TTGaslBputative arylsulfatase regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (116/138);  total (116/138)

Only 100 of 258 total aligned reads displayed.

GAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGAC  >  NC_000913/3980864‑3981043
                                                                                          |                                                                                         
gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTc                                                                                         <  1:1431938/93‑1 (MQ=255)
gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTc                                                                                         >  2:4577285/1‑93 (MQ=255)
gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTc                                                                                         <  1:3901875/93‑1 (MQ=255)
  aggtAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGc                                                                                       >  2:5488419/4‑93 (MQ=255)
   gATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCt                                                                                      <  2:2839570/93‑1 (MQ=255)
   gATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCt                                                                                      <  1:869082/93‑1 (MQ=255)
    aTAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTg                                                                                     <  1:4053758/93‑1 (MQ=255)
     tAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTgg                                                                                    <  2:3089222/93‑1 (MQ=255)
      aaGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGa                                                                                   <  1:105215/93‑1 (MQ=255)
       aGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCACTTTCTTGTTGGTTTGTCGCTGGAt                                                                                  <  2:6072648/93‑1 (MQ=255)
        gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATg                                                                                 <  1:3871594/93‑1 (MQ=255)
         tAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATgg                                                                                >  2:879515/1‑93 (MQ=255)
         tAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATgg                                                                                <  1:4532434/93‑1 (MQ=255)
          aaCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGc                                                                               <  2:5872203/93‑1 (MQ=255)
           aCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGcc                                                                              <  1:2455446/93‑1 (MQ=255)
             aGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCg                                                                            >  1:4168724/1‑93 (MQ=255)
               cTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGcc                                                                          <  1:2198079/93‑1 (MQ=255)
               cTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGcc                                                                          <  2:2257777/93‑1 (MQ=255)
                 tCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTg                                                                        >  2:2966712/1‑93 (MQ=255)
                   cAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGGCCGCCTgag                                                                      >  2:205244/1‑93 (MQ=255)
                    aGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTgaga                                                                     <  1:5636396/93‑1 (MQ=255)
                    aGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTgaga                                                                     >  2:3819656/1‑93 (MQ=255)
                      aCTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATc                                                                   <  1:4494343/93‑1 (MQ=255)
                        tAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCca                                                                 >  2:5848721/1‑93 (MQ=255)
                        tAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCAATTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCca                                                                 >  2:1741516/1‑93 (MQ=255)
                         aaCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCcac                                                                <  1:2233858/93‑1 (MQ=255)
                         aaCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCcac                                                                >  1:4938907/1‑93 (MQ=255)
                          aCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCcaca                                                               >  2:357308/1‑93 (MQ=255)
                           cGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaa                                                              <  2:1703295/93‑1 (MQ=255)
                           cGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaa                                                              <  2:4721964/93‑1 (MQ=255)
                            ggCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaat                                                             <  1:333765/93‑1 (MQ=255)
                             gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatc                                                            >  1:4448541/1‑93 (MQ=255)
                              cGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatca                                                           <  2:3633362/93‑1 (MQ=255)
                              cGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatca                                                           >  1:5648017/1‑93 (MQ=255)
                               gTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatcaa                                                          <  1:5271682/93‑1 (MQ=255)
                                tgctgcTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatcaat                                                         <  1:3503757/93‑1 (MQ=255)
                                 gctgctGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAata                                                        <  2:754148/93‑1 (MQ=255)
                                 gctgctGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAata                                                        <  1:1203009/93‑1 (MQ=255)
                                  ctgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAatat                                                       <  2:2552387/93‑1 (MQ=255)
                                   tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATc                                                      >  1:3121991/1‑93 (MQ=255)
                                    gctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcg                                                     >  1:3445267/1‑93 (MQ=255)
                                       ggATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcgcgg                                                  >  2:5171352/1‑92 (MQ=255)
                                          tGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGAc                                               <  2:1186386/93‑1 (MQ=255)
                                           gACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACt                                              >  2:4958099/1‑93 (MQ=255)
                                              gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaaa                                           >  2:3448598/1‑93 (MQ=255)
                                              gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaaa                                           <  1:1987966/93‑1 (MQ=255)
                                               aaTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAg                                          >  1:2568811/1‑93 (MQ=255)
                                               aaTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAg                                          <  2:3251815/93‑1 (MQ=255)
                                               aaTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAg                                          >  2:633863/1‑93 (MQ=255)
                                                aTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAgg                                         <  2:3888142/93‑1 (MQ=255)
                                                 tggtgcGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGt                                        <  2:2655705/93‑1 (MQ=255)
                                                 tggtgcGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGt                                        <  1:4271880/93‑1 (MQ=255)
                                                    tgcgcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGc                                     >  1:3286067/5‑93 (MQ=255)
                                                     gcgcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCa                                    <  1:2720030/90‑1 (MQ=255)
                                                     gcgcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCa                                    <  1:519879/90‑1 (MQ=255)
                                                         gTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGGCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAcc                                >  2:5691130/1‑93 (MQ=255)
                                                          tttCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAccc                               >  1:440899/1‑93 (MQ=255)
                                                          tttCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAccc                               >  2:426555/1‑93 (MQ=255)
                                                           ttCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCa                              <  2:2172564/93‑1 (MQ=255)
                                                           tcctcGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCa                              <  1:4357119/91‑1 (MQ=255)
                                                            tctcGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCAc                             >  1:2549817/1‑93 (MQ=255)
                                                            tctcGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCAc                             <  1:5910278/93‑1 (MQ=255)
                                                             ctcGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACg                            >  1:377217/1‑93 (MQ=255)
                                                              tcGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGc                           >  1:1594440/1‑93 (MQ=255)
                                                              tcGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGc                           <  1:545005/93‑1 (MQ=255)
                                                                gcgGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCAt                         >  2:5061513/1‑93 (MQ=255)
                                                                gcgGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCAt                         >  1:181236/1‑93 (MQ=255)
                                                                  ggAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATaa                       >  2:573132/1‑93 (MQ=255)
                                                                  ggAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATaa                       >  1:3515187/1‑93 (MQ=255)
                                                                   gAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAg                      >  2:403436/1‑93 (MQ=255)
                                                                   gAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAg                      >  1:1011399/1‑93 (MQ=255)
                                                                   gAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAg                      >  1:6034246/1‑93 (MQ=255)
                                                                    aGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGc                     <  1:2748736/93‑1 (MQ=255)
                                                                     gCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCt                    >  1:6128581/1‑93 (MQ=255)
                                                                      catcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTg                   <  2:2630441/93‑1 (MQ=255)
                                                                      catcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTg                   <  1:2933533/93‑1 (MQ=255)
                                                                      catcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTg                   <  1:1482046/93‑1 (MQ=255)
                                                                      catcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGGGGCAGACCCACGCATAAGCTg                   >  2:1493218/1‑93 (MQ=255)
                                                                       atcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTgg                  <  1:2103082/93‑1 (MQ=255)
                                                                       atcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTgg                  <  2:3713156/93‑1 (MQ=255)
                                                                        tcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGt                 <  2:5615234/93‑1 (MQ=255)
                                                                         catTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATTCACAATCAATATCGCGTGACTAAAGGATGCAGGCAGACGCAAAAGCTtgtg                >  2:4406905/1‑93 (MQ=255)
                                                                         catTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTg                >  1:5952070/1‑93 (MQ=255)
                                                                          atTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGa               >  2:1250686/1‑93 (MQ=255)
                                                                          atTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGa               >  2:4395545/1‑93 (MQ=255)
                                                                          atTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGa               <  2:6221156/93‑1 (MQ=255)
                                                                          atTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGa               >  2:1333978/1‑93 (MQ=255)
                                                                               cTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGt          <  2:4168724/93‑1 (MQ=255)
                                                                               cTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGt          >  2:1234713/1‑93 (MQ=255)
                                                                                ttgttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTg         <  2:812205/93‑1 (MQ=255)
                                                                                  gttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGcc       >  1:6271068/1‑93 (MQ=255)
                                                                                   ttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGccc      <  1:1277173/93‑1 (MQ=255)
                                                                                   ttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGccc      <  2:2371341/93‑1 (MQ=255)
                                                                                   ttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGccc      >  1:4787003/1‑93 (MQ=255)
                                                                                   ttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGccc      <  2:51320/93‑1 (MQ=255)
                                                                                    tgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCt     <  2:3469170/93‑1 (MQ=255)
                                                                                     ggTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTg    >  1:4041139/1‑93 (MQ=255)
                                                                                      gTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGa   <  2:144305/93‑1 (MQ=255)
                                                                                       tttGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGAc  >  1:1861627/1‑93 (MQ=255)
                                                                                       tttGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGAc  <  2:1016763/93‑1 (MQ=255)
                                                                                          |                                                                                         
GAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGAC  >  NC_000913/3980864‑3981043

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