BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	963,286	A→G	H79R (CAC→CGC)	himD →	integration host factor beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	963,286	0	A	G	100.0%	494.1 / NA	154	H79R (CAC→CGC)	himD	integration host factor beta subunit
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (65/89); total (65/89)

Only 100 of 159 total aligned reads displayed.

GCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCG  >  NC_000913/963196‑963375

gcggTTTCTCTTTGCACTACCGCGCACGATGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCt                                                                                         <  1:2284761/90‑1 (MQ=255)

  aGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCttt                                                                                       <  2:1656670/93‑1 (MQ=255)

  aGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCttt                                                                                       >  2:1350442/1‑93 (MQ=255)

  aGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCttt                                                                                       <  1:2694656/93‑1 (MQ=255)

   gTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAATTAGAACGGGAGGGAAAAAACGTTCTTCGCTTTa                                                                                      >  2:3361921/1‑93 (MQ=255)

   gTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTa                                                                                      >  2:3168434/1‑93 (MQ=255)

   gTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTa                                                                                      >  2:2810531/1‑93 (MQ=255)

     ttCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTaaa                                                                                    >  1:2220704/1‑93 (MQ=255)

      tctcTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAACATACGTTCCTCGCTTTAAAc                                                                                   >  1:1361826/1‑93 (MQ=255)

       ttctTTGCACTCCCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAAcc                                                                                  <  2:2002492/92‑1 (MQ=255)

       ctctTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGTAGGAAAATACGTTCCTCGTTTTatagc                                                                                  >  2:3100727/1‑89 (MQ=255)

        tctTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCt                                                                                 <  2:771849/93‑1 (MQ=255)

         ctTTGCGCTGCCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTg                                                                                <  2:1492863/93‑1 (MQ=255)

         ctTTGCACTACCGCGCCCCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGACCTGGAAGGAAAATACGTTCCTCGCTTTAAACCTg                                                                                <  1:506616/93‑1 (MQ=255)

            tGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTa                                                                             <  2:2032830/93‑1 (MQ=255)

             gCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTaa                                                                            >  2:1425986/1‑93 (MQ=255)

             gCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTaa                                                                            >  1:1941968/1‑93 (MQ=255)

              cACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTaaa                                                                           <  1:1042596/93‑1 (MQ=255)

              cACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTaaa                                                                           <  2:3725784/93‑1 (MQ=255)

              cACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTaaa                                                                           <  1:4145602/93‑1 (MQ=255)

               aCTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAg                                                                          >  2:2614640/1‑93 (MQ=255)

                cTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGa                                                                         >  2:4670492/1‑93 (MQ=255)

                cTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGa                                                                         >  2:2519225/1‑93 (MQ=255)

                cTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGa                                                                         >  1:608666/1‑93 (MQ=255)

                  aCCGCGCACCACGTACCGGACGTAATCGGAAGACCGGCGAAAACGGCGAACCGGAAGGAAAATACGTTCCTAGCTGTAAACCTAGTAAAgtct                                                                       >  2:4317610/1‑90 (MQ=255)

                  aCCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAAc                                                                       <  2:1156928/93‑1 (MQ=255)

                  aCCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAAc                                                                       <  2:362051/93‑1 (MQ=255)

                   ccGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACt                                                                      <  2:1114879/93‑1 (MQ=255)

                    cgcgCACCACGTACCGGACGTACTCCGAAGGCTAGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTTAACCAGGTTAAGAACta                                                                     >  2:3872010/1‑92 (MQ=255)

                    cgcgCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTg                                                                     <  2:4471988/93‑1 (MQ=255)

                    cgcgCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTg                                                                     <  2:3928747/93‑1 (MQ=255)

                     gcgcACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTgc                                                                    <  2:2458556/93‑1 (MQ=255)

                      cgcACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTgcg                                                                   <  2:2712282/93‑1 (MQ=255)

                         accacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGa                                                                >  1:370058/1‑93 (MQ=255)

                         accacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGa                                                                >  1:1901541/1‑93 (MQ=255)

                         accacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGa                                                                >  2:2412285/1‑93 (MQ=255)

                         accacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGa                                                                <  1:1837218/93‑1 (MQ=255)

                          ccacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGAt                                                               <  2:4121465/93‑1 (MQ=255)

                          ccacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGAt                                                               >  1:2428653/1‑93 (MQ=255)

                          ccacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGAt                                                               <  2:405512/93‑1 (MQ=255)

                          ccacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGAt                                                               >  1:1658869/1‑93 (MQ=255)

                           cacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATc                                                              >  2:2551015/1‑93 (MQ=255)

                           cacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATc                                                              >  1:3807662/1‑93 (MQ=255)

                           cacGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATc                                                              <  2:4091744/93‑1 (MQ=255)

                              gTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATcgcg                                                           >  1:3313652/1‑93 (MQ=255)

                               tACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCgcgc                                                          >  1:84244/1‑93 (MQ=255)

                                 ccGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGGGATCGCGCCa                                                        >  1:3073128/1‑93 (MQ=255)

                                  cGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCaa                                                       >  1:3040684/1‑93 (MQ=255)

                                  cGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCaa                                                       <  1:3089859/93‑1 (MQ=255)

                                   ggACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAat                                                      >  2:3015635/1‑93 (MQ=255)

                                   ggACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAat                                                      <  1:3575686/93‑1 (MQ=255)

                                    gACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAata                                                     <  2:2766640/93‑1 (MQ=255)

                                     aCGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAatat                                                    <  1:2398617/93‑1 (MQ=255)

                                      cGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATAtt                                                   <  2:3116936/93‑1 (MQ=255)

                                      cGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATAtt                                                   <  1:2081553/93‑1 (MQ=255)

                                         aaTCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTAc                                                <  2:1647270/93‑1 (MQ=255)

                                          aTCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACg                                               >  1:2087765/1‑93 (MQ=255)

                                           tCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACgg                                              >  1:4460225/1‑93 (MQ=255)

                                            ccGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGt                                             <  1:2705819/93‑1 (MQ=255)

                                             cGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGtt                                            <  2:2515348/93‑1 (MQ=255)

                                                 gACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGt                                        >  2:3545516/1‑93 (MQ=255)

                                                  aCTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGtt                                       <  2:1775468/93‑1 (MQ=255)

                                                  aCTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGtt                                       >  1:4369301/1‑93 (MQ=255)

                                                   cTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGGGCGATCGCGCCAATATTTACGGTTAAGttt                                      >  2:3203775/1‑93 (MQ=255)

                                                   cTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGttt                                      >  2:3856571/1‑93 (MQ=255)

                                                   cTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGttt                                      <  1:1338218/93‑1 (MQ=255)

                                                    tGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGtttt                                     >  1:1251013/1‑93 (MQ=255)

                                                     ggCGATAAAGTAGATCTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGttttt                                    <  1:151901/93‑1 (MQ=255)

                                                      gagATAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGtttttt                                   <  1:950928/91‑1 (MQ=255)

                                                         aTAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACt                                <  2:2012495/93‑1 (MQ=255)

                                                          tAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTc                               >  2:326862/1‑93 (MQ=255)

                                                          tAAAGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTc                               <  1:1655765/93‑1 (MQ=255)

                                                            aaGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCaa                             >  2:1676288/1‑93 (MQ=255)

                                                             aGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCaaa                            <  1:4373588/93‑1 (MQ=255)

                                                             aGTAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCaaa                            <  1:1112796/93‑1 (MQ=255)

                                                               tAGAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACt                          <  2:523033/93‑1 (MQ=255)

                                                                 gAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTg                        >  2:3740805/1‑93 (MQ=255)

                                                                 gAACTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTg                        >  1:2668717/1‑93 (MQ=255)

                                                                  aaCTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGa                       <  1:1753719/93‑1 (MQ=255)

                                                                   aCTGGAAGGAAAATGCGTTCCTCGCTTTAGACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGaa                      <  1:1083731/93‑1 (MQ=255)

                                                                   aCTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGaa                      >  2:2825353/1‑93 (MQ=255)

                                                                    cTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAAc                     >  2:4145669/1‑93 (MQ=255)

                                                                    cTGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAAc                     <  2:1266963/93‑1 (MQ=255)

                                                                     tGGAAGGAAAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACg                    >  2:2597541/1‑93 (MQ=255)

                                                                      ggaaggaaAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACga                   <  1:4638828/93‑1 (MQ=255)

                                                                      ggaaggaaAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACga                   >  2:2248019/1‑93 (MQ=255)

                                                                      ggaaggaaAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACga                   <  1:3515985/93‑1 (MQ=255)

                                                                        aaggaaAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACgaga                 <  1:1754418/93‑1 (MQ=255)

                                                                        aaggaaAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACgaga                 >  1:3991990/1‑93 (MQ=255)

                                                                          ggaaAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACgagaga               <  2:1344409/93‑1 (MQ=255)

                                                                          ggaaAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACgagaga               >  1:2899673/1‑93 (MQ=255)

                                                                          ggaaAATACGTTCCTCGCTTTAAACCTGGCAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACgagaga               >  1:605920/1‑93 (MQ=255)

                                                                           gaaAATACGTTCCTCGCTTTAAACCTGGTAAAGAACTGACCGATAGCGCCAATATTTACGGTTAAGGTTTTTACTCAAACTTGAACGAGAGaa              >  2:2987603/1‑93 (MQ=255)

                                                                              aaTACGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAg           <  1:2145691/93‑1 (MQ=255)

                                                                               aTACGTTCCTCGCTTTAATCCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGc          <  2:1499940/93‑1 (MQ=255)

                                                                                  cGTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCAcc       <  2:4190393/93‑1 (MQ=255)

                                                                                   gTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCt      >  1:4134449/1‑93 (MQ=255)

                                                                                   gTTCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCt      <  2:187186/93‑1 (MQ=255)

                                                                                     tCCTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGt    >  2:2727773/1‑93 (MQ=255)

                                                                                      ccTCGCTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTc   <  2:2342493/93‑1 (MQ=255)

GCAGTTTCTCTTTGCACTACCGCGCACCACGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAAAGCACCTGTCG  >  NC_000913/963196‑963375

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