BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,899,975	G→A	A452V (GCA→GTA)	bglB ←	phospho‑beta‑glucosidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,899,975	0	G	A	100.0%	953.0 / NA	258	A452V (GCA→GTA)	bglB	phospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (111/147); total (111/147)

Only 100 of 265 total aligned reads displayed.

CGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑G‑CCCCCGCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTAGCGCTTA  >  NC_000913/3899886‑3900064

   atcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATAc                                                                                      >  2:2743827/1‑93 (MQ=255)

   atcCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATAc                                                                                      >  2:4664455/1‑93 (MQ=255)

     ccAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCa                                                                                    >  1:6140671/1‑93 (MQ=255)

     ccAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCa                                                                                    <  1:5490127/93‑1 (MQ=255)

     ccAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCa                                                                                    >  1:652232/1‑93 (MQ=255)

     ccAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCa                                                                                    >  2:4282232/1‑93 (MQ=255)

       aGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATc                                                                                  <  2:2471539/93‑1 (MQ=255)

       aGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATc                                                                                  <  2:3033171/93‑1 (MQ=255)

       aGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATc                                                                                  <  2:322201/93‑1 (MQ=255)

         aCTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCg                                                                                <  1:4195586/93‑1 (MQ=255)

          cTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGa                                                                               <  1:5993986/93‑1 (MQ=255)

          cTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGa                                                                               <  1:4507710/93‑1 (MQ=255)

          cTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGa                                                                               <  1:3772702/93‑1 (MQ=255)

           tGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGaa                                                                              <  1:1099620/93‑1 (MQ=255)

             ttCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGc                                                                            <  2:4132154/93‑1 (MQ=255)

              tCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCt                                                                           <  1:1032921/93‑1 (MQ=255)

               cTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCCCCTCTACATGCCATCCGAAGCtt                                                                          <  1:3458301/93‑1 (MQ=255)

               cTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCtt                                                                          >  1:1885149/1‑93 (MQ=255)

                tGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCttt                                                                         <  2:4527366/93‑1 (MQ=255)

                 gAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCtttt                                                                        <  1:5895369/93‑1 (MQ=255)

                 gAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCtttt                                                                        <  1:5174413/93‑1 (MQ=255)

                  aaTGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTc                                                                       >  2:1491303/1‑93 (MQ=255)

                   aTGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGACGT‑C‑CCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCt                                                                      <  2:500107/93‑1 (MQ=255)

                   aTGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCt                                                                      <  1:1420408/93‑1 (MQ=255)

                   aTGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCt                                                                      <  1:490728/93‑1 (MQ=255)

                   aTGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTGAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCt                                                                      <  2:3818483/93‑1 (MQ=255)

                    tGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCtt                                                                     <  2:3934260/93‑1 (MQ=255)

                    tGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCtt                                                                     >  2:5867235/1‑93 (MQ=255)

                      cgacgaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGc                                                                   <  2:3672332/93‑1 (MQ=255)

                      cgacgaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTACTTGc                                                                   <  1:5643584/93‑1 (MQ=255)

                       gacgaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCg                                                                  <  2:1034197/93‑1 (MQ=255)

                          gaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCgtgt                                                               >  2:255303/1‑93 (MQ=255)

                          gaTAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGGCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTCCgtgt                                                               >  1:4096154/1‑93 (MQ=255)

                              aTTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCtt                                                           <  1:4748764/93‑1 (MQ=255)

                              aTTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCtt                                                           >  1:2059434/1‑93 (MQ=255)

                               ttAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTtg                                                          <  1:6202491/93‑1 (MQ=255)

                                tAAGGTTGTTTAATGGTTTTTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTtgt                                                         <  2:5602877/93‑1 (MQ=255)

                                tAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTtgt                                                         <  2:1375607/93‑1 (MQ=255)

                                  aGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGa                                                       <  2:1194330/93‑1 (MQ=255)

                                  aGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGa                                                       <  1:1166974/93‑1 (MQ=255)

                                  aGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGa                                                       <  2:1757665/93‑1 (MQ=255)

                                   ggTGCTTTAATGGTTATTTTTTTTAATGAC‑AGG‑CCCCCGCGTCCTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGGTCTTTTGAg                                                      >  2:4829121/1‑93 (MQ=255)

                                   ggTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAg                                                      <  2:2329583/93‑1 (MQ=255)

                                   agggCTTTTAAGGGTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAg                                                      <  2:1907876/90‑1 (MQ=255)

                                      gCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGct                                                   <  2:751267/93‑1 (MQ=255)

                                      gCTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGct                                                   >  1:5157260/1‑93 (MQ=255)

                                       cTTTAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCCTGGTGACCTCCTGATACCATCCGGAGCTCTTTTTGCGTGTTCTTGTggggtgt                                                  >  2:2196320/1‑90 (MQ=255)

                                        tttAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttc                                                 >  1:3294696/1‑93 (MQ=255)

                                         ttAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcc                                                >  1:4858943/1‑93 (MQ=255)

                                          tAATGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcct                                               <  2:4497235/93‑1 (MQ=255)

                                           aaTGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcctt                                              <  2:2201849/93‑1 (MQ=255)

                                           aaTGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcctt                                              <  1:2567652/93‑1 (MQ=255)

                                           aaTGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcctt                                              <  1:942263/93‑1 (MQ=255)

                                           aaTGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttcctt                                              >  1:5462941/1‑93 (MQ=255)

                                            aTGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttccttc                                             >  2:5454129/1‑93 (MQ=255)

                                            aTGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGcttccttc                                             <  1:5341690/93‑1 (MQ=255)

                                             tGGTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCg                                            >  1:5198732/1‑93 (MQ=255)

                                               gTTATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGcc                                          >  2:603343/1‑93 (MQ=255)

                                                ttATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCa                                         >  1:2230325/1‑93 (MQ=255)

                                                 tATTTTTTTTAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCAt                                        <  2:4896629/93‑1 (MQ=255)

                                                   ttttttttAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCCTTa                                      >  2:5194043/1‑93 (MQ=255)

                                                   ttttttttAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTa                                      <  1:5760034/93‑1 (MQ=255)

                                                    tttttttAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTat                                     <  1:3122026/93‑1 (MQ=255)

                                                     ttttttAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTatc                                    <  1:3938338/93‑1 (MQ=255)

                                                      tttttAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTatca                                   >  2:5357175/1‑93 (MQ=255)

                                                      tttttAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTatca                                   >  1:2539142/1‑93 (MQ=255)

                                                      tttttAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTatca                                   <  2:2001208/93‑1 (MQ=255)

                                                         ttAATGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTAtcatca                                >  2:2463133/1‑93 (MQ=255)

                                                          tAATGAC‑A‑GCCCCCCGCGGCTTGAGCACCTCTACATACCATCCGAAGGTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCCATATCAGCAc                               >  2:3526380/1‑93 (MQ=255)

                                                            aTGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGa                             >  2:2860093/1‑93 (MQ=255)

                                                            aTGAC‑A‑GCCCCCCGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGa                             <  1:2908852/93‑1 (MQ=255)

                                                                   agccccccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCcac                        >  2:1678051/4‑93 (MQ=255)

                                                                   agccccccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCcac                        <  1:4690595/90‑1 (MQ=255)

                                                                    gccccccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCcaca                       <  2:4465608/91‑1 (MQ=255)

                                                                      cccccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACata                     >  1:4867964/1‑93 (MQ=255)

                                                                       ccccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACataa                    >  1:2871498/1‑93 (MQ=255)

                                                                        cccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACataaa                   >  2:6233125/1‑93 (MQ=255)

                                                                        cccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACataaa                   <  2:1315518/93‑1 (MQ=255)

                                                                        cccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACataaa                   <  1:3317171/93‑1 (MQ=255)

                                                                         ccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACataaat                  >  2:211924/1‑93 (MQ=255)

                                                                         ccGCGTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACataaat                  >  1:5268945/1‑93 (MQ=255)

                                                                           gcgTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACataaataa                >  1:4280113/1‑93 (MQ=255)

                                                                           gcgTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACataaataa                <  1:1266899/93‑1 (MQ=255)

                                                                            cgTCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACAtaaataaa               <  1:828356/93‑1 (MQ=255)

                                                                              tCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGc             >  1:2708608/1‑93 (MQ=255)

                                                                              tCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGc             <  1:3868446/93‑1 (MQ=255)

                                                                              tCTTGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGc             <  1:4177405/93‑1 (MQ=255)

                                                                               tttGATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGcc            <  1:1453987/92‑1 (MQ=255)

                                                                                  taTCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTa         <  1:1504504/92‑1 (MQ=255)

                                                                                  gATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTa         >  1:3690512/1‑93 (MQ=255)

                                                                                  gATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTa         >  1:172314/1‑93 (MQ=255)

                                                                                  gATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTa         <  2:4844778/93‑1 (MQ=255)

                                                                                  gATCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTa         <  1:5507867/93‑1 (MQ=255)

                                                                                    tCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTAgc       >  2:5864761/1‑93 (MQ=255)

                                                                                    tCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTAgc       >  2:3000247/1‑93 (MQ=255)

                                                                                    tCACCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTAgc       >  2:265843/1‑93 (MQ=255)

                                                                                      aCCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTAgcgc     >  2:5662587/1‑93 (MQ=255)

                                                                                      aCCTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTAgcgc     >  2:1780209/1‑93 (MQ=255)

                                                                                       ccTCTACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTAGCGCt    <  2:3756091/93‑1 (MQ=255)

                                                                                         tctACATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTAGCGCTTa  <  2:702538/93‑1 (MQ=255)

CGCATCCAGACTGTTCTGAATGCGACGATAATTAAGGTGCTTTAATGGTTATTTTTTTTAATGAC‑A‑G‑CCCCCGCGTCTTGATCACCTCTGCATACCATCCGAAGCTTTTCTTGCGTGTTCTTGTGAGGCTTCCTTCGCCATTATCATCACGATCCACATAAATAAAGCCGTAGCGCTTA  >  NC_000913/3899886‑3900064

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: `‑`