BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,847,139	A→G	G42G (GGT→GGC)	hycC ←	membrane‑spanning protein of hydrogenase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,847,139	0	A	G	100.0%	892.3 / NA	273	G42G (GGT→GGC)	hycC	membrane‑spanning protein of hydrogenase 3
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (123/150); total (123/150)

Only 100 of 285 total aligned reads displayed.

CACATCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTACCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACAAACCA  >  NC_000913/2847052‑2847228

cacaTCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAAc                                                                                      <  1:2238331/93‑1 (MQ=255)

 ccaTCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAGCCGTGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAAcc                                                                                     <  1:2442164/92‑1 (MQ=255)

 acaTCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAAcc                                                                                     <  1:3514436/93‑1 (MQ=255)

  caTCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACcg                                                                                    <  1:4617583/93‑1 (MQ=255)

   aTCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACcgc                                                                                   >  2:4520555/1‑93 (MQ=255)

   aTCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACcgc                                                                                   >  1:2377414/1‑93 (MQ=255)

    tCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACcgcg                                                                                  >  1:896497/1‑93 (MQ=255)

    tCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACcgcg                                                                                  >  2:3302045/1‑93 (MQ=255)

      gTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACAGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGcgcc                                                                                >  1:4184813/1‑93 (MQ=255)

           ttcccAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCaa                                                                           <  2:4286740/90‑1 (MQ=255)

           ttACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCaa                                                                           >  2:5302230/1‑93 (MQ=255)

           ttACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCaa                                                                           >  1:5247757/1‑93 (MQ=255)

             aCCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAAtt                                                                         <  1:2011493/93‑1 (MQ=255)

              ccAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTc                                                                        >  2:5420801/1‑93 (MQ=255)

                aGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACCGCCAACCGCGCCGCCAATTccc                                                                      >  1:4695458/1‑92 (MQ=255)

                 gCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAg                                                                     <  1:2840682/93‑1 (MQ=255)

                      aGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATc                                                                >  2:3432482/1‑93 (MQ=255)

                      aGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATc                                                                <  2:728618/93‑1 (MQ=255)

                      aGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATc                                                                >  2:5909931/1‑93 (MQ=255)

                         gcACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAg                                                             >  2:5786479/1‑93 (MQ=255)

                          cACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGc                                                            <  2:5640626/93‑1 (MQ=255)

                             cGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCAc                                                         >  2:4156341/1‑93 (MQ=255)

                             cGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCAc                                                         >  1:4431071/1‑93 (MQ=255)

                             cGCTCACGCCAACCGCGACAGTCAGTACAGGGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCAc                                                         <  1:3166973/93‑1 (MQ=255)

                               cTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACtt                                                       >  1:1020245/1‑93 (MQ=255)

                                tCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTa                                                      >  1:2337902/1‑93 (MQ=255)

                                tCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTa                                                      <  1:1625419/93‑1 (MQ=255)

                                tCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTa                                                      >  1:3815534/1‑93 (MQ=255)

                                ccACGCAAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAACCACTTa                                                      <  2:6197679/92‑1 (MQ=255)

                                    gCCACCCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGc                                                  <  2:933643/93‑1 (MQ=255)

                                    gCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTCTCCAGCCACTTAACGc                                                  >  1:1118520/1‑93 (MQ=255)

                                    gCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGc                                                  <  2:4077033/93‑1 (MQ=255)

                                     ccAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCt                                                 <  1:1889101/93‑1 (MQ=255)

                                      cAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCtt                                                <  1:4113011/93‑1 (MQ=255)

                                        aCCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCtttt                                              >  1:2733725/1‑93 (MQ=255)

                                        aCCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCtttt                                              <  1:3622836/93‑1 (MQ=255)

                                         ccGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCttttt                                             >  2:2494029/1‑93 (MQ=255)

                                          cgcgCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTg                                            <  1:4048625/93‑1 (MQ=255)

                                             gcCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaa                                         <  1:5630078/93‑1 (MQ=255)

                                              ccAGTCAGTACAGTGAAGCCCGCGTCGGCCGTATACAGACGGCCAACCGCGCGGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaa                                        <  2:158895/93‑1 (MQ=255)

                                              ccAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaa                                        >  1:3234598/1‑93 (MQ=255)

                                              ccAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaa                                        <  1:4394482/93‑1 (MQ=255)

                                               cagtcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaag                                       >  2:3507381/1‑93 (MQ=255)

                                                 gtcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaagaa                                     >  2:334597/1‑93 (MQ=255)

                                                 gtcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaagaa                                     <  2:5279653/93‑1 (MQ=255)

                                                  tcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaagaaa                                    <  1:5396767/93‑1 (MQ=255)

                                                  tcagtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTgaaaagaaa                                    <  1:5538014/93‑1 (MQ=255)

                                                    agtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGaaaagaaaag                                  <  2:3331808/93‑1 (MQ=255)

                                                     gtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAaaagaaaaga                                 <  1:1595242/93‑1 (MQ=255)

                                                     gtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAaaagaaaaga                                 <  2:4784991/93‑1 (MQ=255)

                                                     gtACAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAaaagaaaaga                                 <  2:1188654/93‑1 (MQ=255)

                                                       aCAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAgaga                               <  1:1398610/93‑1 (MQ=255)

                                                       aCAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAgaga                               <  1:3977908/93‑1 (MQ=255)

                                                        cAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGaa                              >  2:4505369/1‑93 (MQ=255)

                                                        cAGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGaa                              <  1:3103724/93‑1 (MQ=255)

                                                        cAGTGAAGCCCGCGCCTGCCGTATACAGCCTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGaa                              <  2:2237041/93‑1 (MQ=255)

                                                         aGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCTAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGaaa                             <  1:3372483/93‑1 (MQ=255)

                                                         aGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGaaa                             <  1:614498/93‑1 (MQ=255)

                                                         aGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGaaa                             <  2:4350953/93‑1 (MQ=255)

                                                         aGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGaaa                             >  2:2565118/1‑93 (MQ=255)

                                                         aGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGaaa                             <  2:5225359/93‑1 (MQ=255)

                                                         aGTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGaaa                             <  1:6255986/93‑1 (MQ=255)

                                                          gTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAAt                            >  2:3177031/1‑93 (MQ=255)

                                                          gTGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAAt                            <  2:5985284/93‑1 (MQ=255)

                                                           tGAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATg                           <  2:6254821/93‑1 (MQ=255)

                                                           ggAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATg                           <  2:2158796/92‑1 (MQ=255)

                                                            gAAGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGc                          <  2:1962315/93‑1 (MQ=255)

                                                              aGCCCGCGGCTGCCGTATACAGACTGCCACCCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCa                        <  1:2672425/93‑1 (MQ=255)

                                                              aGCCCGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCa                        <  2:4667234/93‑1 (MQ=255)

                                                                cccGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGa                      <  1:5190586/93‑1 (MQ=255)

                                                                 ccGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGaa                     <  2:4887098/93‑1 (MQ=255)

                                                                 ccGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGaa                     <  1:6040446/93‑1 (MQ=255)

                                                                 ccGCGGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGaa                     <  1:2764352/93‑1 (MQ=255)

                                                                  cgcgGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAAc                    >  2:6070559/1‑93 (MQ=255)

                                                                   gcgGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACa                   >  2:2689398/1‑93 (MQ=255)

                                                                   gcgGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACa                   <  1:858185/93‑1 (MQ=255)

                                                                   gcgGCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACa                   <  1:4976983/93‑1 (MQ=255)

                                                                     ggCTGCCGTATACAGACTGCCAACCGGGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgc                 <  2:4117832/93‑1 (MQ=255)

                                                                     ggCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgc                 <  1:12092/93‑1 (MQ=255)

                                                                     ggCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgc                 <  2:1841117/93‑1 (MQ=255)

                                                                     ggCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgc                 <  2:1253903/93‑1 (MQ=255)

                                                                      gCTGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcg                <  2:457437/93‑1 (MQ=255)

                                                                        tGCCGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcggc              <  1:257681/93‑1 (MQ=255)

                                                                          ccGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcggcgg            <  2:815/93‑1 (MQ=255)

                                                                          ccGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcggcgg            >  2:365984/1‑93 (MQ=255)

                                                                          ccGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcggccg            <  2:3974254/93‑3 (MQ=255)

                                                                           cGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcggcggc           <  1:4351519/93‑1 (MQ=255)

                                                                           cGTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAgcggcggc           <  2:2445104/93‑1 (MQ=255)

                                                                            gTATACAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGACAGAACAgcggcggcg          >  1:5004303/1‑93 (MQ=255)

                                                                             tataCAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGa         >  1:4245722/1‑93 (MQ=255)

                                                                             tataCAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGa         >  1:4268154/1‑93 (MQ=255)

                                                                             tataCAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGa         >  1:1938006/1‑93 (MQ=255)

                                                                             tataCAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGa         >  2:2857286/1‑93 (MQ=255)

                                                                              ataCAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGAc        >  1:5751545/1‑93 (MQ=255)

                                                                               taCAGACTGCAAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACa       >  2:3856005/1‑93 (MQ=255)

                                                                                 cAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACaaa     >  1:5794026/1‑93 (MQ=255)

                                                                                 cAGACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACaaa     >  2:3755710/1‑93 (MQ=255)

                                                                                  aGACTGCCACCCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACAAAc    <  1:1267423/93‑1 (MQ=255)

                                                                                   gACTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACAAAcc   >  1:3770300/1‑93 (MQ=255)

                                                                                    aCTGCCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACAAAcca  >  1:3285208/1‑93 (MQ=255)

CACATCGTAGCTTACCAGCGACAGCGCACCGCTCACGCCAACCGCGCCAGTCAGTACAGTGAAGCCCGCGGCTGCCGTATACAGACTACCAACCGCGCCGCCAATTCCAGCTATCCAGCCACTTAACGCTTTTTGAAAAGAAAAGAGAAATGCCAGAACAGCGGCGGCGACAAACCA  >  NC_000913/2847052‑2847228

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