BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,901,321	T→G	A3A (GCA→GCC)	bglB ←	phospho‑beta‑glucosidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,901,321	0	T	G	100.0%	959.6 / NA	290	A3A (GCA→GCC)	bglB	phospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand): ref base T (0/0); new base G (155/135); total (155/135)

Only 100 of 298 total aligned reads displayed.

TGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTAATAC  >  NC_000913/3901231‑3901410

tGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGc                                                                                         >  1:3852472/1‑93 (MQ=255)

 gAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCt                                                                                        <  1:1423861/93‑1 (MQ=255)

 gAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCt                                                                                        >  1:5221298/1‑93 (MQ=255)

 gAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCt                                                                                        <  2:5985411/93‑1 (MQ=255)

  aGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGGCACCCCAAAGAAATGTTTCTGGAAAGGCtt                                                                                       >  1:2515290/1‑93 (MQ=255)

  aGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCtt                                                                                       <  1:829339/93‑1 (MQ=255)

    gTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTc                                                                                     <  1:3768569/93‑1 (MQ=255)

     tCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCa                                                                                    <  1:3983780/93‑1 (MQ=255)

     tCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCa                                                                                    <  1:5284819/93‑1 (MQ=255)

       gagaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATa                                                                                  <  1:6432074/93‑1 (MQ=255)

        aggtCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATaa                                                                                 <  2:2599497/90‑1 (MQ=255)

        agaTGCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATaa                                                                                 <  1:438453/93‑1 (MQ=255)

        agaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATaa                                                                                 <  2:653305/93‑1 (MQ=255)

        agaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATaa                                                                                 <  1:3526564/93‑1 (MQ=255)

        agaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATaa                                                                                 <  2:4938440/93‑1 (MQ=255)

         gaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAAt                                                                                <  2:2339057/93‑1 (MQ=255)

         gaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAAt                                                                                >  1:171197/1‑93 (MQ=255)

         gaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAAt                                                                                >  2:1025182/1‑93 (MQ=255)

         gaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAAt                                                                                >  1:1364354/1‑93 (MQ=255)

          aTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCCAAGAAATGTTTCTGGGAAGGGTTTCATAAtt                                                                               >  2:5543128/1‑93 (MQ=255)

          aTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAAtt                                                                               >  1:5745128/1‑93 (MQ=255)

          aTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAAtt                                                                               <  2:2127178/93‑1 (MQ=255)

           tCCCTTTGCCATCTTCCTGCCAGGCACCTTTCACCTGATTGGCAGCTGTTGCGGCACCCCAAAGAAAAGTTTCTGGAAAGGCTTTCATAAATa                                                                              >  1:3572215/1‑93 (MQ=255)

           tCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTa                                                                              <  2:5556065/93‑1 (MQ=255)

           tCCCTTT‑CCATCTTCCTGCCAGGCGCCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTaa                                                                             <  1:159578/93‑1 (MQ=255)

           acccTTTTCCATCTTCCTGCCAGGCACCTTCCCCCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTa                                                                              <  2:5250783/92‑1 (MQ=255)

            cccTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTaa                                                                             <  1:6329453/93‑1 (MQ=255)

            cccTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTaa                                                                             <  2:6134646/93‑1 (MQ=255)

            cccTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTaa                                                                             <  2:822188/93‑1 (MQ=255)

             ccTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAAc                                                                            <  1:1895001/93‑1 (MQ=255)

             ccTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAAc                                                                            >  1:2213526/1‑93 (MQ=255)

             ccTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAAc                                                                            >  1:2741076/1‑93 (MQ=255)

              cTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACt                                                                           >  2:4356660/1‑93 (MQ=255)

              cTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACt                                                                           <  2:5860474/93‑1 (MQ=255)

              cTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACt                                                                           <  2:1250165/93‑1 (MQ=255)

                ttGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTcc                                                                         <  2:5156220/93‑1 (MQ=255)

                  gCCATCTCCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCCCCCCAAAGAAATGTTTCTGCAAAGGCTTTCATAATTAACTCCtt                                                                       <  1:2314129/93‑1 (MQ=255)

                    cATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCtttt                                                                     >  1:2676657/1‑93 (MQ=255)

                     aTCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTa                                                                    >  1:1975889/1‑93 (MQ=255)

                     aTCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTa                                                                    >  2:4745403/1‑93 (MQ=255)

                        ttCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCg                                                                 >  2:3931983/1‑93 (MQ=255)

                          ccTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGtt                                                               <  1:5412261/93‑1 (MQ=255)

                                aGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGaa                                                         <  1:2927505/93‑1 (MQ=255)

                                 ggCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTTCTTTTATCGTTAGCGAat                                                        >  2:3459257/1‑93 (MQ=255)

                                 ggCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAat                                                        >  2:1635625/1‑93 (MQ=255)

                                 ggCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAat                                                        >  2:1822703/1‑93 (MQ=255)

                                 ggCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAat                                                        >  2:5037123/1‑93 (MQ=255)

                                  gCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAatg                                                       >  1:2786731/1‑93 (MQ=255)

                                  gCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAatg                                                       >  2:292150/1‑93 (MQ=255)

                                    aCCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTGCTTTTATCGTTAGCGAatgat                                                     >  1:5468874/1‑93 (MQ=255)

                                      cTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATgg                                                   >  1:4295502/1‑93 (MQ=255)

                                       ttCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGa                                                  >  1:4897014/1‑93 (MQ=255)

                                          cACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATaa                                               <  2:897244/93‑1 (MQ=255)

                                          cACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATaa                                               >  1:1513256/1‑93 (MQ=255)

                                           aCCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAAc                                              >  2:673187/1‑93 (MQ=255)

                                           aCCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAAc                                              <  1:1228607/93‑1 (MQ=255)

                                           aCCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAAc                                              >  1:5874217/1‑93 (MQ=255)

                                            ccTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACa                                             >  1:784219/1‑93 (MQ=255)

                                            ccTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACa                                             <  2:4922693/93‑1 (MQ=255)

                                            ccTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACa                                             >  1:4194367/1‑93 (MQ=255)

                                              tGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGc                                           <  1:3104780/93‑1 (MQ=255)

                                                aTTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCgg                                         >  1:3681699/1‑93 (MQ=255)

                                                 ttGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGt                                        <  2:6447541/93‑1 (MQ=255)

                                                  tGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGtt                                       <  1:21443/93‑1 (MQ=255)

                                                   ggCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTc                                      >  2:4543339/1‑93 (MQ=255)

                                                     cAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCAc                                    >  1:4760524/1‑93 (MQ=255)

                                                     cAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCAc                                    <  2:241777/93‑1 (MQ=255)

                                                      aGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCAcc                                   <  1:5104042/93‑1 (MQ=255)

                                                      aGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCAcc                                   >  2:2833478/1‑93 (MQ=255)

                                                       gCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCt                                  <  2:2003977/93‑1 (MQ=255)

                                                             gcgcCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCtt                            <  2:1465835/93‑1 (MQ=255)

                                                             gcgcCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCtt                            <  1:6120488/93‑1 (MQ=255)

                                                             gcgcCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCtt                            >  1:4622436/1‑93 (MQ=255)

                                                             gcgcCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCtt                            <  1:5645551/93‑1 (MQ=255)

                                                              cgcCACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTa                           <  2:4194888/93‑1 (MQ=255)

                                                                ccACCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATc                         >  2:5340998/1‑93 (MQ=255)

                                                                 cACCCCAAGGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCt                        >  1:6569943/1‑93 (MQ=255)

                                                                  aCCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTg                       >  1:118306/1‑93 (MQ=255)

                                                                  aCCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTg                       >  2:3689494/1‑93 (MQ=255)

                                                                  aCCCCAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTg                       >  1:5056107/1‑93 (MQ=255)

                                                                   ccccAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTgc                      >  2:3029283/1‑93 (MQ=255)

                                                                      cAAAGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgcc                   <  2:3832251/93‑1 (MQ=255)

                                                                       aaaGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccg                  <  1:5995961/93‑1 (MQ=255)

                                                                       aaaGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccg                  >  2:1916588/1‑93 (MQ=255)

                                                                        aaGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccgt                 <  2:6600952/93‑1 (MQ=255)

                                                                         aGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccgtg                >  1:1156540/1‑93 (MQ=255)

                                                                         aGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccgtg                >  2:3119980/1‑93 (MQ=255)

                                                                         aGAAATGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccgtg                >  2:3388963/1‑93 (MQ=255)

                                                                             aTGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccgtgccgt            >  2:2662828/1‑93 (MQ=255)

                                                                             aTGTTTCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCgccgtgccgt            <  1:4103232/93‑1 (MQ=255)

                                                                                tttCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTggg         >  1:6710458/1‑93 (MQ=255)

                                                                                 ttCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTgggg        >  1:4462626/1‑93 (MQ=255)

                                                                                  tCTGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGt       <  2:5452073/93‑1 (MQ=255)

                                                                                    tGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTaa     >  2:1757430/1‑93 (MQ=255)

                                                                                    tGGAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTaa     >  1:6504350/1‑93 (MQ=255)

                                                                                     ggAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTAAt    <  2:2234791/93‑1 (MQ=255)

                                                                                      gAAAGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTAATa   >  1:5528602/1‑93 (MQ=255)

                                                                                       aaaGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTAATAc  <  1:6193486/93‑1 (MQ=255)

                                                                                       aaaGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTAATAc  <  1:5775063/93‑1 (MQ=255)

                                                                                       aaaGGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTAATAc  >  1:5141034/1‑93 (MQ=255)

TGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCGCTTATCTGCGCCGTGCCGTGGGGTAATAC  >  NC_000913/3901231‑3901410

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