BRESEQ :: Evidence

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_000913	3,980,954	0	A	G	100.0%	769.7 / NA	244	L128L (TTA→TTG)	aslB	putative arylsulfatase regulator
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (111/133); total (111/133)

Only 100 of 251 total aligned reads displayed.

GAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGAC  >  NC_000913/3980864‑3981043

gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTc                                                                                         <  2:160719/93‑1 (MQ=255)

gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTc                                                                                         <  1:1444298/93‑1 (MQ=255)

 aaGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCg                                                                                        <  1:2229688/93‑1 (MQ=255)

  aGATAAGTAACAGCTTCCAGACTAGCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGc                                                                                       <  2:3894882/93‑1 (MQ=255)

  aGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGc                                                                                       <  1:3603993/93‑1 (MQ=255)

  aGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGc                                                                                       <  1:3831734/93‑1 (MQ=255)

  aGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGc                                                                                       >  2:2198436/1‑93 (MQ=255)

     tAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTgg                                                                                    >  1:1603603/1‑93 (MQ=255)

      aaGTAACATCATCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGa                                                                                   <  1:350646/93‑1 (MQ=255)

       aGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGAt                                                                                  >  2:4424794/1‑93 (MQ=255)

        gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGATGGTTg                                                                                 >  2:1654734/1‑93 (MQ=255)

        gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATg                                                                                 <  2:177751/93‑1 (MQ=255)

        gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATg                                                                                 <  2:3025579/93‑1 (MQ=255)

          aaCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCCTTTTCTTGTTGGTTTGTCGCTGGATGGc                                                                               >  1:2420051/1‑93 (MQ=255)

          aaCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGc                                                                               >  2:1443281/1‑93 (MQ=255)

          aaCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGc                                                                               <  1:58081/93‑1 (MQ=255)

              gCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGc                                                                           >  1:4494287/1‑93 (MQ=255)

               cTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGcc                                                                          <  1:1168843/93‑1 (MQ=255)

               cTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGcc                                                                          <  1:4404233/93‑1 (MQ=255)

                 tCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTg                                                                        >  2:75948/1‑93 (MQ=255)

                 tCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTg                                                                        >  1:647367/1‑93 (MQ=255)

                  ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTga                                                                       <  2:3438178/93‑1 (MQ=255)

                  ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTga                                                                       <  1:842828/93‑1 (MQ=255)

                  ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTga                                                                       <  2:2929705/93‑1 (MQ=255)

                   cAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTgag                                                                      <  1:4088841/93‑1 (MQ=255)

                    aGACTAACGGCGTGCTGCTGGATGACGAATGGGGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTgaga                                                                     <  2:4361957/93‑1 (MQ=255)

                     gACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGAt                                                                    >  2:2365922/1‑93 (MQ=255)

                     gACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGAt                                                                    <  1:1262509/93‑1 (MQ=255)

                      aCTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATc                                                                   <  2:3528959/93‑1 (MQ=255)

                        tAACGGCGTTCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCca                                                                 <  1:2998680/93‑1 (MQ=255)

                        tAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCca                                                                 <  1:4002882/93‑1 (MQ=255)

                         aaCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCcac                                                                >  1:2924988/1‑93 (MQ=255)

                          aCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCcaca                                                               >  2:3829742/1‑93 (MQ=255)

                          aCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCcaca                                                               >  1:1275379/1‑93 (MQ=255)

                            ggCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaat                                                             <  1:1657099/93‑1 (MQ=255)

                             gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGCGATCCAcaatc                                                            >  1:1721260/1‑93 (MQ=255)

                             gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatc                                                            <  2:4193096/93‑1 (MQ=255)

                               gTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatcaa                                                          >  1:2076672/1‑93 (MQ=255)

                               gTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatcaa                                                          >  2:2863149/1‑93 (MQ=255)

                                tgctgcTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatcaat                                                         >  2:4540507/1‑93 (MQ=255)

                                 gctgctGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAata                                                        <  2:2307614/93‑1 (MQ=255)

                                 gctgctGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAata                                                        <  2:1933205/93‑1 (MQ=255)

                                  ctgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAatat                                                       >  1:3439912/1‑93 (MQ=255)

                                   tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATc                                                      <  1:4052224/93‑1 (MQ=255)

                                   tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATc                                                      <  1:2839598/93‑1 (MQ=255)

                                   tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATc                                                      <  1:3206699/93‑1 (MQ=255)

                                   tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATc                                                      <  2:606257/93‑1 (MQ=255)

                                      tgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcgcg                                                   >  2:2147581/1‑93 (MQ=255)

                                       ggATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGt                                                  >  1:250392/1‑93 (MQ=255)

                                        gATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTg                                                 >  1:740633/1‑93 (MQ=255)

                                         aTGACGAATGGTGCGCGTTTCTCGCGGGGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGa                                                >  1:1008353/1‑93 (MQ=255)

                                         aTGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGa                                                <  1:1439579/93‑1 (MQ=255)

                                         aTGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGa                                                >  2:1755867/1‑93 (MQ=255)

                                          tGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGAc                                               >  2:4323158/1‑93 (MQ=255)

                                          tGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGAc                                               >  2:2824482/1‑93 (MQ=255)

                                           gACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACt                                              >  2:3497442/1‑93 (MQ=255)

                                           gACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACt                                              <  1:1938961/93‑1 (MQ=255)

                                             cGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaa                                            >  2:1579587/1‑93 (MQ=255)

                                             cGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaa                                            <  2:482342/93‑1 (MQ=255)

                                             cGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaa                                            >  2:1837180/1‑93 (MQ=255)

                                              gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaaa                                           <  2:163693/93‑1 (MQ=255)

                                              gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaaa                                           >  1:11309/1‑93 (MQ=255)

                                               aaTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAg                                          <  2:1213148/93‑1 (MQ=255)

                                                aTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAgg                                         <  1:3563481/93‑1 (MQ=255)

                                                aTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAgg                                         >  1:1642896/1‑93 (MQ=255)

                                                     gcgcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCa                                    >  2:2726117/4‑93 (MQ=255)

                                                     gcgcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCa                                    <  1:4598855/90‑1 (MQ=255)

                                                      ggcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAg                                   <  1:253017/91‑1 (MQ=255)

                                                       gcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGa                                  <  1:933681/92‑1 (MQ=255)

                                                       gcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGa                                  >  2:3945388/2‑93 (MQ=255)

                                                        cGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAc                                 <  1:4585069/93‑1 (MQ=255)

                                                         gTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGGGGCAGaac                                >  2:2532436/1‑91 (MQ=255)

                                                         gTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGGGGCAGAcc                                >  1:4171308/1‑93 (MQ=255)

                                                              tcGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGc                           >  1:3864162/1‑93 (MQ=255)

                                                              tcGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGc                           <  1:2027241/93‑1 (MQ=255)

                                                                ccgGAGCATCATCGTCTTGTTCGTTTGTCGCTGGATGGTCCGACTGAGAGCCGCAATCAATAGCGCGTGTCTAAAGGTGGCAGACCCACGCag                         <  2:538001/92‑2 (MQ=255)

                                                                   gAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAg                      <  1:3237471/93‑1 (MQ=255)

                                                                    aGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGc                     >  1:3071955/1‑93 (MQ=255)

                                                                      catcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGTTCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTg                   >  1:3297585/1‑93 (MQ=255)

                                                                      catcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTg                   <  2:862447/93‑1 (MQ=255)

                                                                        tcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGt                 >  1:3049746/1‑93 (MQ=255)

                                                                        tcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGt                 >  1:2591043/1‑93 (MQ=255)

                                                                          atTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAGCCACGCATAAGCTGGTGa               >  2:2007606/1‑93 (MQ=255)

                                                                          atTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGa               <  1:178058/93‑1 (MQ=255)

                                                                          atTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGa               <  2:3162161/93‑1 (MQ=255)

                                                                          atTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCCGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGa               <  1:1498795/93‑1 (MQ=255)

                                                                           ttttCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGCCCCACGCATAAGCTGGTGAt              <  1:174371/93‑1 (MQ=255)

                                                                           ttttCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGAt              >  2:3617708/1‑93 (MQ=255)

                                                                           ttttCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGAt              <  2:1537291/93‑1 (MQ=255)

                                                                           ttttCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGAt              <  1:3945963/93‑1 (MQ=255)

                                                                            tttCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATg             >  2:2326325/1‑93 (MQ=255)

                                                                             ttCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGc            <  2:1333213/93‑1 (MQ=255)

                                                                              tCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCg           <  2:4625758/93‑1 (MQ=255)

                                                                                  gttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGcc       <  1:551875/93‑1 (MQ=255)

                                                                                  gttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGcc       >  2:57399/1‑93 (MQ=255)

                                                                                   ttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGccc      >  2:1382621/1‑93 (MQ=255)

                                                                                    tgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCt     <  1:1091382/93‑1 (MQ=255)

                                                                                      gTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGa   <  2:1486205/93‑1 (MQ=255)

                                                                                       tttGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGAc  <  1:895637/93‑1 (MQ=255)

                                                                                       tttGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGAc  <  2:17165/93‑1 (MQ=255)

GAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGAC  >  NC_000913/3980864‑3981043

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