BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,295,417	G→A	P510L (CCG→CTG)	fdhF ←	selenopolypeptide subunit of formate dehydrogenase H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,295,417	0	G	A	100.0%	800.1 / NA	219	P510L (CCG→CTG)	fdhF	selenopolypeptide subunit of formate dehydrogenase H
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (119/100); total (119/100)

Only 100 of 227 total aligned reads displayed.

GTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATT  >  NC_000913/4295327‑4295506

 tCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGc                                                                                        >  2:3342423/1‑93 (MQ=255)

   cccTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGcac                                                                                      >  1:996086/1‑93 (MQ=255)

   cccTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGcac                                                                                      >  1:4325745/1‑93 (MQ=255)

    ccTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGCAATCCAGGcaca                                                                                     >  2:807480/1‑93 (MQ=255)

     cTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAg                                                                                    >  2:2270819/1‑93 (MQ=255)

     cTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAg                                                                                    >  2:3555062/1‑93 (MQ=255)

     cTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAg                                                                                    <  1:2034935/93‑1 (MQ=255)

     cTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAg                                                                                    >  1:574572/1‑93 (MQ=255)

      tGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGa                                                                                   <  1:2502175/93‑1 (MQ=255)

      tGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGa                                                                                   <  2:4496518/93‑1 (MQ=255)

      tGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGa                                                                                   >  1:4475781/1‑93 (MQ=255)

       gATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGAt                                                                                  >  1:4192770/1‑93 (MQ=255)

        aTCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATg                                                                                 >  2:1835089/1‑93 (MQ=255)

         tCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGa                                                                                <  1:3092442/93‑1 (MQ=255)

          cGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGAc                                                                               <  1:3932890/93‑1 (MQ=255)

          cGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGAc                                                                               <  2:2567847/93‑1 (MQ=255)

            gCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGc                                                                             >  2:1971526/1‑93 (MQ=255)

            gCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGc                                                                             >  2:2300388/1‑93 (MQ=255)

              aTCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCaa                                                                           >  1:1975140/1‑93 (MQ=255)

               tCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGGCGCAAc                                                                          >  2:3420620/1‑93 (MQ=255)

               tCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAAc                                                                          >  1:3246501/1‑93 (MQ=255)

               tCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAAc                                                                          >  1:1312415/1‑93 (MQ=255)

                 tGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTc                                                                        <  2:1315922/93‑1 (MQ=255)

                 tGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTc                                                                        >  1:50895/1‑93 (MQ=255)

                  gAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCa                                                                       >  1:3406401/1‑93 (MQ=255)

                    aGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATc                                                                     <  1:2617222/93‑1 (MQ=255)

                    aGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATc                                                                     >  2:3696937/1‑93 (MQ=255)

                      tATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATccc                                                                   >  2:2658183/1‑93 (MQ=255)

                       aTCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCa                                                                  >  1:1130716/1‑93 (MQ=255)

                          gcgGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAt                                                               <  1:94013/93‑1 (MQ=255)

                          gcgGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAt                                                               >  2:4405434/1‑93 (MQ=255)

                           cgGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAtc                                                              <  2:1237472/93‑1 (MQ=255)

                           cgGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATACCAGAtc                                                              >  1:1621948/1‑93 (MQ=255)

                            ggCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAtct                                                             >  2:3093620/1‑93 (MQ=255)

                             gCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAtctc                                                            <  2:878135/93‑1 (MQ=255)

                             gCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGAtctc                                                            >  2:736178/1‑93 (MQ=255)

                              cAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTcc                                                           >  1:3839417/1‑93 (MQ=255)

                                aGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTg                                                         >  1:2819364/1‑93 (MQ=255)

                                 ggCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTgg                                                        >  2:4551987/1‑93 (MQ=255)

                                 ggCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTgg                                                        >  1:3182744/1‑93 (MQ=255)

                                    cACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgtg                                                     >  1:3468867/1‑93 (MQ=255)

                                    cACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgtg                                                     <  1:823555/93‑1 (MQ=255)

                                     aCTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCATCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgtgt                                                    >  2:1448372/1‑93 (MQ=255)

                                     aCTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGgtgt                                                    >  2:433399/1‑93 (MQ=255)

                                        gaatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgt                                                 >  1:4404821/1‑93 (MQ=255)

                                        gaatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgt                                                 <  2:3440996/93‑1 (MQ=255)

                                        gaatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgt                                                 >  1:2732650/1‑93 (MQ=255)

                                         actaaaGCCCAGTTCGCCTATTTTCTCGTAAGTAGCACCGTAGAAATCCCGGCACAGATGACGCAACTCATCCCAGATCTCCTGGTcgttgtt                                                >  2:596885/5‑93 (MQ=255)

                                         aatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgtt                                                >  2:2840198/1‑93 (MQ=255)

                                         aatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgtt                                                >  2:2350345/1‑93 (MQ=255)

                                         aatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgtt                                                <  1:2599122/93‑1 (MQ=255)

                                         aatgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgtt                                                >  1:1818750/1‑93 (MQ=255)

                                           tgaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGtgttgttgt                                              <  1:4319940/93‑1 (MQ=255)

                                            gaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAGATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTa                                             <  2:4643693/93‑1 (MQ=255)

                                            gaaGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTa                                             <  1:3949378/93‑1 (MQ=255)

                                              aGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGt                                           <  2:2018841/93‑1 (MQ=255)

                                                  cAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCAt                                       <  2:3329418/93‑1 (MQ=255)

                                                   aGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATc                                      <  1:1340520/93‑1 (MQ=255)

                                                     ttCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCgg                                    >  1:1013577/1‑93 (MQ=255)

                                                      tCGCCCATTTTCTCGTAAGTCGCTCCGTAGAAATCCAGGCACAGATGACGCATCTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGa                                   <  2:4560874/93‑1 (MQ=255)

                                                      tCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGa                                   <  2:2159662/93‑1 (MQ=255)

                                                       cGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGAt                                  <  2:1497676/93‑1 (MQ=255)

                                                         cccATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATaa                                >  2:3306035/1‑93 (MQ=255)

                                                          ccATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAAc                               <  1:862196/93‑1 (MQ=255)

                                                           cATTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGCTCTCCTGGGTGTTGTTGTAGTGCATCGGATAAcc                              >  2:3022613/1‑93 (MQ=255)

                                                            aTTTTCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAAccc                             <  2:1685989/93‑1 (MQ=255)

                                                             ttttCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCa                            >  1:2351179/1‑93 (MQ=255)

                                                             ttttCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCa                            <  1:4599438/93‑1 (MQ=255)

                                                              tttCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCAt                           <  2:2942222/93‑1 (MQ=255)

                                                               ttCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCAtc                          >  1:3598491/1‑92 (MQ=255)

                                                               ttCTCGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATa                          >  1:2050561/1‑93 (MQ=255)

                                                                 ctcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACg                        <  2:359176/93‑1 (MQ=255)

                                                                  tcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACgg                       >  2:1723628/1‑93 (MQ=255)

                                                                  tcGTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACgg                       <  1:4111607/93‑1 (MQ=255)

                                                                    gTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACgggg                     >  1:4008917/1‑92 (MQ=255)

                                                                    gTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGt                     <  1:1863375/93‑1 (MQ=255)

                                                                    gTAAGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGt                     >  2:4213528/1‑93 (MQ=255)

                                                                       aGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGc                  <  1:1709544/93‑1 (MQ=255)

                                                                       aGTCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGc                  <  2:1957367/93‑1 (MQ=255)

                                                                         tCGCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGa                >  1:1391280/1‑93 (MQ=255)

                                                                           gCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGAtt              >  1:251863/1‑93 (MQ=255)

                                                                           gCACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGAGTCGAtt              >  1:1502860/1‑93 (MQ=255)

                                                                            cACCGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGAttt             >  2:2662235/1‑93 (MQ=255)

                                                                              ccGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCa           >  2:2626001/1‑93 (MQ=255)

                                                                               cGTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCAc          <  2:497932/93‑1 (MQ=255)

                                                                                gTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACt         <  1:2659499/93‑1 (MQ=255)

                                                                                gTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACt         >  1:2451363/1‑93 (MQ=255)

                                                                                gTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACt         >  2:3921031/1‑93 (MQ=255)

                                                                                gTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACt         >  2:3659739/1‑93 (MQ=255)

                                                                                gTAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACACATACGGGAGGGGATTTCatt         >  2:4009936/1‑91 (MQ=255)

                                                                                 tAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtg        >  2:4201741/1‑93 (MQ=255)

                                                                                 tAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtg        >  2:4118693/1‑93 (MQ=255)

                                                                                 tAGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtg        >  2:3252684/1‑93 (MQ=255)

                                                                                  aGAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtga       <  1:1118770/93‑1 (MQ=255)

                                                                                   gAAATCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtgat      >  2:640188/1‑93 (MQ=255)

                                                                                     aaTCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACtgatga    >  2:1694585/1‑93 (MQ=255)

                                                                                      aTCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTgatgat   >  2:3184157/1‑93 (MQ=255)

                                                                                      aTCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTgatgat   <  2:1531029/93‑1 (MQ=255)

                                                                                       tCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGAtt  >  2:3853266/1‑93 (MQ=255)

                                                                                       tCCAGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGAtt  <  1:1167789/93‑1 (MQ=255)

GTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATT  >  NC_000913/4295327‑4295506

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