BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,958,877	2 bp→CG	coding (586‑587/2022 nt)	rep →	rep helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,958,877	0	G	C	100.0%	1888.3 / NA	540	A196R (GCC→CGC)	rep	rep helicase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base C (287/253); total (287/253)
*	NC_000913	3,958,878	0	C	G	100.0%	1886.4 / NA	539	A196R (GCC→CGC)	rep	rep helicase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base G (285/254); total (285/254)

Only 100 of 553 total aligned reads displayed.

TGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAAGCCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATG  >  NC_000913/3958790‑3958964

||

tGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAAt                                                                                    >  2:4726520/1‑93 (MQ=255)

  ttATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATga                                                                                  >  2:3342067/1‑93 (MQ=255)

   taggggCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaa                                                                                 <  2:1756628/90‑1 (MQ=255)

   tATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaa                                                                                 <  2:5053259/93‑1 (MQ=255)

   tATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaa                                                                                 <  2:5365266/93‑1 (MQ=255)

    agaggCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaag                                                                                <  1:3958164/90‑1 (MQ=255)

    aTGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaag                                                                                >  2:3202910/1‑93 (MQ=255)

     tGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaaga                                                                               >  1:1150759/1‑93 (MQ=255)

     tGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaaga                                                                               <  1:3228015/93‑1 (MQ=255)

     tGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaaga                                                                               <  1:6527080/93‑1 (MQ=255)

     tGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaaga                                                                               <  1:2173047/93‑1 (MQ=255)

      gggCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATgaagaa                                                                              <  1:6034535/93‑1 (MQ=255)

         cTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTc                                                                           <  2:3897812/93‑1 (MQ=255)

         cTGTATGATGCACAACTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTc                                                                           <  1:7589584/93‑1 (MQ=255)

           gTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCg                                                                         <  2:4482061/93‑1 (MQ=255)

           gTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCg                                                                         >  2:6202955/1‑93 (MQ=255)

            tATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGc                                                                        >  2:6916084/1‑93 (MQ=255)

            tATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGc                                                                        >  2:5517376/1‑93 (MQ=255)

            tATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGc                                                                        >  1:3103170/1‑93 (MQ=255)

            tATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGc                                                                        <  1:8019239/93‑1 (MQ=255)

              tgatgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCaa                                                                      <  1:4437669/93‑1 (MQ=255)

               gatgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAg                                                                     <  2:665188/93‑1 (MQ=255)

               gatgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAg                                                                     <  1:6223769/93‑1 (MQ=255)

               gatgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAg                                                                     <  2:8462770/93‑1 (MQ=255)

                 tgCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcg                                                                   <  1:1581906/93‑1 (MQ=255)

                  gCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAgcgc                                                                  >  1:2007554/1‑93 (MQ=255)

                      acCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCGGCAAGCGCCGGc                                                              >  1:5194208/1‑93 (MQ=255)

                      acCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGc                                                              <  2:2732000/93‑1 (MQ=255)

                      acCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGc                                                              >  1:5096797/1‑93 (MQ=255)

                         tGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGa                                                           <  2:1912209/93‑1 (MQ=255)

                         tGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTACCGAAGTTGCTGCTGCAACGCCATGAAGAAGTCCCCAAGCGCTGGCAGa                                                           >  2:5811375/1‑93 (MQ=255)

                          gAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCGGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGaa                                                          <  1:375490/93‑1 (MQ=255)

                            aaGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACa                                                        >  2:2886517/1‑93 (MQ=255)

                            aaGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACa                                                        >  1:7209688/1‑93 (MQ=255)

                            aaGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACa                                                        >  2:2600454/1‑93 (MQ=255)

                              gCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACaaa                                                      >  2:528832/1‑93 (MQ=255)

                                 tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAtt                                                   >  1:3886358/1‑93 (MQ=255)

                                 tGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAAtt                                                   >  1:1003616/1‑93 (MQ=255)

                                  gTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTc                                                  >  1:2436327/1‑93 (MQ=255)

                                    aaCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGc                                                <  2:3723547/93‑1 (MQ=255)

                                     aCGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCt                                               >  2:1203245/1‑93 (MQ=255)

                                      cGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTa                                              <  1:1914442/93‑1 (MQ=255)

                                       gTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTAt                                             >  2:5184443/1‑93 (MQ=255)

                                        ttCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATc                                            >  1:7059420/1‑93 (MQ=255)

                                        ttCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATc                                            <  1:8451137/93‑1 (MQ=255)

                                          ctcGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctg                                          <  2:8551862/93‑1 (MQ=255)

                                          ctcGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGGCCGCAAGCGCTGGCAGAACAAAATTCGCTATctg                                          >  2:5637915/1‑93 (MQ=255)

                                           tcGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATctgc                                         >  1:2979652/1‑93 (MQ=255)

                                              aCTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtgg                                      >  2:7084029/1‑93 (MQ=255)

                                               cTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggt                                     >  1:3497810/1‑93 (MQ=255)

                                               cTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggg                                     >  2:1611487/1‑92 (MQ=255)

                                               cTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggg                                     >  1:690329/1‑92 (MQ=255)

                                                ttCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtg                                    >  2:2126603/1‑93 (MQ=255)

                                                ttCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtg                                    <  2:4198935/93‑1 (MQ=255)

                                                 tCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCtggtgg                                   <  2:80144/93‑1 (MQ=255)

                                                  cGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGa                                  <  1:8093706/93‑1 (MQ=255)

                                                  cGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGa                                  >  1:7684621/1‑93 (MQ=255)

                                                  cGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGa                                  <  1:5728285/93‑1 (MQ=255)

                                                   gatgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGAt                                 >  1:8289110/1‑93 (MQ=255)

                                                   gatgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGAt                                 >  2:2991230/1‑93 (MQ=255)

                                                   gatgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGAt                                 <  1:4867485/93‑1 (MQ=255)

                                                    atgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATg                                <  1:2742268/93‑1 (MQ=255)

                                                    atgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATg                                >  1:735287/1‑93 (MQ=255)

                                                    atgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATg                                <  1:2375051/93‑1 (MQ=255)

                                                    atgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATg                                <  2:4574098/93‑1 (MQ=255)

                                                     tgatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGa                               >  1:5029301/1‑93 (MQ=255)

                                                      gatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAg                              >  1:2802016/1‑93 (MQ=255)

                                                      gatCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAg                              >  2:7208595/1‑93 (MQ=255)

                                                        tCTGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTa                            >  1:4515576/1‑93 (MQ=255)

                                                          tGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATc                          >  2:6734295/1‑93 (MQ=255)

                                                          tGATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATc                          >  2:8761963/1‑93 (MQ=255)

                                                           gATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCa                         <  2:6345087/93‑1 (MQ=255)

                                                           gATTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGATGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCa                         <  1:131432/93‑1 (MQ=255)

                                                            aTTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAg                        >  1:864783/1‑93 (MQ=255)

                                                            aTTTTATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAg                        <  1:8170334/93‑1 (MQ=255)

                                                             ttttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAgg                       <  1:6798499/93‑1 (MQ=255)

                                                              gttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGa                      <  1:4309790/92‑1 (MQ=255)

                                                               ttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGac                     <  1:1420729/93‑1 (MQ=255)

                                                               ttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGac                     <  2:7383163/93‑1 (MQ=255)

                                                               ttATTGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGac                     <  2:5995840/93‑1 (MQ=255)

                                                                 aTTGCCGACGTGGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacac                   <  2:3799862/93‑1 (MQ=255)

                                                                  ttGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacacc                  >  1:4727790/1‑93 (MQ=255)

                                                                  ttGCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaca                  >  2:2357031/1‑92 (MQ=255)

                                                                    gCCGACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaa                <  2:3573371/93‑1 (MQ=255)

                                                                       gACGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacac             <  1:4129171/93‑1 (MQ=255)

                                                                        aCGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCATAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacacc            >  2:830542/1‑93 (MQ=255)

                                                                        aCGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacacc            <  1:5938337/93‑1 (MQ=255)

                                                                         cgctgctgctgcAACGCAATGAAGAAGTCCGCAAGCGCAGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacacca           <  1:7881513/90‑1 (MQ=255)

                                                                         cGTTGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGacaccaacacca           <  2:1050983/93‑1 (MQ=255)

                                                                          gTTGCTGCTGCAACGCAATGAAGAAGGCCGCAAGCGCTGGCAGAACAAAATGCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccag          >  2:7840844/1‑93 (MQ=255)

                                                                           ttGCTGCTGCAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagc         >  1:5328434/1‑93 (MQ=255)

                                                                           ttGCTGCTGCAACGCAATGAAGAAGTCCGCAAGAGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGCATCAGGACACCAACAccagc         >  2:1453087/1‑93 (MQ=255)

                                                                            tgctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagcc        <  2:2245144/93‑1 (MQ=255)

                                                                            tgctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagcc        >  2:7609267/1‑93 (MQ=255)

                                                                            tgctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagcc        <  1:2309891/93‑1 (MQ=255)

                                                                            tgctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagcc        <  1:1649937/93‑1 (MQ=255)

                                                                             gctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagcca       >  1:5516047/1‑93 (MQ=255)

                                                                             gctgctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACAccagcca       <  2:5707445/93‑1 (MQ=255)

                                                                                gctgcAACGCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTa    <  2:4889294/93‑1 (MQ=255)

                                                                                  tgcAACGCAATGAAGAAGGCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATg  <  2:4375254/93‑1 (MQ=255)

||

TGTTATGGGCTGTATGATGCACACCTGAAAGCCTGTAACGTTCTCGACTTCGATGATCTGATTTTATTGCCGACGTTGCTGCTGCAAGCCAATGAAGAAGTCCGCAAGCGCTGGCAGAACAAAATTCGCTATCTGCTGGTGGATGAGTATCAGGACACCAACACCAGCCAGTATG  >  NC_000913/3958790‑3958964

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