BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,871,084	A→G	T234A (ACT→GCT)	cdgI →	putative c‑di‑GMP binding protein CdgI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,871,084	0	A	G	100.0%	1610.8 / NA	375	T234A (ACT→GCT)	cdgI	putative c‑di‑GMP binding protein CdgI
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (164/211); total (164/211)

Only 100 of 381 total aligned reads displayed.

TGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACT  >  NC_000913/1870986‑1871181

 gATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTg                                                                                                <  1:3426024/101‑1 (MQ=255)

  aTCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCAtgctgc                                                                                               >  1:3426877/1‑101 (MQ=255)

  aTCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCAtgctgc                                                                                               <  2:3426388/101‑1 (MQ=255)

  aTCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCAtgctgc                                                                                               <  2:3426033/101‑1 (MQ=255)

   tCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCg                                                                                              >  2:3426881/1‑101 (MQ=255)

     ccTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGAc                                                                                            >  2:3426886/1‑101 (MQ=255)

     ccTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGAc                                                                                            <  2:3426525/101‑1 (MQ=255)

      cTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTTTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGAcc                                                                                           >  2:3426892/1‑101 (MQ=255)

      cTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGAcc                                                                                           >  1:3426891/1‑101 (MQ=255)

      cTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGAcc                                                                                           <  2:3426468/101‑1 (MQ=255)

       ttttttAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCt                                                                                          <  2:3426352/101‑1 (MQ=255)

        tttttAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTg                                                                                         <  1:3426377/101‑1 (MQ=255)

        tttttAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTg                                                                                         >  1:3426896/1‑101 (MQ=255)

        tttttAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTg                                                                                         <  1:3426155/101‑1 (MQ=255)

         ttttAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTgg                                                                                        <  1:3426287/101‑1 (MQ=255)

             aGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTg                                                                                                >  2:3426902/1‑89 (MQ=255)

             aGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGaa                                                                                    <  1:3425947/101‑1 (MQ=255)

             aGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGaa                                                                                    <  1:3425853/101‑1 (MQ=255)

              gTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAAt                                                                                   >  2:3426905/1‑101 (MQ=255)

              gTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAAt                                                                                   <  1:3426478/101‑1 (MQ=255)

               tttttCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATc                                                                                  <  1:3426486/101‑1 (MQ=255)

                ttttCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCa                                                                                 <  2:3426505/101‑1 (MQ=255)

                ttttCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCa                                                                                 >  1:3426907/1‑101 (MQ=255)

                 tttCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCaa                                                                                >  2:3426911/1‑101 (MQ=255)

                    ccccTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACta                                                                             >  2:3426916/1‑101 (MQ=255)

                     cccTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACtat                                                                            >  1:3426919/1‑101 (MQ=255)

                      ccTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACtata                                                                           <  2:3426337/101‑1 (MQ=255)

                      ccTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACtata                                                                           <  2:3426183/101‑1 (MQ=255)

                       cTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATAc                                                                          <  2:3426291/101‑1 (MQ=255)

                        ttCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACa                                                                         <  2:3426090/101‑1 (MQ=255)

                        ttCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACa                                                                         <  2:3426274/101‑1 (MQ=255)

                          cTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACaaa                                                                       >  2:3426924/1‑101 (MQ=255)

                             gCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAat                                                                    <  1:3426644/101‑1 (MQ=255)

                              cacacaATCGGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAata                                                                   >  1:3426929/1‑101 (MQ=255)

                               acacaATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAatat                                                                  <  1:3426345/101‑1 (MQ=255)

                               acacaATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAatat                                                                  >  2:3426931/1‑101 (MQ=255)

                               acacaATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAatat                                                                  <  1:3426088/101‑1 (MQ=255)

                                cacaATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATAtt                                                                 <  1:3426421/101‑1 (MQ=255)

                                  caATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTgg                                                               >  1:3426938/1‑101 (MQ=255)

                                  caATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTgg                                                               >  2:3426937/1‑101 (MQ=255)

                                    aTCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGtt                                                             >  2:3426942/1‑101 (MQ=255)

                                     tCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGttt                                                            >  1:3426944/1‑101 (MQ=255)

                                       tGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAAt                                                                                   >  2:3426950/1‑76 (MQ=255)

                                       tGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGt                                                          <  1:3426353/101‑1 (MQ=255)

                                         aGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCt                                                        >  1:3426953/1‑101 (MQ=255)

                                           cAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCtgt                                                      <  2:3426370/101‑1 (MQ=255)

                                            aGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCtgtg                                                     >  2:3426957/1‑101 (MQ=255)

                                             gTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTgg                                                    <  2:3426565/101‑1 (MQ=255)

                                              ttATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTggg                                                   <  2:3425978/101‑1 (MQ=255)

                                               tataACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGc                                                  <  1:3426596/101‑1 (MQ=255)

                                                 taACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCAt                                                >  2:3426963/1‑101 (MQ=255)

                                                 taACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCAt                                                <  2:3426323/101‑1 (MQ=255)

                                                  aaggCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCAtt                                               >  1:3426966/4‑101 (MQ=255)

                                                    cGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCAtttt                                             <  1:3426385/101‑1 (MQ=255)

                                                     gCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCAttttt                                            <  1:3426454/101‑1 (MQ=255)

                                                      cTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTa                                           <  2:3426435/101‑1 (MQ=255)

                                                       taaCTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGAc                                                                                            <  2:3426208/49‑1 (MQ=255)

                                                       tGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTAc                                          >  2:3426968/1‑101 (MQ=255)

                                                       tGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTAc                                          <  1:3426218/101‑1 (MQ=255)

                                                       tGACTATTCTTGGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTAc                                          >  2:3426972/1‑101 (MQ=255)

                                                        gACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACt                                         <  2:3426450/101‑1 (MQ=255)

                                                           tATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGtt                                      <  1:3426066/101‑1 (MQ=255)

                                                           tATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGtt                                      >  1:3426976/1‑101 (MQ=255)

                                                           tATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGtt                                      <  2:3426623/101‑1 (MQ=255)

                                                            aTTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTc                                     >  2:3426979/1‑101 (MQ=255)

                                                             ttCTTTGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCt                                    <  1:3426326/101‑1 (MQ=255)

                                                                tttGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTtta                                 >  2:3426984/1‑101 (MQ=255)

                                                                tttGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTtta                                 <  2:3426256/101‑1 (MQ=255)

                                                                 ttGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTttat                                <  2:3426419/101‑1 (MQ=255)

                                                                  tGTATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTttatt                               >  1:3426987/1‑101 (MQ=255)

                                                                    tATGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTtattat                             <  2:3426591/101‑1 (MQ=255)

                                                                     aTGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATc                            >  1:3426994/1‑101 (MQ=255)

                                                                      tGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCa                           >  2:3426996/1‑101 (MQ=255)

                                                                      tGTCGTCGATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCa                           <  1:3426248/101‑1 (MQ=255)

                                                                        tcgtcgATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATg                         <  1:3426306/101‑1 (MQ=255)

                                                                         cgtcgATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATgc                        >  1:3427002/1‑101 (MQ=255)

                                                                         cgtcgATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATgc                        >  2:3427001/1‑101 (MQ=255)

                                                                          gtcgATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATgcg                       <  1:3426400/101‑1 (MQ=255)

                                                                          gtcgATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATgcg                       >  2:3427006/1‑101 (MQ=255)

                                                                           tcgATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATgcgc                      <  2:3426485/101‑1 (MQ=255)

                                                                             gATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATgcg                       <  1:3426754/98‑1 (MQ=255)

                                                                             gATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGcac                    <  1:3426136/101‑1 (MQ=255)

                                                                             gATTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGcac                    >  2:3427012/1‑101 (MQ=255)

                                                                              aTTACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGcaca                   >  1:3427013/1‑101 (MQ=255)

                                                                               ttACTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCacac                  >  1:3427019/1‑101 (MQ=255)

                                                                                 aCTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGcac                    <  1:3426204/97‑1 (MQ=255)

                                                                                 aCTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGa                >  2:3427023/1‑101 (MQ=255)

                                                                                  cTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGAt               <  1:3426275/101‑1 (MQ=255)

                                                                                  cTGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGAt               <  2:3426651/101‑1 (MQ=255)

                                                                                   tGTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGAtt              <  1:3426362/101‑1 (MQ=255)

                                                                                    gTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTa             >  1:3427028/1‑101 (MQ=255)

                                                                                    gTCCGGACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTa             <  2:3426539/101‑1 (MQ=255)

                                                                                        ggACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCa         >  2:3427032/1‑101 (MQ=255)

                                                                                        ggACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCa         <  1:3426798/101‑1 (MQ=255)

                                                                                        ggACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCa         <  1:3426496/101‑1 (MQ=255)

                                                                                         gACAACCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAg        >  1:3427034/1‑101 (MQ=255)

                                                                                            aaCCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGa     >  1:3427039/1‑101 (MQ=255)

                                                                                            aaCCATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGa     >  2:3427040/1‑101 (MQ=255)

                                                                                              ccATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAAc   <  1:3426537/101‑1 (MQ=255)

                                                                                              ccATGCTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAAc   <  2:3426428/101‑1 (MQ=255)

TGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGACCTGGGGAATCAACTATACAAAAATATTGGTTTGTCTGTGGGCATTTTTACTGTTCTTTATTATCATGCGCACACGATTAGCCAGCGAACT  >  NC_000913/1870986‑1871181

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: `‑`