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,9540AG99.6% 1461.7 / ‑3.6 460L128L (TTA→TTGaslBputative arylsulfatase regulator
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (215/243);  total (216/244)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

Only 100 of 470 total aligned reads displayed.

GAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGA  >  NC_000913/3980864‑3981042
                                                                                          |                                                                                        
gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTc                                                                                        <  1:101662/93‑1 (MQ=255)
gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTc                                                                                        >  1:6706677/1‑93 (MQ=255)
gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTc                                                                                        >  2:3595320/1‑93 (MQ=255)
gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTc                                                                                        <  2:9666336/93‑1 (MQ=255)
    aTAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTg                                                                                    <  2:5571612/93‑1 (MQ=255)
    aTAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTg                                                                                    <  2:7964010/93‑1 (MQ=255)
    aTAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTg                                                                                    <  2:7650581/93‑1 (MQ=255)
     tAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTgg                                                                                   <  2:1615226/93‑1 (MQ=255)
      aaGTAACAGCTTCCAGACTAACGGCGTGTTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGa                                                                                  <  1:2551913/93‑1 (MQ=255)
      aaGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGa                                                                                  <  2:4235750/93‑1 (MQ=255)
      aaGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGa                                                                                  <  1:8648010/93‑1 (MQ=255)
      aaGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGa                                                                                  <  1:886100/93‑1 (MQ=255)
       aGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGAt                                                                                 <  2:8302968/93‑1 (MQ=255)
       aGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGAt                                                                                 <  1:8428053/93‑1 (MQ=255)
       aGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGAt                                                                                 <  1:3277620/93‑1 (MQ=255)
        gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATg                                                                                >  1:7754981/1‑93 (MQ=255)
        gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATg                                                                                <  1:7948437/93‑1 (MQ=255)
        gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATg                                                                                >  2:6609427/1‑93 (MQ=255)
        gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATg                                                                                >  2:7061761/1‑93 (MQ=255)
         tAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATgg                                                                               <  2:7500268/93‑1 (MQ=255)
           aCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGcc                                                                             <  1:480175/93‑1 (MQ=255)
           aCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGcc                                                                             >  2:4515862/1‑93 (MQ=255)
            cAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGccc                                                                            <  1:9505455/93‑1 (MQ=255)
              gCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGc                                                                          <  1:1771795/93‑1 (MQ=255)
                 tCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTg                                                                       >  1:8567408/1‑93 (MQ=255)
                    aGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTgaga                                                                    >  1:5213838/1‑93 (MQ=255)
                     gACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGAt                                                                   <  2:3051536/93‑1 (MQ=255)
                        tAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCca                                                                >  2:1437648/1‑93 (MQ=255)
                        tAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCca                                                                >  2:8182898/1‑93 (MQ=255)
                        tAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCca                                                                <  1:3998768/93‑1 (MQ=255)
                         aaCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCcac                                                               >  2:2551067/1‑93 (MQ=255)
                          aCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCcaca                                                              >  2:9523780/1‑93 (MQ=255)
                             gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatc                                                           <  2:1952811/93‑1 (MQ=255)
                                tgctgcTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAgaatcaat                                                        >  1:2678368/1‑93 (MQ=255)
                                tgctgcTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatcaat                                                        >  2:459890/1‑93 (MQ=255)
                                tgctgcTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatcaat                                                        >  1:3729843/1‑93 (MQ=255)
                                 gctgctGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAata                                                       <  1:2959014/93‑1 (MQ=255)
                                 gctgctGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAata                                                       >  2:3355596/1‑93 (MQ=255)
                                   tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATc                                                     <  1:1905834/93‑1 (MQ=255)
                                       ggATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGt                                                 >  1:7601376/1‑93 (MQ=255)
                                        gATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTg                                                >  2:1697558/1‑93 (MQ=255)
                                         aTGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGa                                               <  1:3444041/93‑1 (MQ=255)
                                          tggcGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGAc                                              <  1:4019107/90‑1 (MQ=255)
                                           gACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACt                                             <  2:2136734/93‑1 (MQ=255)
                                              gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaaa                                          <  1:5921546/93‑1 (MQ=255)
                                                aTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATTGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAgg                                        >  2:5796726/1‑93 (MQ=255)
                                                aTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAgg                                        <  2:3810069/93‑1 (MQ=255)
                                                aTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAgg                                        <  1:6426952/93‑1 (MQ=255)
                                                 tggtgcGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGt                                       >  1:4451425/1‑93 (MQ=255)
                                                 tggtgcGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGt                                       >  1:6067831/1‑93 (MQ=255)
                                                    tgcgcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGc                                    >  1:3151126/5‑93 (MQ=255)
                                                     gcgcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCa                                   >  2:3712518/4‑93 (MQ=255)
                                                     gcgcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCa                                   <  2:5036155/90‑1 (MQ=255)
                                                        cGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAc                                <  2:7778456/93‑1 (MQ=255)
                                                         gTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAcc                               >  2:5125640/1‑93 (MQ=255)
                                                          tttCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAccc                              >  1:1056040/1‑93 (MQ=255)
                                                          tttCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAccc                              >  2:1000880/1‑93 (MQ=255)
                                                          ttcctcGCGTAGCATGATTTTCTTGTGTGTTGGTCGCTGGATGGCCCGCCGGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAccc                              <  1:1507284/90‑1 (MQ=255)
                                                           ttCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCa                             <  1:2227058/93‑1 (MQ=255)
                                                           ttCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCa                             <  1:4129421/93‑1 (MQ=255)
                                                            tctcGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCAc                            <  1:1388685/93‑1 (MQ=255)
                                                                gcgGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCAACGCAt                        <  2:5432297/93‑1 (MQ=255)
                                                                 cgGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATa                       <  2:6111728/93‑1 (MQ=255)
                                                                  ggAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATaa                      >  2:941606/1‑93 (MQ=255)
                                                                  ggAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATaa                      >  2:7415304/1‑93 (MQ=255)
                                                                  ggAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATaa                      <  1:5092681/93‑1 (MQ=255)
                                                                  ggAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATaa                      >  1:7158094/1‑93 (MQ=255)
                                                                   gAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGGTCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAg                     >  2:2374318/1‑93 (MQ=255)
                                                                   gAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAg                     >  1:6352917/1‑93 (MQ=255)
                                                                   gAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAg                     >  2:6255643/1‑93 (MQ=255)
                                                                   gAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGACCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACACACGCATAAg                     >  2:3195494/1‑93 (MQ=255)
                                                                   gAGCATCATTTTCTTGTTGGCTTGTCGCTGGATGGCTCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAg                     <  2:6925284/93‑1 (MQ=255)
                                                                    aGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGc                    <  1:4858941/93‑1 (MQ=255)
                                                                     gCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCt                   >  2:4698678/1‑93 (MQ=255)
                                                                     gCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCt                   <  1:2070860/93‑1 (MQ=255)
                                                                     gCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCt                   >  2:5922470/1‑93 (MQ=255)
                                                                       atcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTgg                 >  1:8268398/1‑93 (MQ=255)
                                                                       atcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTgg                 >  1:556998/1‑93 (MQ=255)
                                                                       atcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCCCGTGACTAAAGGTGGCACACCCACGCATAAGCTgg                 >  1:7461888/1‑93 (MQ=255)
                                                                        tcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGt                >  2:8569722/1‑93 (MQ=255)
                                                                         catTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTg               <  1:4750792/93‑1 (MQ=255)
                                                                          atTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGa              >  2:9875109/1‑93 (MQ=255)
                                                                          atTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGa              >  1:3558318/1‑93 (MQ=255)
                                                                          atTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGa              <  2:5617899/93‑1 (MQ=255)
                                                                           ttttCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGAt             <  2:4379910/93‑1 (MQ=255)
                                                                           ttttCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGAt             <  2:8722839/93‑1 (MQ=255)
                                                                              tCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCg          <  2:2941412/93‑1 (MQ=255)
                                                                              tCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCg          >  2:2737968/1‑93 (MQ=255)
                                                                               cTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGt         <  2:5191531/93‑1 (MQ=255)
                                                                               cTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGt         <  1:734154/93‑1 (MQ=255)
                                                                               cTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGt         >  1:8749915/1‑93 (MQ=255)
                                                                                  gttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGcc      >  1:6952448/1‑93 (MQ=255)
                                                                                  gttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGcc      >  2:978453/1‑93 (MQ=255)
                                                                                   ttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGccc     <  2:1216222/93‑1 (MQ=255)
                                                                                      gTTTGTCGCTGGATGGCCCGCCTGAGATCCGCAATCAATATCGCGTGACTAAAGGTGGCAGCCCCACGCATAAGCTGGTGATGCGTGCCCTGa  <  1:9048818/93‑1 (MQ=255)
                                                                                      gTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGa  <  1:9766992/93‑1 (MQ=255)
                                                                                      gTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGa  <  1:4302342/93‑1 (MQ=255)
                                                                                      gTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGa  >  2:8906626/1‑93 (MQ=255)
                                                                                      gTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGa  >  1:1242737/1‑93 (MQ=255)
                                                                                      gTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGa  >  1:9215516/1‑93 (MQ=255)
                                                                                          |                                                                                        
GAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGA  >  NC_000913/3980864‑3981042

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