BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,039,960	Δ1 bp	intergenic (+42/‑36)	yihD → / → yihE	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,039,959	0	A	.	100.0%	1555.5 / NA	334	intergenic (+41/‑37)	yihD/yihE	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base . (162/172); total (162/172)

Only 100 of 354 total aligned reads displayed.

AGAAGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGG  >  NC_000913/4039867‑4040047

   agaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGAGGCGCCACCAAAATCATCGTGAAATGATatcc                                                                                       >  2:1245397/1‑90 (MQ=255)

    gaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatcct                                                                                      >  1:5453302/1‑89 (MQ=255)

    gaCTTCAAAACGGCGCTGTTACGAGCTCGCGGCGTAATTAAAGAGCAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatcct                                                                                      >  1:5374225/1‑89 (MQ=255)

      cTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTc                                                                                   <  1:1611888/93‑1 (MQ=255)

      cTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTc                                                                                   <  2:4626773/93‑1 (MQ=255)

      cTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTc                                                                                   <  2:5192942/93‑1 (MQ=255)

       ttCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCg                                                                                  >  1:401501/1‑93 (MQ=255)

        tCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGt                                                                                 <  2:4882189/93‑1 (MQ=255)

         cAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTc                                                                                <  1:4405949/93‑1 (MQ=255)

         cAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTc                                                                                >  2:6702720/1‑93 (MQ=255)

         cAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTc                                                                                <  2:2359343/93‑1 (MQ=255)

          aaaaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCa                                                                               >  2:5367842/1‑93 (MQ=255)

          aaaaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCa                                                                               >  1:3515656/1‑93 (MQ=255)

          aaaaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCa                                                                               >  1:6076068/1‑93 (MQ=255)

          aaaaCGGCGATGTTACGCGCTCGCGGCGTACTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCa                                                                               <  2:4352704/93‑1 (MQ=255)

           aaaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCAt                                                                              >  2:276450/1‑93 (MQ=255)

           aaaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCAt                                                                              <  2:1828276/93‑1 (MQ=255)

            aaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCAtt                                                                             >  1:5277662/1‑93 (MQ=255)

            aaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCAtt                                                                             >  1:1007156/1‑93 (MQ=255)

             aCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTc                                                                            >  1:25771/1‑93 (MQ=255)

             aCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTc                                                                            >  1:5694374/1‑93 (MQ=255)

                gcgcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTa                                                                         >  1:3154116/1‑93 (MQ=255)

                gcgcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTa                                                                         <  1:1753585/93‑1 (MQ=255)

                 ggcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCCCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTaa                                                                        <  2:4699336/92‑1 (MQ=255)

                 cgcTGTTGCGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTaa                                                                        <  2:6348021/93‑1 (MQ=255)

                 cgcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTaa                                                                        <  2:603153/93‑1 (MQ=255)

                    tGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGt                                                                     <  2:5718500/93‑1 (MQ=255)

                    tGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGt                                                                     >  1:1196741/1‑93 (MQ=255)

                      ttACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGttt                                                                   <  2:1401747/93‑1 (MQ=255)

                      ttACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGttt                                                                   <  2:4186241/93‑1 (MQ=255)

                       tACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGtttt                                                                  >  2:2541654/1‑93 (MQ=255)

                        aCGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTc                                                                 <  2:18997/93‑1 (MQ=255)

                        aCGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTc                                                                 <  2:5023504/93‑1 (MQ=255)

                         cgcgCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTcc                                                                <  1:2821525/93‑1 (MQ=255)

                          gcgcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCg                                                               <  2:3601752/93‑1 (MQ=255)

                            gcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGga                                                             >  2:4111666/1‑93 (MQ=255)

                             cccgcgGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGgat                                                            <  2:3453869/91‑1 (MQ=255)

                             cTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGgat                                                            >  2:2199249/1‑93 (MQ=255)

                             cTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGgat                                                            <  2:527392/93‑1 (MQ=255)

                             cTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGgat                                                            <  2:365950/93‑1 (MQ=255)

                              tCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGgatg                                                           >  2:5462976/1‑93 (MQ=255)

                              tCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGgatg                                                           >  2:2054278/1‑93 (MQ=255)

                                gcggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGgatgat                                                         <  1:3613324/93‑1 (MQ=255)

                                  ggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATgg                                                       <  1:4748108/93‑1 (MQ=255)

                                  ggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATgg                                                       >  1:1375163/1‑93 (MQ=255)

                                  ggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATgg                                                       >  2:4805349/1‑93 (MQ=255)

                                    cgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGa                                                     <  1:4823414/93‑1 (MQ=255)

                                     gTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGAt                                                    <  1:3343042/93‑1 (MQ=255)

                                       aaTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGa                                                  >  1:4983265/1‑93 (MQ=255)

                                       aaTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGa                                                  >  1:3771819/1‑93 (MQ=255)

                                        aTTAAAGAGTAAATGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGaa                                                 <  1:5149464/93‑1 (MQ=255)

                                        aTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGaa                                                 <  2:6214644/93‑1 (MQ=255)

                                        aTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGaa                                                 >  2:5734611/1‑93 (MQ=255)

                                        aTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGaa                                                 >  1:5097464/1‑93 (MQ=255)

                                         ttAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAAt                                                <  2:3914600/93‑1 (MQ=255)

                                             agagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACa                                            <  1:184896/93‑1 (MQ=255)

                                             agagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACa                                            >  1:3029415/1‑93 (MQ=255)

                                             agagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACa                                            >  1:4157703/1‑93 (MQ=255)

                                              gagTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAg                                           <  2:2492038/93‑1 (MQ=255)

                                                 tAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcgc                                        >  1:2177613/1‑93 (MQ=255)

                                                  aaaaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCt                                       <  2:5945081/93‑1 (MQ=255)

                                                   aaaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCtt                                      <  1:6689891/93‑1 (MQ=255)

                                                    aaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCttt                                     >  1:2747298/1‑93 (MQ=255)

                                                    aaGCTGGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCttt                                     <  2:1008584/93‑1 (MQ=255)

                                                         tGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACtt                                <  1:1249945/93‑1 (MQ=255)

                                                         tGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACtt                                >  1:4906742/1‑93 (MQ=255)

                                                         tGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACtt                                <  2:3016538/93‑1 (MQ=255)

                                                         tGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACtt                                <  1:198097/93‑1 (MQ=255)

                                                         tGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACtt                                >  1:423120/1‑93 (MQ=255)

                                                           tAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTc                              <  2:6096299/93‑1 (MQ=255)

                                                            aaGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTcc                             >  2:770714/1‑93 (MQ=255)

                                                              gcggcgCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAg                           >  2:6386591/1‑93 (MQ=255)

                                                               cggcgCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGa                          <  2:312063/93‑1 (MQ=255)

                                                               cggcgCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGa                          <  1:5891076/93‑1 (MQ=255)

                                                                 gcgcCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGaca                        >  1:3734398/1‑93 (MQ=255)

                                                                 gcgcCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGaca                        <  2:3330858/93‑1 (MQ=255)

                                                                    ccaccaAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTa                     >  2:1486534/1‑93 (MQ=255)

                                                                    ccaccaAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTa                     >  1:6411639/1‑93 (MQ=255)

                                                                     caccaAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTac                    >  2:3967773/1‑93 (MQ=255)

                                                                     caccaAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTac                    <  2:5568673/93‑1 (MQ=255)

                                                                      accaAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTaca                   >  1:2433685/1‑93 (MQ=255)

                                                                      accaAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTaca                   >  2:664216/1‑93 (MQ=255)

                                                                      accaAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTaca                   >  1:248195/1‑93 (MQ=255)

                                                                       ccaAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTacac                  >  1:731827/1‑93 (MQ=255)

                                                                         aaaaTCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAccc                >  2:2825991/1‑93 (MQ=255)

                                                                          aaaTCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCg               >  1:2267113/1‑93 (MQ=255)

                                                                           aaTCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCgg              >  1:4519367/1‑93 (MQ=255)

                                                                            atcatcGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGa             >  2:3725885/1‑93 (MQ=255)

                                                                             tcatcGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGAt            >  1:963414/1‑93 (MQ=255)

                                                                              catcGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATa           >  2:1665274/1‑93 (MQ=255)

                                                                               atcGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATAc          >  1:2132916/1‑93 (MQ=255)

                                                                               atcGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATAc          >  1:1072658/1‑93 (MQ=255)

                                                                                tcGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATAcc         >  2:1057458/1‑93 (MQ=255)

                                                                                 cGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATACCAGCGC‑TTTACTTTCCAGACACTACACCCGGATACcat       <  1:5618182/93‑1 (MQ=255)

                                                                                   tGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatc      >  1:2018199/1‑93 (MQ=255)

                                                                                   tGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatc      <  1:6515831/93‑1 (MQ=255)

                                                                                    gAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatca     >  2:3824024/1‑93 (MQ=255)

                                                                                      aaTGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATg   >  2:1763433/1‑93 (MQ=255)

                                                                                      aaTGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATg   >  2:2643057/1‑93 (MQ=255)

                                                                                      aaTGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATg   <  2:4257696/93‑1 (MQ=255)

AGAAGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGG  >  NC_000913/4039867‑4040047

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