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%	954.4 / NA	296	L128L (TTA→TTG)	aslB	putative arylsulfatase regulator
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (129/167); total (129/167)

Only 100 of 301 total aligned reads displayed.

GGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGAC  >  NC_000913/3980863‑3981043

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

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

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

    gATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCt                                                                                      <  2:4631429/93‑1 (MQ=255)

      tAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTgg                                                                                    <  2:1246069/93‑1 (MQ=255)

      tAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTgg                                                                                    <  1:2581001/93‑1 (MQ=255)

       aaGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGa                                                                                   <  1:1948583/93‑1 (MQ=255)

        aGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGAt                                                                                  <  2:3459669/93‑1 (MQ=255)

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

         gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATg                                                                                 >  1:3902873/1‑93 (MQ=255)

         gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATg                                                                                 >  1:60389/1‑93 (MQ=255)

          tAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATgg                                                                                >  2:6629144/1‑93 (MQ=255)

          tAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATgg                                                                                <  2:6157062/93‑1 (MQ=255)

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

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

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

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

            aCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGcc                                                                              <  2:976700/93‑1 (MQ=255)

            aCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGcc                                                                              >  2:2037091/1‑93 (MQ=255)

             cAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGccc                                                                             >  2:6202874/1‑93 (MQ=255)

              aGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCg                                                                            <  1:1694883/93‑1 (MQ=255)

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

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

                cTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGcc                                                                          >  2:6546901/1‑93 (MQ=255)

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

                 ttCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCt                                                                         <  1:2282648/93‑1 (MQ=255)

                            cGGCGTGCTGCTGGATGACGAATGGTGCGCATTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaa                                                              <  2:4063727/93‑1 (MQ=255)

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

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

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

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

                                gTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatcaa                                                          <  1:655310/93‑1 (MQ=255)

                                gTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatcaa                                                          <  1:5479231/93‑1 (MQ=255)

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

                                   ctgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAatat                                                       <  1:2949381/93‑1 (MQ=255)

                                   ctgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAatat                                                       <  1:6444460/93‑1 (MQ=255)

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

                                    tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATc                                                      >  2:1755546/1‑93 (MQ=255)

                                    tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATc                                                      >  1:1229380/1‑93 (MQ=255)

                                     gctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcg                                                     <  1:453009/93‑1 (MQ=255)

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

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

                                        ggATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGt                                                  <  2:223854/93‑1 (MQ=255)

                                        ggATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGt                                                  <  1:736650/93‑1 (MQ=255)

                                         gATGACGAATGGTGCGCGTTTCTCGCTGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTg                                                 <  1:5614649/93‑1 (MQ=255)

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

                                           tGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGAc                                               <  2:6076184/93‑1 (MQ=255)

                                           tGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGAc                                               <  1:5079353/93‑1 (MQ=255)

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

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

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

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

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

                                                  tggtgcGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGt                                        <  1:903303/93‑1 (MQ=255)

                                                  tggtgcGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGt                                        <  2:2556801/93‑1 (MQ=255)

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

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

                                                          gTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAcc                                <  2:2508589/93‑1 (MQ=255)

                                                          gTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAcc                                >  1:4878408/1‑93 (MQ=255)

                                                           tttCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAccc                               <  2:5297494/93‑1 (MQ=255)

                                                           tttCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAccc                               <  2:2332651/93‑1 (MQ=255)

                                                            ttCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCa                              <  1:338058/93‑1 (MQ=255)

                                                             tctcGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCAc                             >  1:1433128/1‑93 (MQ=255)

                                                             tctcGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCACGTGACTAAAGGTGGCAGACCCAc                             >  2:4880086/1‑93 (MQ=255)

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

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

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

                                                                cgcgGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCa                          <  1:4066061/93‑1 (MQ=255)

                                                                cgcgGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGGAGCAGACCCACGCa                          >  1:6265852/1‑93 (MQ=255)

                                                                 gcgGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCAt                         <  1:4934197/93‑1 (MQ=255)

                                                                   ggAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATaa                       <  2:1704302/93‑1 (MQ=255)

                                                                   ggAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATaa                       >  1:941612/1‑93 (MQ=255)

                                                                   ggAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATaa                       >  2:3641597/1‑93 (MQ=255)

                                                                    gggCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGT‑GCAGACCCACGCATAAGc                     <  1:2489941/91‑1 (MQ=255)

                                                                    gAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAg                      >  1:4117117/1‑93 (MQ=255)

                                                                    gAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGACTGAGGTCCACAATCAATATCGCGTGACTAAAGGTGGCAGACTTCCGCATacg                      >  2:4754505/1‑91 (MQ=255)

                                                                     aGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATCAGc                     >  1:2382110/1‑93 (MQ=255)

                                                                      gCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCt                    <  2:4411049/93‑1 (MQ=255)

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

                                                                        atcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTgg                  >  2:2125040/1‑93 (MQ=255)

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

                                                                           atTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGa               >  1:1329556/1‑93 (MQ=255)

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

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

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

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

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

                                                                               tCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCg           >  1:1016987/1‑93 (MQ=255)

                                                                                cTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGt          <  2:1321256/93‑1 (MQ=255)

                                                                                cTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGt          >  2:5990994/1‑93 (MQ=255)

                                                                                 ttgttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTg         >  2:574451/1‑93 (MQ=255)

                                                                                 ttgttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTg         <  2:3928232/93‑1 (MQ=255)

                                                                                    ttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGccc      <  1:1539326/93‑1 (MQ=255)

                                                                                    ttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGccc      <  1:3434546/93‑1 (MQ=255)

                                                                                    ttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGccc      <  2:6348765/93‑1 (MQ=255)

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

                                                                                     tgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCt     >  1:3538664/1‑93 (MQ=255)

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

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

                                                                                        tttGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGAc  >  2:4581782/1‑93 (MQ=255)

GGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGAC  >  NC_000913/3980863‑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: `‑`