BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	52,312	G→A	A40V (GCG→GTG)	ksgA ←	S‑adenosylmethionine‑6‑N',N'‑adenosyl (rRNA) dimethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	52,312	0	G	A	100.0%	725.5 / NA	208	A40V (GCG→GTG)	ksgA	S‑adenosylmethionine‑6‑N',N'‑adenosyl (rRNA) dimethyltransferase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (76/132); total (76/132)

Only 100 of 218 total aligned reads displayed.

TCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCT  >  NC_000913/52222‑52401

tCGCGGTGAAGTTCGATGCCCGTCAGCTCGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCAcc                                                                                         <  2:2035830/93‑1 (MQ=255)

tCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCAcc                                                                                         >  1:964924/1‑93 (MQ=255)

tCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCAcc                                                                                         <  2:2134305/93‑1 (MQ=255)

tCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATCCCGCCAGACCGGGGCCGATTTCGACCATCAcc                                                                                         <  2:5759/93‑1 (MQ=255)

 cgcgGTCAAGTTCGATGACCGTCAGCTGATCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCt                                                                                        <  1:1064611/93‑1 (MQ=255)

    ggTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGc                                                                                     <  1:309594/93‑1 (MQ=255)

    ggTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCAGCCAGACCGGGGCCGATTTCGACCATCACCTGGc                                                                                     >  1:1893087/1‑93 (MQ=255)

     gTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGcc                                                                                    >  1:981059/1‑93 (MQ=255)

     gTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCAGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTggac                                                                                    >  2:3195189/1‑91 (MQ=255)

     gTCAAGTTCGATGACCGTCAGCCGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGcc                                                                                    >  2:4347790/1‑93 (MQ=255)

        aaGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCtt                                                                                 >  2:1532000/1‑93 (MQ=255)

         aGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCttt                                                                                <  1:150211/93‑1 (MQ=255)

         aGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCttt                                                                                >  2:2285061/1‑93 (MQ=255)

         aGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCttt                                                                                <  1:3930985/93‑1 (MQ=255)

         aGTTCGATCACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCttt                                                                                <  1:4172115/93‑1 (MQ=255)

         aGTGCGATGACCGTCAGCTGGTCAAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCttt                                                                                <  2:2931628/93‑1 (MQ=255)

          gTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCtttt                                                                               <  1:1397763/93‑1 (MQ=255)

           ttCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTg                                                                              >  1:266486/1‑93 (MQ=255)

            tCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGc                                                                             <  2:3394166/93‑1 (MQ=255)

             cGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCg                                                                            <  2:3461968/93‑1 (MQ=255)

              gATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCgg                                                                           >  2:2183279/1‑93 (MQ=255)

              gATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCgg                                                                           >  1:1155028/1‑93 (MQ=255)

              gATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCgg                                                                           <  1:918142/93‑1 (MQ=255)

               aTGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCggg                                                                          <  1:2481556/93‑1 (MQ=255)

               aTGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCggg                                                                          <  1:1802194/93‑1 (MQ=255)

                ggACCGTCAGCTGGTCCAGGCGTTCGCCGACCGGTTCGGTCAATGCGGCCAAACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGt                                                                         <  2:2408713/92‑1 (MQ=255)

                 gACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGtt                                                                        <  2:2840942/93‑1 (MQ=255)

                  aCCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTa                                                                       <  1:1549461/93‑1 (MQ=255)

                   ccGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTaa                                                                      <  1:1559479/93‑1 (MQ=255)

                   ccGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTaa                                                                      <  2:1858073/93‑1 (MQ=255)

                    cGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAAt                                                                     <  2:2772659/93‑1 (MQ=255)

                    cGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAAt                                                                     <  2:1453159/93‑1 (MQ=255)

                    cGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAAt                                                                     <  2:3541792/93‑1 (MQ=255)

                    cGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAAt                                                                     <  1:242348/93‑1 (MQ=255)

                     gTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGTCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATg                                                                    <  2:3817513/93‑1 (MQ=255)

                     gTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATg                                                                    <  2:1281384/93‑1 (MQ=255)

                     gTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATg                                                                    >  1:4249637/1‑93 (MQ=255)

                     gTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATg                                                                    <  1:1282146/93‑1 (MQ=255)

                      tCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATgg                                                                   >  2:334481/1‑93 (MQ=255)

                      tCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATgg                                                                   <  1:986179/93‑1 (MQ=255)

                       cAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGc                                                                  <  1:2056461/93‑1 (MQ=255)

                        aGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGTGGGTTAATGGCa                                                                 <  1:1939821/93‑1 (MQ=255)

                        aGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCa                                                                 <  2:1716824/93‑1 (MQ=255)

                        aGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCa                                                                 <  1:3960608/93‑1 (MQ=255)

                         gCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAg                                                                >  2:405547/1‑93 (MQ=255)

                          cTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGa                                                               <  2:2996730/93‑1 (MQ=255)

                           tGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGac                                                              <  2:810298/93‑1 (MQ=255)

                              tCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGAcaca                                                           <  1:1468981/93‑1 (MQ=255)

                              tCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGAAGAcaca                                                           <  2:4454639/93‑1 (MQ=255)

                              gccAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGAcaca                                                           <  2:1248264/92‑1 (MQ=255)

                                cAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAAt                                                         >  1:4286584/1‑93 (MQ=255)

                                  gACGTTCGCCGACGGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATAc                                                       <  1:2759359/93‑1 (MQ=255)

                                  gACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATAc                                                       >  1:552104/1‑93 (MQ=255)

                                     gTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGt                                                    >  1:2846511/1‑93 (MQ=255)

                                       tCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCg                                                  >  2:735475/1‑93 (MQ=255)

                                          ccGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATc                                               >  1:4587502/1‑93 (MQ=255)

                                          ccGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATc                                               <  1:742325/93‑1 (MQ=255)

                                           cGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCa                                              <  2:1573084/93‑1 (MQ=255)

                                              ccGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGa                                           <  2:4477966/93‑1 (MQ=255)

                                               cGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGaa                                          <  2:4121977/93‑1 (MQ=255)

                                               cGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGaa                                          <  1:472180/93‑1 (MQ=255)

                                               cGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGaa                                          <  2:2084768/93‑1 (MQ=255)

                                                ggTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAAc                                         <  2:521069/93‑1 (MQ=255)

                                                 gTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACt                                        <  1:826656/93‑1 (MQ=255)

                                                  ttCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTg                                       <  2:3181204/93‑1 (MQ=255)

                                                     ggTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGCACTGATc                                    >  2:2719002/1‑93 (MQ=255)

                                                      gTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCg                                   >  1:3169924/1‑93 (MQ=255)

                                                       tCAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGt                                  >  1:2836310/1‑93 (MQ=255)

                                                        cAATGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGtt                                 >  1:3269393/1‑93 (MQ=255)

                                                         aaTGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTg                                >  1:3012767/1‑93 (MQ=255)

                                                          aTGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTga                               <  1:3156646/93‑1 (MQ=255)

                                                          aTGCCGCCAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTga                               <  1:3736861/93‑1 (MQ=255)

                                                             ccgccAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGaa                            <  2:1975079/93‑1 (MQ=255)

                                                               gccAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAg                          <  1:1672801/93‑1 (MQ=255)

                                                                ccAGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGt                         <  1:4030824/93‑1 (MQ=255)

                                                                  aGACGGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGttt                       <  1:2599952/93‑1 (MQ=255)

                                                                  aGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGttt                       <  1:2891982/93‑1 (MQ=255)

                                                                  aGACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGttt                       <  1:3976712/93‑1 (MQ=255)

                                                                   gACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGtttt                      <  1:3371846/93‑1 (MQ=255)

                                                                   gACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGtttt                      >  1:237888/1‑93 (MQ=255)

                                                                   gACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGtttt                      >  2:4400731/1‑93 (MQ=255)

                                                                   gACCGGGGCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGtttt                      <  1:2023402/93‑1 (MQ=255)

                                                                      cGGGGCCGATTTCGACAATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGcc                   >  1:3748767/1‑93 (MQ=255)

                                                                         ggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGa                <  2:1895716/93‑1 (MQ=255)

                                                                         ggCCGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGa                >  2:2581401/1‑93 (MQ=255)

                                                                            cGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAgt             >  2:3779102/1‑92 (MQ=255)

                                                                            cGATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGc             <  2:1334988/93‑1 (MQ=255)

                                                                             gATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCg            >  1:1096834/1‑93 (MQ=255)

                                                                             gATTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCg            >  1:1868269/1‑93 (MQ=255)

                                                                              aTTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGt           <  1:1018876/93‑1 (MQ=255)

                                                                              aTTTCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGt           >  2:4296066/1‑93 (MQ=255)

                                                                               tttCGACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGtt          >  2:343549/1‑93 (MQ=255)

                                                                                   gACCATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTatg      >  2:1073883/1‑91 (MQ=255)

                                                                                     ccATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACggg    >  2:372922/1‑93 (MQ=255)

                                                                                     ccATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACggg    <  1:3677133/93‑1 (MQ=255)

                                                                                     ccATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACggg    <  1:4425257/93‑1 (MQ=255)

                                                                                     ccATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACggg    >  1:3277423/1‑93 (MQ=255)

                                                                                      cATCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGc   <  1:2275440/93‑1 (MQ=255)

                                                                                       aTCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCt  >  2:4659401/1‑93 (MQ=255)

                                                                                       aTCACCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCt  >  1:2927315/1‑93 (MQ=255)

TCGCGGTCAAGTTCGATGACCGTCAGCTGGTCCAGACGTTCGCCGACCGGTTCGGTCAATGCCGCCAGACCGGGGCCGATTTCGACCATCGCCTGGCCCTTTTGCGGGTTAATGGCAGACACAATACTGTCGATCACGAACTGATCGTTGAGAAAGTTTTGCCCGAAGCGTTTACGGGCT  >  NC_000913/52222‑52401

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