BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,152,410	C→T	R59C (CGC→TGC)	metC →	cystathionine beta‑lyase/L‑cysteine desulfhydrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,152,410	0	C	T	100.0%	1598.0 / NA	385	R59C (CGC→TGC)	metC	cystathionine beta‑lyase/L‑cysteine desulfhydrase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (165/220); total (165/220)

Only 100 of 391 total aligned reads displayed.

GCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGC  >  NC_000913/3152312‑3152507

gCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGgtgc                                                                                                 <  2:5720255/101‑1 (MQ=255)

gCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGgtgc                                                                                                 >  1:5720907/1‑101 (MQ=255)

  gTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCgg                                                                                               >  1:5720912/1‑101 (MQ=255)

    gATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGaa                                                                                             <  1:5720292/101‑1 (MQ=255)

    gATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGaa                                                                                             >  2:5720915/1‑101 (MQ=255)

     aTTCACCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAGCACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAAc                                                                                            <  1:5720398/101‑1 (MQ=255)

      ttCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACg                                                                                           <  1:5720539/101‑1 (MQ=255)

       tCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGt                                                                                          <  2:5720182/101‑1 (MQ=255)

        cAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGtt                                                                                         <  2:5720336/101‑1 (MQ=255)

        cAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGtt                                                                                         <  1:5720326/101‑1 (MQ=255)

         aGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTa                                                                                        <  2:5720151/101‑1 (MQ=255)

         aGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTa                                                                                        >  2:5720922/1‑101 (MQ=255)

         aGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTa                                                                                        >  1:5720923/1‑101 (MQ=255)

          gcgcgcTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTaa                                                                                       <  1:5720482/101‑1 (MQ=255)

            gcgcTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAAcc                                                                                     <  1:5720505/101‑1 (MQ=255)

             cgcTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAAccc                                                                                    >  2:5720930/1‑101 (MQ=255)

              gcTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCa                                                                                   >  1:5720931/1‑101 (MQ=255)

              gcTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCa                                                                                   <  1:5720254/101‑1 (MQ=255)

               cttcttCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCAt                                                                                  <  1:5720586/101‑1 (MQ=255)

               cttcttCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCAt                                                                                  >  1:5720936/1‑101 (MQ=255)

               cttcttCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCAt                                                                                  <  2:5720310/101‑1 (MQ=255)

                ttcttcGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCAtt                                                                                 >  2:5720940/1‑101 (MQ=255)

                 tcttcGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCAttt                                                                                <  2:5719969/101‑1 (MQ=255)

                     cGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTcc                                                                            <  1:5720347/101‑1 (MQ=255)

                     cGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTcc                                                                            >  1:5720946/1‑101 (MQ=255)

                      gCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCt                                                                           >  2:5720947/1‑101 (MQ=255)

                       cTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCtt                                                                          <  2:5720322/101‑1 (MQ=255)

                       cTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCtt                                                                          <  2:5720522/101‑1 (MQ=255)

                         ggTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTac                                                                        <  1:5720384/101‑1 (MQ=255)

                          gTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTaca                                                                       >  1:5720955/1‑101 (MQ=255)

                             tttGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTacaaca                                                                    >  1:5720958/1‑101 (MQ=255)

                               tGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAg                                                                  >  2:5720961/1‑101 (MQ=255)

                                gACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTa                                                                         <  1:5720427/93‑1 (MQ=255)

                                gACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGc                                                                 >  1:5720963/1‑101 (MQ=255)

                                  cAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGa                                                               <  2:5720503/101‑1 (MQ=255)

                                   aGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGAt                                                              <  1:5720194/101‑1 (MQ=255)

                                   aGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGAt                                                              <  1:5720523/101‑1 (MQ=255)

                                      gtAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGAtgtg                                                           >  1:5720970/1‑101 (MQ=255)

                                       tAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGAtgtgt                                                          >  2:5720973/1‑101 (MQ=255)

                                         gAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGa                                                        >  1:5720977/1‑101 (MQ=255)

                                          aaGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGaa                                                       <  2:5720242/101‑1 (MQ=255)

                                           aGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAAc                                                      >  2:5720979/1‑101 (MQ=255)

                                             ccAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGTGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTg                                                    >  2:5720986/1‑101 (MQ=255)

                                             ccAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTtc                                                                               <  2:5720220/74‑1 (MQ=255)

                                             ccAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTg                                                    >  2:5720983/1‑101 (MQ=255)

                                               aaaaaaCACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGAtgt                                                            <  1:5720641/91‑1 (MQ=255)

                                               aaaaaaCACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGAtg                                                             <  2:5720533/90‑1 (MQ=255)

                                               aaaaaaCACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGa                                                        <  1:5719985/95‑1 (MQ=255)

                                               aaaaaaCACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAAc                                                      <  2:5720380/97‑1 (MQ=255)

                                                aaaaaCACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGaa                                                 >  2:5720994/1‑101 (MQ=255)

                                                aaaaaCACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGaa                                                 >  1:5720992/1‑101 (MQ=255)

                                                  aaaCACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAgg                                               <  1:5720356/101‑1 (MQ=255)

                                                   aaCACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGt                                              >  1:5720999/1‑101 (MQ=255)

                                                   aaCACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGt                                              <  1:5720065/101‑1 (MQ=255)

                                                    acacGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTg                                             <  2:5720619/101‑1 (MQ=255)

                                                       cgcgACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGgcg                                          <  2:5720439/101‑1 (MQ=255)

                                                        gcgACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTtc                                                                               <  2:5720129/63‑1 (MQ=255)

                                                        gcgACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGgcgc                                         >  2:5721003/1‑101 (MQ=255)

                                                          gACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAg                                       <  2:5720073/101‑1 (MQ=255)

                                                          gACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAg                                       <  1:5720376/101‑1 (MQ=255)

                                                            cacGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGc                                     >  1:5721008/1‑101 (MQ=255)

                                                             acGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCt                                    <  1:5720337/101‑1 (MQ=255)

                                                              cGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTg                                   >  2:5721011/1‑101 (MQ=255)

                                                               gTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGc                                  <  1:5720456/101‑1 (MQ=255)

                                                                tAATCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCg                                 >  1:5721016/1‑101 (MQ=255)

                                                                  aTCGCGCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTg                               >  2:5721019/1‑101 (MQ=255)

                                                                    cgcgCCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCt                             <  1:5720765/101‑1 (MQ=255)

                                                                      cgcCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTAt                           <  1:5720444/101‑1 (MQ=255)

                                                                       gcCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTAtt                          <  2:5720702/101‑1 (MQ=255)

                                                                       gcCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTAtt                          >  1:5721023/1‑101 (MQ=255)

                                                                       gcCAATGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTAtt                          <  1:5720556/101‑1 (MQ=255)

                                                                         cAATGGAGAGTTGATCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGgc                                           <  2:5720651/82‑1 (MQ=255)

                                                                          aaTGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTcc                       >  2:5721029/1‑101 (MQ=255)

                                                                           attgagaGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTccc                      >  1:5721032/4‑101 (MQ=255)

                                                                           aTGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGt                                <  2:5720490/91‑1 (MQ=255)

                                                                           aTGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTccc                      <  2:5720417/101‑1 (MQ=255)

                                                                            tGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCt                     <  1:5720223/101‑1 (MQ=255)

                                                                            tGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCt                     >  2:5721034/1‑101 (MQ=255)

                                                                            tGGAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCt                     <  2:5720756/101‑1 (MQ=255)

                                                                             ggAGAGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTg                    <  1:5720369/101‑1 (MQ=255)

                                                                              gagaGTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGc                   >  1:5721039/1‑101 (MQ=255)

                                                                               agagTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCg                  <  2:5720835/101‑1 (MQ=255)

                                                                               agagTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCg                  <  2:5720365/101‑1 (MQ=255)

                                                                               agagTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCg                  >  2:5721043/1‑101 (MQ=255)

                                                                                gagTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCgg                 <  2:5720473/101‑1 (MQ=255)

                                                                                 agTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCggg                <  2:5720281/101‑1 (MQ=255)

                                                                                 agTTGTTCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCggg                >  1:5721046/1‑101 (MQ=255)

                                                                                    tgttCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGggcg             >  1:5721050/1‑101 (MQ=255)

                                                                                     gttCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGggcgg            <  2:5720598/101‑1 (MQ=255)

                                                                                      ttCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTgg                                            <  1:5720495/68‑1 (MQ=255)

                                                                                      ttCTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGggcggc           >  2:5721054/1‑101 (MQ=255)

                                                                                        cTATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAg         >  1:5721058/1‑101 (MQ=255)

                                                                                         tATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAg         <  2:5720548/100‑1 (MQ=255)

                                                                                         tATGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGc        <  2:5720458/101‑1 (MQ=255)

                                                                                          aTGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCg       >  2:5721062/1‑101 (MQ=255)

                                                                                           tGGACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCgg      >  1:5721065/1‑101 (MQ=255)

                                                                                             gACGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGtt    >  2:5721069/1‑101 (MQ=255)

                                                                                              aCGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTg   <  1:5720308/101‑1 (MQ=255)

                                                                                              aCGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTg   <  2:5720449/101‑1 (MQ=255)

                                                                                               cGGTGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGggc              <  2:5720404/89‑1 (MQ=255)

GCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGC  >  NC_000913/3152312‑3152507

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`