BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,448,541	+C	intergenic (+46/‑92)	ytfQ → / → ytfR	putative LACI‑type transcriptional regulator/YtfR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,448,541	1	.	C	100.0%	343.7 / NA	118	intergenic (+46/‑92)	ytfQ/ytfR	putative LACI‑type transcriptional regulator/YtfR
Reads supporting (aligned to +/- strand): ref base . (0/0); new base C (65/53); total (65/53)

Only 100 of 132 total aligned reads displayed.

TACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCG‑G‑C‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCG  >  NC_000913/4448450‑4448623

tACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgc                                                                                       <  1:2856142/93‑2 (MQ=255)

tACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgc                                                                                       <  1:3714334/93‑2 (MQ=255)

tACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgc                                                                                       >  2:3785179/1‑92 (MQ=255)

 aCCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcg                                                                                      >  2:488097/1‑91 (MQ=255)

   cGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcgcc                                                                                    <  2:1868052/93‑5 (MQ=255)

           aagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGTTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCgg                                                                          <  1:4487507/93‑1 (MQ=25)

           aagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCgg                                                                          >  2:1356555/1‑93 (MQ=37)

             gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCg                                                                        >  1:4074328/1‑93 (MQ=37)

             gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCg                                                                        >  2:2989429/1‑93 (MQ=37)

              aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGa                                                                       <  2:1856431/93‑1 (MQ=37)

              aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGa                                                                       >  1:2896154/1‑93 (MQ=37)

               aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGAt                                                                      >  2:3571215/1‑93 (MQ=37)

                ttAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCTGCGATg                                                                     >  1:830001/1‑93 (MQ=37)

                ttAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATg                                                                     >  1:528330/1‑93 (MQ=37)

                ttAGAGAAGAAGAAAAATATGGGGAATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATg                                                                     >  2:3267301/1‑93 (MQ=25)

                 tAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATgg                                                                    <  1:3158658/93‑1 (MQ=37)

                 tAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATgg                                                                    <  1:2705302/93‑1 (MQ=37)

                  agcgaagtagaaAACTATGGGGAATTGAAGGGTGCTATGCCTGATGCCGATTAGTAGGCCGGATAAGGCGCTCGCG‑C‑CCAGTCCGGCGATggg                                                                   >  2:3758891/4‑92 (MQ=19)

                  agagAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGt                                                                   >  2:3758892/1‑93 (MQ=37)

                  agagAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGt                                                                   >  2:2838956/1‑93 (MQ=37)

                   gagaAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTg                                                                  >  1:4369816/1‑93 (MQ=37)

                      aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCAc                                                               >  1:2555288/1‑93 (MQ=37)

                      aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCAc                                                               >  2:3369546/1‑93 (MQ=37)

                       agaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACt                                                              <  1:277201/93‑1 (MQ=37)

                       agaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACt                                                              >  2:2132841/1‑93 (MQ=37)

                       agaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACt                                                              <  1:3966012/93‑1 (MQ=37)

                        gaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTg                                                             >  2:1965260/1‑93 (MQ=37)

                        gaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTg                                                             >  2:741458/1‑93 (MQ=37)

                         aagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGa                                                            <  1:1438459/93‑1 (MQ=37)

                          agaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGaa                                                           >  1:3242036/1‑93 (MQ=37)

                          agaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGaa                                                           >  1:3242038/1‑93 (MQ=37)

                           gaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAg                                                          <  2:2472126/93‑1 (MQ=37)

                            aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGc                                                         <  1:547422/93‑1 (MQ=37)

                             aaaaTATGGGGTATTGAGGGTTGCTGTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGcc                                                        <  2:2292360/93‑1 (MQ=25)

                             aaaaTATGGGGTATGGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGcc                                                        <  2:2292361/93‑1 (MQ=25)

                              aaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCt                                                       <  2:425608/93‑1 (MQ=37)

                                atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGa                                                     >  1:3772710/1‑93 (MQ=37)

                                atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGACAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGa                                                     >  1:3772709/1‑93 (MQ=37)

                                 tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGAt                                                    <  2:786810/93‑1 (MQ=37)

                                 tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGAt                                                    <  2:2988587/93‑1 (MQ=37)

                                 tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGAt                                                    <  2:1401802/93‑1 (MQ=37)

                                 tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGAt                                                    >  1:4170326/1‑93 (MQ=37)

                                  atGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATg                                                   <  1:532068/93‑1 (MQ=37)

                                  atGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATg                                                   <  2:1995073/93‑1 (MQ=37)

                                  atGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATg                                                   >  1:1105365/1‑93 (MQ=25)

                                   tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGc                                                  <  1:2061757/93‑1 (MQ=25)

                                     gggTATTGAGTTTTGCTATTCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGa                                                <  1:2257442/93‑1 (MQ=21)

                                     gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGa                                                <  1:2037490/93‑1 (MQ=25)

                                        tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                             >  1:3540096/1‑93 (MQ=25)

                                         aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                                            >  1:939264/1‑93 (MQ=25)

                                         aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                                            >  1:3986171/1‑93 (MQ=25)

                                          ttGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                           <  1:567366/93‑1 (MQ=25)

                                          ttGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                           >  1:286449/1‑93 (MQ=25)

                                            gAGGGTTTCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCG‑CTGCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAc                                         <  2:2648785/93‑1 (MQ=37)

                                             aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                                        >  1:892464/1‑93 (MQ=25)

                                              gggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcg                                       >  2:2621671/1‑93 (MQ=25)

                                                gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcg                                     <  2:165793/93‑1 (MQ=25)

                                                 ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                                    >  1:3709662/1‑93 (MQ=25)

                                                 ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                                    <  2:4309331/93‑1 (MQ=25)

                                                 ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                                    >  1:2570799/1‑93 (MQ=25)

                                                  tGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTc                                   >  1:2980682/1‑93 (MQ=25)

                                                  tGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTc                                   <  1:2450460/93‑1 (MQ=25)

                                                   tcTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCTCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                  <  1:1906558/92‑1 (MQ=25)

                                                   gCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                  <  1:1561812/93‑1 (MQ=25)

                                                    cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCtt                                 <  1:3204192/93‑1 (MQ=25)

                                                     tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                >  1:253266/1‑93 (MQ=25)

                                                     tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                >  1:4649264/1‑93 (MQ=21)

                                                        gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             >  2:4581493/1‑93 (MQ=21)

                                                        gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             >  2:2236947/1‑93 (MQ=21)

                                                         ccTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt                            <  1:2083303/93‑1 (MQ=21)

                                                          cTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                           >  1:1585595/1‑93 (MQ=21)

                                                          cTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                           >  2:2912015/1‑93 (MQ=21)

                                                           tGATGCCGATTCGTAGGCGGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc                          <  2:1015942/93‑1 (MQ=21)

                                                           tGATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc                          <  1:1773001/93‑1 (MQ=21)

                                                            gATGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                         <  2:3782634/93‑1 (MQ=21)

                                                             aTGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCt                        <  2:468634/93‑1 (MQ=21)

                                                             aTGCCGATTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCt                        <  2:1259964/93‑1 (MQ=21)

                                                                  gATTCGTAGGCGGGATAGGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTgg                   <  2:1085094/93‑1 (MQ=21)

                                                                   aTTCGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTggg                  <  1:130751/93‑1 (MQ=21)

                                                                      cGTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGc               <  2:4004477/93‑1 (MQ=21)

                                                                       gTAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCa              >  1:4335397/1‑93 (MQ=21)

                                                                        tAGGCGGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc             <  2:3555035/93‑1 (MQ=21)

                                                                        tAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc             <  1:3996655/93‑1 (MQ=21)

                                                                        tAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc             <  1:4079938/93‑1 (MQ=21)

                                                                        tAGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc             <  2:1644583/93‑1 (MQ=21)

                                                                         aGGCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACg            >  2:780588/1‑93 (MQ=21)

                                                                          ggCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGc           <  1:4531120/93‑1 (MQ=21)

                                                                           gCCGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCt          <  2:3733495/93‑1 (MQ=21)

                                                                             cGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCttt        <  1:1745065/93‑1 (MQ=21)

                                                                             cGGATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCttt        >  1:4247790/1‑93 (MQ=21)

                                                                             cGGATAAGGCCCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCttt        <  2:3615565/93‑1 (MQ=21)

                                                                              ggATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTa       <  1:126720/93‑1 (MQ=21)

                                                                              ggATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTa       <  2:1172559/93‑1 (MQ=25)

                                                                               gATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac      >  2:4277633/1‑93 (MQ=21)

                                                                               gATAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac      >  2:1717578/1‑93 (MQ=25)

                                                                                aTAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTaca     >  1:3571987/1‑93 (MQ=25)

                                                                                 tAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac    <  2:1619910/93‑1 (MQ=25)

                                                                                 tAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac    >  1:1299632/1‑93 (MQ=25)

                                                                                 tAAGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac    >  1:2408161/1‑93 (MQ=25)

                                                                                  aaGGCGCTCGCGCC‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACAcc   <  2:323595/93‑1 (MQ=25)

TACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCG‑G‑C‑GCA‑TCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCG  >  NC_000913/4448450‑4448623

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