BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	446,355	T→C	N192S (AAC→AGC)	cyoE ←	protoheme IX farnesyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	446,355	0	T	C	100.0%	577.5 / NA	181	N192S (AAC→AGC)	cyoE	protoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (95/86); total (95/86)

Only 100 of 190 total aligned reads displayed.

TCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCT  >  NC_000913/446266‑446444

tCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTg                                                                                        >  1:2496683/1‑93 (MQ=255)

tCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTg                                                                                        >  2:505950/1‑93 (MQ=255)

  aGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTgcc                                                                                      <  1:3741434/93‑1 (MQ=255)

  aGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTgcc                                                                                      <  1:1925637/93‑1 (MQ=255)

     gTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTgccgcc                                                                                   <  2:1376877/93‑1 (MQ=255)

     gTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTgccgcc                                                                                   >  1:4454243/1‑93 (MQ=255)

     gTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTgccgcc                                                                                   <  1:1849686/93‑1 (MQ=255)

     gTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTgccgcc                                                                                   >  2:4054617/1‑93 (MQ=255)

     gTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTgccgcc                                                                                   <  2:558985/93‑1 (MQ=255)

      tGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCt                                                                                  <  2:3409724/93‑1 (MQ=255)

      tGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCt                                                                                  >  1:574600/1‑93 (MQ=255)

       ggCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTg                                                                                 >  1:3705227/1‑93 (MQ=255)

        gCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTgg                                                                                <  1:2455965/93‑1 (MQ=255)

         cAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGt                                                                               <  1:2804261/93‑1 (MQ=255)

          aaCGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCTGAATGCTTGCCGCCTGGTa                                                                              <  2:4573269/93‑1 (MQ=255)

           aCGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTaa                                                                             >  1:1305089/1‑93 (MQ=255)

            cGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAAt                                                                            >  2:510499/1‑93 (MQ=255)

              gCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATcc                                                                          <  2:3755179/93‑1 (MQ=255)

              gCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATcc                                                                          >  2:3623835/1‑93 (MQ=255)

              gCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATcc                                                                          <  2:3751773/93‑1 (MQ=255)

                aaaGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCCATACCGGAATGCTTGCCGCCTGGTAATCAtt                                                                        >  2:1585236/1‑93 (MQ=255)

                aaaGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCtt                                                                        <  2:1548955/93‑1 (MQ=255)

                 aaGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTa                                                                       >  1:3047013/1‑93 (MQ=255)

                  aGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTaa                                                                      <  2:2158957/93‑1 (MQ=255)

                   ggCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTaaa                                                                     <  1:4224328/93‑1 (MQ=255)

                   ggCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTaaa                                                                     <  1:690611/93‑1 (MQ=255)

                   ggCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTaaa                                                                     <  2:1106967/93‑1 (MQ=255)

                      gatgatATACCGCGGGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCg                                                                  >  2:2870118/1‑93 (MQ=255)

                      gatgatATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCg                                                                  >  2:2047128/1‑93 (MQ=255)

                      gatgatATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCg                                                                  >  1:1601299/1‑93 (MQ=255)

                      gatgatATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCg                                                                  >  1:1295537/1‑93 (MQ=255)

                       atgatATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGTgg                                                                 >  2:2645111/1‑93 (MQ=255)

                       atgatATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCgg                                                                 >  1:819329/1‑93 (MQ=255)

                       atgatATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCgg                                                                 <  1:2725709/93‑1 (MQ=255)

                       atgatATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCgg                                                                 >  2:1597638/1‑93 (MQ=255)

                          atatACAGCGGGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGGCGCCTGGTAATCCTTCAACCGGaca                                                              >  2:1718685/1‑91 (MQ=255)

                           tataCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGaaaa                                                             <  2:2394921/93‑1 (MQ=255)

                            ataCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAAt                                                            <  2:2496683/93‑1 (MQ=255)

                             taCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATg                                                           >  1:1438296/1‑93 (MQ=255)

                             taCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATg                                                           <  1:1826806/93‑1 (MQ=255)

                             taCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACGGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATg                                                           <  1:1258009/93‑1 (MQ=255)

                              aCAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATgg                                                          <  2:399039/93‑1 (MQ=255)

                               cAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGc                                                         <  1:296712/93‑1 (MQ=255)

                               cAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGc                                                         >  2:4625267/1‑93 (MQ=255)

                                aGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATGCCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCg                                                        <  1:2245403/93‑1 (MQ=255)

                                 gCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGa                                                       <  1:2086951/93‑1 (MQ=255)

                                  cGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGAt                                                      <  1:2373162/93‑1 (MQ=255)

                                   gtgatgtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATg                                                     >  2:2348251/1‑93 (MQ=255)

                                   gtgatgtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATg                                                     <  1:3897645/93‑1 (MQ=255)

                                   gtgatgtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATg                                                     >  1:3071399/1‑93 (MQ=255)

                                    ggatgtgatTCCTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATgg                                                    <  1:3235246/92‑1 (MQ=255)

                                     gatgtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATgga                                                   >  1:341197/1‑92 (MQ=255)

                                        gtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATa                                                >  2:773286/1‑93 (MQ=255)

                                        gtgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATa                                                <  1:2200188/93‑1 (MQ=255)

                                         tgatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAg                                               <  1:1867455/93‑1 (MQ=255)

                                          gatTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAgg                                              >  2:1472682/1‑93 (MQ=255)

                                           atTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGa                                             <  2:4379866/93‑1 (MQ=255)

                                            ttcttcGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAg                                            >  1:3364425/1‑93 (MQ=255)

                                            ttcttcGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAg                                            <  2:4445719/93‑1 (MQ=255)

                                            ttcttcGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAg                                            >  2:1894741/1‑93 (MQ=255)

                                            ttcttcGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAg                                            <  1:2041302/93‑1 (MQ=255)

                                             tcttcGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGt                                           >  1:3634240/1‑93 (MQ=255)

                                               ttcGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGa                                         <  2:416810/93‑1 (MQ=255)

                                                tcGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAg                                        >  2:1282592/1‑93 (MQ=255)

                                                   ccACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCa                                     <  2:1028375/93‑1 (MQ=255)

                                                    cACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCAt                                    >  2:3721540/1‑93 (MQ=255)

                                                    cACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCAt                                    >  2:4274377/1‑93 (MQ=255)

                                                    cACCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCAt                                    >  1:2135249/1‑93 (MQ=255)

                                                     aCCGAAATGCCTTTTACCACTGGCAATCCCGGAATGCGTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATc                                   <  1:1029884/93‑1 (MQ=255)

                                                     aCCGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATc                                   >  1:1077004/1‑93 (MQ=255)

                                                      ccGAAATGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCt                                  <  1:2263652/93‑1 (MQ=255)

                                                         aaaTGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGcc                               >  1:342813/1‑93 (MQ=255)

                                                          aaTGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCca                              <  1:4335055/93‑1 (MQ=255)

                                                            tGCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCcaca                            <  2:1854387/93‑1 (MQ=255)

                                                             gCCTTTTACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAg                           >  1:2941452/1‑93 (MQ=255)

                                                                   tACCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGaa                     <  1:3464349/93‑1 (MQ=255)

                                                                    aCCACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAg                    >  2:3212961/1‑93 (MQ=255)

                                                                     ccACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGa                   <  2:4611484/93‑1 (MQ=255)

                                                                      cACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGAt                  <  2:1332921/93‑1 (MQ=255)

                                                                      cACTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGAt                  >  1:2002112/1‑93 (MQ=255)

                                                                       aCTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATa                 >  2:1261051/1‑93 (MQ=255)

                                                                        cTGGCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAg                <  1:1017314/93‑1 (MQ=255)

                                                                           gCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCca             <  2:974898/93‑1 (MQ=255)

                                                                           gCAATACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCca             >  2:834464/1‑93 (MQ=255)

                                                                             aaTACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCcagc           >  2:428796/1‑93 (MQ=255)

                                                                             aaTACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGAGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCcagc           >  1:857640/1‑93 (MQ=255)

                                                                               tACCGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCcagcag         >  2:4508073/1‑93 (MQ=255)

                                                                                 ccGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGa       <  1:2330947/93‑1 (MQ=255)

                                                                                 ccGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGa       <  1:16683/93‑1 (MQ=255)

                                                                                  cGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGAt      <  1:3314337/93‑1 (MQ=255)

                                                                                  cGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGAt      <  2:2238084/93‑1 (MQ=255)

                                                                                  cGGAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTCCCACCGGCTGAAGATAGCCAGCAGGAt      >  2:1237101/1‑93 (MQ=255)

                                                                                   ggAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATc     <  1:1135013/93‑1 (MQ=255)

                                                                                   ggAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATc     >  2:676928/1‑93 (MQ=255)

                                                                                   ggAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATc     <  2:3026171/93‑1 (MQ=255)

                                                                                    gAATGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCg    >  1:1841565/1‑93 (MQ=255)

                                                                                     aaTGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGTATCGc   <  2:332655/93‑1 (MQ=255)

                                                                                     aaTGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGc   <  2:1080547/93‑1 (MQ=255)

                                                                                     aaTGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGc   <  2:2758304/93‑1 (MQ=255)

                                                                                      aTGCTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCt  <  1:3608957/93‑1 (MQ=255)

TCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCT  >  NC_000913/446266‑446444

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