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

Only 100 of 287 total aligned reads displayed.

AGAAGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGG  >  NC_000913/4039867‑4040047

agaagaCTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATa                                                                                          <  1:1303609/93‑1 (MQ=255)

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

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

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

  aagaCTGCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATatc                                                                                        <  1:1768978/93‑3 (MQ=255)

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

   agaATTCAAAACGCCGCTGTTACGCGCTCGCGGCTTAATTAAAGAGTAAAAGCTTTTAAGCGGCGCCACCAAAATCATCGTGAAATGATatcc                                                                                       <  1:2272765/93‑4 (MQ=255)

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

       ttCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCg                                                                                  <  2:4494070/93‑1 (MQ=255)

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

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

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

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

            aaCGGCGCTGTTACGCGCTCGCGGCGTAGTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCAtt                                                                             >  1:821454/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

                    tGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCCGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGt                                                                     >  1:3479861/1‑93 (MQ=255)

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

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

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

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

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

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

                               cgcgGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGgatga                                                          <  1:621539/93‑1 (MQ=255)

                               cgcgGCGTAAGTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGgatga                                                          <  1:1895238/93‑1 (MQ=255)

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

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

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

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

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

                                   gcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATggg                                                      >  1:582236/1‑93 (MQ=255)

                                   gcgTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATggg                                                      >  1:1575979/1‑93 (MQ=255)

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

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

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

                                       aaTTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGa                                                  <  1:1494510/93‑1 (MQ=255)

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

                                         ttAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAAt                                                >  2:1694653/1‑93 (MQ=255)

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

                                         ttAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAAt                                                >  1:3232525/1‑93 (MQ=255)

                                          tAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATa                                               <  2:3551227/93‑1 (MQ=255)

                                           aaaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATaa                                              <  1:3768500/93‑1 (MQ=255)

                                           aaaGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATaa                                              <  2:910906/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                      gCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTa                                   <  2:1475391/93‑1 (MQ=255)

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

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

                                                        ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACt                                 >  1:4288381/1‑93 (MQ=255)

                                                        ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACt                                 >  2:735634/1‑93 (MQ=255)

                                                        ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACt                                 >  1:1423031/1‑93 (MQ=255)

                                                        ttGTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACt                                 >  1:2072856/1‑93 (MQ=255)

                                                          gTAAGCGGCGCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACttt                               >  1:3323246/1‑93 (MQ=255)

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

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

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

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

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

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

                                                                ggcgCCACCAAAATCATCGTGAAATGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGac                         <  1:2326762/93‑1 (MQ=255)

                                                                  cgcCACCAAAATCATCGTGAAATGAT‑ATCATTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGacac                       <  1:3080113/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                     aaaTGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACcatcat    <  1:3374290/93‑1 (MQ=255)

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

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

                                                                                       aTGAT‑ATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATgg  <  2:3406765/93‑1 (MQ=255)

AGAAGACTTCAAAACGGCGCTGTTACGCGCTCGCGGCGTAATTAAAGAGTAAAAGCTTGTAAGCGGCGCCACCAAAATCATCGTGAAATGATAATCCTTCGTCATTCGTAATGTTTTCCGGATGATGGGATGAATAACAGCGCTTTTACTTTCCAGACACTACACCCGGATACCATCATGG  >  NC_000913/4039867‑4040047

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