BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,842,033	G→T	R67R (CGC→CGA)	hycG ←	hydrogenase activity

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,842,033	0	G	T	100.0%	912.1 / NA	243	R67R (CGC→CGA)	hycG	hydrogenase activity
Reads supporting (aligned to +/- strand): ref base G (0/0); new base T (109/134); total (109/134)

Only 100 of 247 total aligned reads displayed.

CAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAAGCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAGGC  >  NC_000913/2841943‑2842122

 agcgcagGGGATCGCATTGCACGTGTGACCGCGCCGGTAAACAGAAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGt                                                                                        <  1:571850/93‑1 (MQ=255)

  gcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCggt                                                                                       >  2:5068880/1‑91 (MQ=255)

   cgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTc                                                                                      <  1:4980136/93‑1 (MQ=255)

   cgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTc                                                                                      <  2:6331061/93‑1 (MQ=255)

    gcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCt                                                                                     >  2:4629929/1‑93 (MQ=255)

    gcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCt                                                                                     <  2:5343188/93‑1 (MQ=255)

    gcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCt                                                                                     >  1:4747976/1‑93 (MQ=255)

     cagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGtttct                                                                                    <  1:3893506/93‑4 (MQ=255)

     cagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTg                                                                                    <  2:6061758/93‑1 (MQ=255)

      agGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGc                                                                                   <  1:5724734/93‑1 (MQ=255)

      agGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGc                                                                                   >  2:806080/1‑93 (MQ=255)

       ggggATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCa                                                                                  <  2:4240298/93‑1 (MQ=255)

       ggggATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCa                                                                                  >  1:6212497/1‑93 (MQ=255)

        gggATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCAt                                                                                 <  2:1716685/93‑1 (MQ=255)

         ggATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATc                                                                                >  2:1920084/1‑93 (MQ=255)

               cATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAg                                                                          >  2:1586650/1‑93 (MQ=255)

                 ttGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcg                                                                        <  2:2751707/93‑1 (MQ=255)

                 ttGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcg                                                                        <  1:4496217/93‑1 (MQ=255)

                  tGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcgg                                                                       <  2:3164466/93‑1 (MQ=255)

                   gCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcggc                                                                      <  2:2716346/93‑1 (MQ=255)

                   gCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcggc                                                                      >  1:5840723/1‑93 (MQ=255)

                       gggTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGaaa                                                                  <  1:1854144/93‑1 (MQ=255)

                        ggTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAg                                                                 <  1:5592624/93‑1 (MQ=255)

                        ggTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAg                                                                 <  1:6048642/93‑1 (MQ=255)

                         gTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGc                                                                >  2:5264322/1‑93 (MQ=255)

                          tGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCg                                                               <  2:3270216/93‑1 (MQ=255)

                          tGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCg                                                               <  1:3683317/93‑1 (MQ=255)

                           gACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGt                                                              <  2:4600484/93‑1 (MQ=255)

                           gACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGt                                                              >  1:3379639/1‑93 (MQ=255)

                            aCCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTg                                                             >  1:1005291/1‑93 (MQ=255)

                             ccGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTgg                                                            >  1:4352257/1‑93 (MQ=255)

                              cgcgGCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGc                                                           <  2:2446372/93‑1 (MQ=255)

                              cgcgCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGc                                                           <  2:2352630/93‑1 (MQ=255)

                                cgcCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGa                                                         <  2:4140898/93‑1 (MQ=255)

                                  ccGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGaaa                                                       >  1:3822392/1‑93 (MQ=255)

                                   cGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGaaaa                                                      <  2:4699132/93‑1 (MQ=255)

                                   cGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGaaaa                                                      <  2:2292625/93‑1 (MQ=255)

                                   cGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGaaaa                                                      <  2:867317/93‑1 (MQ=255)

                                    ggTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAat                                                     <  2:3651790/93‑1 (MQ=255)

                                     gTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatt                                                    <  1:2446709/93‑1 (MQ=255)

                                      tAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAattt                                                   >  1:6197282/1‑93 (MQ=255)

                                       aaaCAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttc                                                  >  2:5207635/1‑93 (MQ=255)

                                       aaaCAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttc                                                  >  1:1435176/1‑93 (MQ=255)

                                        aaCAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttct                                                 >  1:2630953/1‑92 (MQ=255)

                                        aaCAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcg                                                 >  1:1286262/1‑93 (MQ=255)

                                         aCAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcga                                                <  2:6022100/93‑1 (MQ=255)

                                          cAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgat                                               >  1:2873197/1‑93 (MQ=255)

                                           aGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgatt                                              <  1:5151876/93‑1 (MQ=255)

                                           aGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgatt                                              >  1:1241392/1‑93 (MQ=255)

                                             tAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgatttc                                            <  1:4066603/93‑1 (MQ=255)

                                             tAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgatttc                                            >  1:3664314/1‑93 (MQ=255)

                                              aaaaTATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgatttcg                                           >  1:2474582/1‑93 (MQ=255)

                                               aaaTATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgatttcga                                          >  1:2210098/1‑92 (MQ=255)

                                                aaTATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCa                                         >  1:1374664/1‑93 (MQ=255)

                                                aaTATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCa                                         >  2:1845084/1‑93 (MQ=255)

                                                 atatCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCaa                                        >  2:2241532/1‑93 (MQ=255)

                                                  tatCCGCGTGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAAc                                       <  2:4951921/93‑1 (MQ=255)

                                                  tatCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAAc                                       <  2:6228839/93‑1 (MQ=255)

                                                  tatCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAAc                                       <  1:5421404/93‑1 (MQ=255)

                                                    tCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCg                                     <  2:1479534/93‑1 (MQ=255)

                                                     ccGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGt                                    >  2:3737903/1‑93 (MQ=255)

                                                     ccGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGt                                    <  2:6713137/93‑1 (MQ=255)

                                                      cGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGtt                                   <  2:2108390/93‑1 (MQ=255)

                                                      cGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGtt                                   <  1:6172459/93‑1 (MQ=255)

                                                       gCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTg                                  >  1:876138/1‑93 (MQ=255)

                                                        cATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGc                                 <  1:6329242/93‑1 (MQ=255)

                                                         aTGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCa                                >  1:1704297/1‑93 (MQ=255)

                                                          tGACGCGGGGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAg                               >  2:1701036/1‑93 (MQ=255)

                                                           gACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgc                              <  2:5449080/93‑1 (MQ=255)

                                                            aCGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgcc                             <  1:762575/93‑1 (MQ=255)

                                                                ggTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGccgccg                         >  1:1627440/1‑93 (MQ=255)

                                                                ggTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGccgccg                         >  2:4045711/1‑93 (MQ=255)

                                                                  tGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCa                       <  2:6510238/93‑1 (MQ=255)

                                                                   gAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAg                      <  1:500607/93‑1 (MQ=255)

                                                                    aaGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAgc                     >  1:1739295/1‑92 (MQ=255)

                                                                    aaGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGt                     >  2:5831100/1‑93 (MQ=255)

                                                                     aGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTc                    <  2:5674518/93‑1 (MQ=255)

                                                                      ggAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTcc                   <  2:4875750/93‑1 (MQ=255)

                                                                      ggAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTcc                   >  1:103895/1‑93 (MQ=255)

                                                                      ggAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTcc                   >  1:6651440/1‑93 (MQ=255)

                                                                      ggAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTcc                   >  1:2509660/1‑93 (MQ=255)

                                                                       gAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCa                  <  2:3421427/93‑1 (MQ=255)

                                                                       gAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCa                  <  2:1508332/93‑1 (MQ=255)

                                                                        aaCGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTcaaa                 >  2:4420542/1‑90 (MQ=255)

                                                                        aaCGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCAc                 <  1:3075657/93‑1 (MQ=255)

                                                                          cGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCAcgc               <  2:6651266/93‑1 (MQ=255)

                                                                           gACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCAcgcg              <  2:3333012/93‑1 (MQ=255)

                                                                            aCTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCgg             >  1:5092819/1‑93 (MQ=255)

                                                                            aCTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCCTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAACCCACCCgg             >  1:1942649/1‑93 (MQ=38)

                                                                               ttAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTaa          >  1:429334/1‑93 (MQ=255)

                                                                               ttAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCAACGCGGTaa          >  2:1020370/1‑93 (MQ=255)

                                                                                tAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTaaa         >  1:2082623/1‑93 (MQ=255)

                                                                                 aaTGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAAc        <  2:550607/93‑1 (MQ=255)

                                                                                  aTGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACa       >  2:5699836/1‑93 (MQ=255)

                                                                                  aTGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACa       <  1:2422613/93‑1 (MQ=255)

                                                                                  aTGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACa       <  1:2686526/93‑1 (MQ=255)

                                                                                     ccGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAg    <  2:1246975/93‑1 (MQ=255)

                                                                                     ccGAATCGTTCTGCAACAAGCAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAg    <  2:3451396/93‑1 (MQ=255)

                                                                                      cGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAgg   <  1:4213783/93‑1 (MQ=255)

                                                                                      cGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAgg   <  1:4395826/93‑1 (MQ=255)

CAGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAAGCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAGGC  >  NC_000913/2841943‑2842122

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