BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,934,276	C→T	H125Y (CAC→TAC)	rbsB →	D‑ribose periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,934,276	0	C	T	100.0%	969.4 / NA	259	H125Y (CAC→TAC)	rbsB	D‑ribose periplasmic binding protein
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (119/140); total (119/140)

Only 100 of 265 total aligned reads displayed.

GCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTG  >  NC_000913/3934186‑3934365

   gTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACAtt                                                                                      >  1:863647/1‑93 (MQ=255)

    tGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTg                                                                                     >  1:1534695/1‑93 (MQ=255)

      ggTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCt                                                                                   <  1:2089005/93‑1 (MQ=255)

        tAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTc                                                                                 >  1:4287366/1‑93 (MQ=255)

        tAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTc                                                                                 >  1:2008397/1‑93 (MQ=255)

        tAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTc                                                                                 <  2:1299633/93‑1 (MQ=255)

         aaTGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCt                                                                                >  1:2794200/1‑93 (MQ=255)

             cTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATa                                                                            <  2:669139/93‑1 (MQ=255)

             cTGTGAAGATGGCTAACCAGGCGAACATCCCGGTGATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATa                                                                            <  2:151549/93‑1 (MQ=255)

              tgtgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATaa                                                                           >  1:5261363/1‑93 (MQ=255)

               gtgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGGTTCTGATAAc                                                                          <  1:2866380/93‑1 (MQ=255)

               gtgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAAc                                                                          <  1:1644345/93‑1 (MQ=255)

                 gAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGt                                                                        <  2:6196191/93‑1 (MQ=255)

                 gAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGt                                                                        >  2:3821078/1‑93 (MQ=255)

                  aaGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTa                                                                       <  2:5457595/93‑1 (MQ=255)

                   aGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTAc                                                                      <  1:122158/93‑1 (MQ=255)

                   aGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTAc                                                                      <  2:3548614/93‑1 (MQ=255)

                    gATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACt                                                                     >  1:3715227/1‑93 (MQ=255)

                      tGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTgg                                                                   <  1:4782917/93‑1 (MQ=255)

                          tAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGggcgg                                                               >  1:560701/1‑93 (MQ=255)

                           aaCCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGggcggc                                                              >  1:503188/1‑93 (MQ=255)

                              cAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAATGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCaaa                                                           >  2:5133875/1‑93 (MQ=255)

                              cAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCaaa                                                           >  1:5720781/1‑93 (MQ=255)

                              cAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCaaa                                                           >  2:3697573/1‑93 (MQ=255)

                               aGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCaaaa                                                          >  1:1719549/1‑93 (MQ=255)

                               aGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCaaaa                                                          <  1:5781941/93‑1 (MQ=255)

                               aGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCaaaa                                                          <  2:4490191/93‑1 (MQ=255)

                                   gAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGc                                                      >  2:4936118/1‑93 (MQ=255)

                                   gAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGc                                                      <  1:2594294/93‑1 (MQ=255)

                                    aaCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCt                                                     >  1:4218602/1‑93 (MQ=255)

                                      cATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTgg                                                   >  1:328737/1‑93 (MQ=255)

                                      cATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTgg                                                   >  1:1393416/1‑93 (MQ=255)

                                        tCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTg                                                 >  1:1107506/1‑93 (MQ=255)

                                         cccGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGa                                                <  2:2100029/93‑1 (MQ=255)

                                         cccGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGa                                                >  2:556914/1‑93 (MQ=255)

                                         cccGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGa                                                <  1:5433614/93‑1 (MQ=255)

                                         cccGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGa                                                >  1:3133544/1‑93 (MQ=255)

                                          ccGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGAt                                               >  2:2488541/1‑93 (MQ=255)

                                          ccGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGAt                                               >  1:1510342/1‑93 (MQ=255)

                                             gTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTAc                                            >  2:3618528/1‑93 (MQ=255)

                                              ttATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTGCTGGGCGGCAAAATCGCTGGTGATTACa                                           >  1:3755434/1‑93 (MQ=255)

                                              ttATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACa                                           <  2:2428606/93‑1 (MQ=255)

                                              ttATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACa                                           >  1:3656663/1‑93 (MQ=255)

                                               tATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACAt                                          <  2:3373047/93‑1 (MQ=255)

                                               tATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACAt                                          <  2:3970846/93‑1 (MQ=255)

                                               tATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACAt                                          >  1:5231908/1‑93 (MQ=255)

                                                ttCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTTGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATc                                         <  2:791933/92‑1 (MQ=255)

                                                aTCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATc                                         <  2:1841114/93‑1 (MQ=255)

                                                  cgctctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATcgc                                       <  2:4781119/91‑1 (MQ=255)

                                                  cACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATcgc                                       <  1:603465/93‑1 (MQ=255)

                                                   aCTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATcgcg                                      >  1:5520713/1‑93 (MQ=255)

                                                       ttGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCgaaga                                  >  2:3329922/1‑93 (MQ=255)

                                                        tGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCgaagaa                                 <  1:2384560/93‑1 (MQ=255)

                                                        tGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCgaagaa                                 <  2:6070358/93‑1 (MQ=255)

                                                         gACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGaaa                                <  1:4115971/93‑1 (MQ=255)

                                                          aCCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAg                               >  1:1336919/1‑93 (MQ=255)

                                                           ccGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGc                              <  1:5713014/93‑1 (MQ=255)

                                                           ccGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGc                              <  1:5170028/93‑1 (MQ=255)

                                                              ccAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCggg                           <  2:4592150/93‑1 (MQ=255)

                                                               cAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCgggg                          >  1:1820889/1‑92 (MQ=255)

                                                                aGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTg                         <  2:5958089/93‑1 (MQ=255)

                                                                aGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTg                         >  1:2104476/1‑93 (MQ=255)

                                                                 ggCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGa                        >  2:2736052/1‑93 (MQ=255)

                                                                   cAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAg                      <  2:5762116/93‑1 (MQ=255)

                                                                    aaCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAgg                     >  1:3374398/1‑93 (MQ=255)

                                                                     aCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGt                    <  1:4186613/93‑1 (MQ=255)

                                                                     aCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGt                    <  1:362472/93‑1 (MQ=255)

                                                                     aCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGt                    >  1:1495007/1‑93 (MQ=255)

                                                                     aCGAAAGGTGAAGTGGTGAGCTACATTACTTCTGATACCGTACTGGGCGCCAAAATCGCGCGTGATTCCATCGCGAAGAAAGCGGGTGAAGGt                    <  2:2052105/93‑1 (MQ=255)

                                                                      cGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTg                   <  2:1063444/93‑1 (MQ=255)

                                                                      cGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTg                   >  1:3335447/1‑93 (MQ=255)

                                                                       gAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGc                  >  1:2195689/1‑93 (MQ=255)

                                                                       gAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGc                  <  1:643041/93‑1 (MQ=255)

                                                                         aaGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCa                >  2:2580583/1‑93 (MQ=255)

                                                                          aGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCaa               >  2:4700705/1‑93 (MQ=255)

                                                                          aGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCaa               >  2:5598689/1‑93 (MQ=255)

                                                                          aGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCaa               <  2:4321655/93‑1 (MQ=255)

                                                                           ggTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCaaa              <  2:4439560/93‑1 (MQ=255)

                                                                           ggTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCaaa              >  1:900416/1‑93 (MQ=255)

                                                                           ggTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCaaa              >  1:2444041/1‑93 (MQ=255)

                                                                             tGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGt            >  2:1365610/1‑93 (MQ=255)

                                                                             tGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGt            >  2:1601242/1‑93 (MQ=255)

                                                                              gAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGtt           >  2:738876/1‑93 (MQ=255)

                                                                              gAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGtt           <  2:3241191/93‑1 (MQ=255)

                                                                               aaGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTa          >  2:343732/1‑93 (MQ=255)

                                                                               aaGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTa          <  1:39933/93‑1 (MQ=255)

                                                                               aaGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTa          <  1:4716892/93‑1 (MQ=255)

                                                                                aGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAACATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTAt         >  1:4434467/1‑93 (MQ=255)

                                                                                aGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTAt         <  2:3796176/93‑1 (MQ=255)

                                                                                 gtggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATc        <  2:2971772/93‑1 (MQ=255)

                                                                                 gtggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATc        <  1:5337658/93‑1 (MQ=255)

                                                                                 gtggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATc        >  2:1185358/1‑93 (MQ=255)

                                                                                 gtggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATc        <  1:6257531/93‑1 (MQ=255)

                                                                                  tggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCg       >  1:6187191/1‑93 (MQ=255)

                                                                                   ggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGa      <  2:1919240/93‑1 (MQ=255)

                                                                                     tgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGc    <  1:2922304/93‑1 (MQ=255)

                                                                                      gAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCt   <  2:2231419/93‑1 (MQ=255)

                                                                                       aGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTg  >  2:4232853/1‑93 (MQ=255)

                                                                                       aGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTg  <  2:52544/93‑1 (MQ=255)

                                                                                       aGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTg  >  2:170249/1‑93 (MQ=255)

GCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTG  >  NC_000913/3934186‑3934365

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