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%	1508.3 / NA	325	intergenic (+41/‑37)	yihD/yihE	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base . (150/175); total (150/175)

Only 100 of 347 total aligned reads displayed.

AGAAGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGG  >  NC_000913/4039867‑4040047

 gaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATat                                                                                         <  2:6860612/93‑2 (MQ=255)

 gaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATat                                                                                         >  2:2105499/1‑92 (MQ=255)

  aagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatc                                                                                        <  2:1299972/93‑3 (MQ=255)

  aagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatc                                                                                        <  2:2903884/93‑3 (MQ=255)

   agaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatcc                                                                                       <  1:4687199/93‑4 (MQ=255)

   agaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatcc                                                                                       >  1:6911225/1‑90 (MQ=255)

     aCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCtt                                                                                    <  1:1599537/93‑1 (MQ=255)

     aCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCtt                                                                                    <  2:3850028/93‑1 (MQ=255)

     aCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCtt                                                                                    >  2:3083193/1‑93 (MQ=255)

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

          aaaaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCa                                                                               <  1:1758387/93‑1 (MQ=255)

          aaaaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCa                                                                               <  2:6386241/93‑1 (MQ=255)

          aaaaCGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCa                                                                               <  1:3393707/93‑1 (MQ=255)

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

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

              cGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCg                                                                           >  1:2695446/1‑93 (MQ=255)

                  gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAAt                                                                       >  1:1124121/1‑93 (MQ=255)

                  gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAAt                                                                       <  2:4127790/93‑1 (MQ=255)

                  gcTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAAt                                                                       >  2:2641544/1‑93 (MQ=255)

                   cTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTTTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATg                                                                      >  1:7215629/1‑93 (MQ=255)

                   cTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATg                                                                      >  2:194761/1‑93 (MQ=255)

                   cTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATg                                                                      >  1:4058785/1‑93 (MQ=255)

                   cTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATg                                                                      >  2:3774792/1‑93 (MQ=255)

                       tACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGtttt                                                                  <  2:2334705/93‑1 (MQ=255)

                        nngcgcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTc                                                                 <  1:3492867/91‑1 (MQ=255)

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

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

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

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

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

                            gcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGga                                                             <  2:5018080/93‑1 (MQ=255)

                            gcTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGga                                                             <  2:4710968/93‑1 (MQ=255)

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

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

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

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

                                 cggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGatgatg                                                        <  1:5990240/93‑1 (MQ=255)

                                 cggcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGatgatg                                                        <  1:7069471/93‑1 (MQ=255)

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

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

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

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

                                   gcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATggg                                                      <  1:6707531/93‑1 (MQ=255)

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

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

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

                                      tAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATg                                                   <  1:5298305/93‑1 (MQ=255)

                                      tAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATg                                                   >  1:1208511/1‑93 (MQ=255)

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

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

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

                                            aaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAAc                                             <  1:2246407/93‑1 (MQ=255)

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

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

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

                                               agTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgc                                          <  1:4258905/93‑1 (MQ=255)

                                                gTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcg                                         >  2:1823224/1‑93 (MQ=255)

                                                gTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAgcg                                         >  1:2409214/1‑93 (MQ=255)

                                                  aaaaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCt                                       >  1:650437/1‑93 (MQ=255)

                                                   aaaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCtt                                      >  2:5248287/1‑93 (MQ=255)

                                                    aaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGNttt                                     >  1:6064825/1‑93 (MQ=255)

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

                                                    aaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCttt                                     <  2:4916108/93‑1 (MQ=255)

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

                                                    aaGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCttt                                     <  1:7485107/93‑1 (MQ=255)

                                                      gCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTa                                   >  1:3166429/1‑93 (MQ=255)

                                                       cTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTAc                                  <  2:5225907/93‑1 (MQ=255)

                                                        ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACt                                 <  2:743058/93‑1 (MQ=255)

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

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

                                                           tAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTc                              >  2:3364043/1‑93 (MQ=255)

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

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

                                                             aGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCa                            <  2:2434718/93‑1 (MQ=255)

                                                             aGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCa                            >  2:6156417/1‑93 (MQ=255)

                                                             aGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCa                            <  1:4449060/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                        caAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAcc                 >  1:153524/1‑93 (MQ=255)

                                                                        caAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACAcc                 >  2:6548672/1‑93 (MQ=255)

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

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

                                                                          aaaTCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCg               <  2:3468592/93‑1 (MQ=255)

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

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

                                                                            atcatcGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGa             <  2:7327496/93‑1 (MQ=255)

                                                                             tcatcGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGAt            <  2:5954331/93‑1 (MQ=255)

                                                                              catcGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATa           <  1:1949474/93‑1 (MQ=255)

                                                                              catcGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATa           <  2:5168955/93‑1 (MQ=255)

                                                                              catcGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATa           <  1:2104042/93‑1 (MQ=255)

                                                                                  gTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcat       <  2:7216696/93‑1 (MQ=255)

                                                                                  gTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcat       >  2:1353655/1‑93 (MQ=255)

                                                                                    gAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatca     <  1:3978378/93‑1 (MQ=255)

                                                                                    gAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatca     <  2:5403383/93‑1 (MQ=255)

                                                                                    gAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatca     <  2:4248318/93‑1 (MQ=255)

AGAAGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGG  >  NC_000913/4039867‑4040047

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