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%	602.1 / NA	165	R67R (CGC→CGA)	hycG	hydrogenase activity
Reads supporting (aligned to +/- strand): ref base G (0/0); new base T (85/80); total (85/80)

Only 100 of 170 total aligned reads displayed.

AGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAAGCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAGGC  >  NC_000913/2841944‑2842122

agcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGt                                                                                        <  1:3210340/93‑1 (MQ=255)

 gcgcagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGtt                                                                                       >  1:2044704/1‑93 (MQ=255)

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

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

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

    cagGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTg                                                                                    >  2:974581/1‑93 (MQ=255)

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

        ggATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATc                                                                                <  2:508506/93‑1 (MQ=255)

           tCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCaaa                                                                             <  1:495087/93‑1 (MQ=255)

            cGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAAc                                                                            >  2:1075534/1‑93 (MQ=255)

             gCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATCACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACa                                                                           <  2:2323667/93‑1 (MQ=255)

              cATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAg                                                                          <  2:3296520/93‑1 (MQ=255)

               aTTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgc                                                                         >  2:2818358/1‑93 (MQ=255)

               aTTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgc                                                                         <  1:1147812/93‑1 (MQ=255)

               aTTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgc                                                                         >  1:2965385/1‑93 (MQ=255)

                 tGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAgcgg                                                                       >  2:1837423/1‑93 (MQ=255)

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

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

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

                       ggTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAg                                                                 >  1:1061577/1‑93 (MQ=255)

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

                        gTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGc                                                                <  1:1732275/93‑1 (MQ=255)

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

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

                           aCCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTg                                                             <  1:3158104/93‑1 (MQ=255)

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

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

                              gcgcCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCg                                                          >  1:3559131/1‑93 (MQ=255)

                                gcCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAGAGCGTGGCGaa                                                        >  2:2254353/1‑93 (MQ=255)

                                gcCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCCAGAGCCAGGCGaa                                                        >  2:2254352/1‑93 (MQ=255)

                                  cGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGaaaa                                                      >  1:3069383/1‑93 (MQ=255)

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

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

                                   ggTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAat                                                     >  2:4149669/1‑93 (MQ=255)

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

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

                                      aaaCAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttc                                                  <  2:1300372/93‑1 (MQ=255)

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

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

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

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

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

                                           gTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgattt                                             >  1:4563993/1‑93 (MQ=255)

                                           gTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAatttcgattt                                             >  1:3496774/1‑93 (MQ=255)

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

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

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

                                              aaaTATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGc                                          >  2:3505964/1‑93 (MQ=255)

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

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

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

                                                 tatCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAAc                                       >  2:1175513/1‑93 (MQ=255)

                                                   tCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCg                                     >  1:329356/1‑93 (MQ=255)

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

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

                                                     cGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGtt                                   >  2:731766/1‑93 (MQ=255)

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

                                                     cGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGtt                                   >  2:2946317/1‑93 (MQ=255)

                                                     cGCATGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGtt                                   >  2:1125577/1‑93 (MQ=255)

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

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

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

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

                                                        aTGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCa                                <  2:4378982/93‑1 (MQ=255)

                                                         tGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAg                               >  2:598427/1‑93 (MQ=255)

                                                         tGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAg                               <  1:2450938/93‑1 (MQ=255)

                                                         tGACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAg                               >  2:2633949/1‑93 (MQ=255)

                                                          gACGCGGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgc                              >  1:4021119/1‑93 (MQ=255)

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

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

                                                             gcgGTGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAgccgc                           >  1:3859896/1‑93 (MQ=255)

                                                                 tGAAGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCa                       >  2:2500436/1‑93 (MQ=255)

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

                                                                    aGGAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTc                    >  1:3052941/1‑93 (MQ=255)

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

                                                                      gaccGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCa                  <  2:211768/90‑1 (MQ=255)

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

                                                                      gAACGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCa                  >  2:3505895/1‑93 (MQ=255)

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

                                                                        aCGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCAcg                <  1:1005782/93‑1 (MQ=255)

                                                                        aCGACTTTAATGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCAcg                <  1:1582784/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                 aTGCCGAATCGTTCTACATCAAACGGCGGCGAAAGCATGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACa       <  1:2209690/93‑1 (MQ=38)

                                                                                  tGCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACAt      <  2:3824338/93‑1 (MQ=255)

                                                                                   gCCGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATa     <  1:1366490/93‑1 (MQ=255)

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

                                                                                    ccGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAg    >  2:238186/1‑93 (MQ=255)

                                                                                    ccGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAg    >  1:139533/1‑93 (MQ=255)

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

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

                                                                                     cGAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAgg   >  1:1958789/1‑93 (MQ=255)

                                                                                      gAATCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAGGc  >  1:2713561/1‑93 (MQ=255)

AGCGCAGGGGATCGCATTGCACGGGTGACCGCGCCGGTAAACAGTAAAATATCCGCATGACGCGGTGAAGGAACGACTTTAATGCCGAAGCGTTCTGCATCAAACAGCGGCGAAAGCGTGGCGAAAATTTCGATTTCGCAACCGTTGCAGCCGCCGCAGTCCACGCGGTAAACATAGGC  >  NC_000913/2841944‑2842122

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