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%	602.1 / NA	212	intergenic (+46/‑92)	ytfQ/ytfR	putative LACI‑type transcriptional regulator/YtfR
Reads supporting (aligned to +/- strand): ref base . (0/0); new base C (130/82); total (130/82)

Only 100 of 232 total aligned reads displayed.

ATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCG‑‑G‑CGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCG  >  NC_000913/4448449‑4448623

aTACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCg                                                                                       <  2:1285754/93‑1 (MQ=255)

  aCCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcg                                                                                     <  1:2257312/93‑3 (MQ=255)

    cGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcgcc                                                                                   >  2:4039379/1‑89 (MQ=255)

            aagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCgg                                                                          <  2:2620303/93‑1 (MQ=37)

            aagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCgg                                                                          <  1:2496790/93‑1 (MQ=37)

             agaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGc                                                                         <  2:1926489/93‑1 (MQ=37)

             agaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGc                                                                         >  2:3341962/1‑93 (MQ=37)

              gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCg                                                                        >  2:2810640/1‑93 (MQ=37)

              gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCg                                                                        >  1:835105/1‑93 (MQ=37)

               aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGa                                                                       <  1:1251979/93‑1 (MQ=37)

               aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGATTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGa                                                                       <  1:3645422/93‑1 (MQ=25)

                aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGAt                                                                      <  1:4916553/93‑1 (MQ=37)

                 ttAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATg                                                                     >  1:5694943/1‑93 (MQ=37)

                  tAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATgg                                                                    <  1:3883312/93‑1 (MQ=37)

                   agagAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATggg                                                                   >  2:6088448/1‑92 (MQ=37)

                   agagAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATggg                                                                   >  2:3301381/1‑92 (MQ=37)

                   agagAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGt                                                                   >  1:2920332/1‑93 (MQ=37)

                    gagaAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTg                                                                  >  2:1802245/1‑93 (MQ=37)

                    gagaAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTg                                                                  >  1:5971206/1‑93 (MQ=37)

                      gaagaagaaAAATATGGGTTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCa                                                                <  2:1454739/93‑1 (MQ=25)

                      gaagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCa                                                                >  2:1343554/1‑93 (MQ=37)

                       aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCAc                                                               >  2:2344801/1‑93 (MQ=37)

                       aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCAc                                                               <  1:5917917/93‑1 (MQ=37)

                       aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCAc                                                               >  2:5307301/1‑93 (MQ=37)

                       aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCAc                                                               >  2:5548940/1‑93 (MQ=37)

                        agaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACt                                                              >  2:5894308/1‑93 (MQ=37)

                         gaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTg                                                             >  2:2028768/1‑93 (MQ=37)

                         gaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTg                                                             >  1:3672460/1‑93 (MQ=37)

                         gaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTg                                                             <  1:2363741/93‑1 (MQ=37)

                          aagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGa                                                            <  1:6186135/93‑1 (MQ=37)

                             aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGCAGc                                                         >  2:3163104/1‑93 (MQ=25)

                             aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGc                                                         <  2:1741064/93‑1 (MQ=37)

                              aaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGcc                                                        <  1:4011937/93‑1 (MQ=37)

                               aaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGccg                                                       >  2:1212983/1‑92 (MQ=37)

                               aaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCt                                                       <  2:4649313/93‑1 (MQ=37)

                               aaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCt                                                       <  2:5111350/93‑1 (MQ=37)

                                  tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    >  1:193551/1‑93 (MQ=25)

                                  gatGGGGTTGTGAGGGTTGTTATGCTTGATCCCTTTTGGTGGGCGGGATAGGGCGCCCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    <  2:4355407/92‑1 (MQ=9)

                                   atGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATg                                                   >  1:2784107/1‑93 (MQ=25)

                                    tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGc                                                  >  2:3280858/1‑93 (MQ=25)

                                    tGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCAGGCGATGGTGCACTGAAGCCTGATGc                                                  >  2:1601616/1‑93 (MQ=37)

                                      ggttATGGAGGTTTGCTATGCTGGATGCCGATTCGTAGGCGGGATAGGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGa                                                <  2:4473549/90‑1 (MQ=18)

                                       ggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGAc                                               <  1:3606299/93‑1 (MQ=25)

                                        gTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACg                                              >  1:4594421/1‑93 (MQ=25)

                                        gTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACg                                              >  1:3486385/1‑93 (MQ=25)

                                         tATTGAGGTTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                             <  2:4562120/93‑1 (MQ=25)

                                         tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                             >  2:5147601/1‑93 (MQ=25)

                                         tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCACCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                             <  2:2980254/93‑1 (MQ=25)

                                          aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                                            >  1:5543690/1‑93 (MQ=25)

                                          aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                                            >  2:3531564/1‑93 (MQ=25)

                                          aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                                            >  2:1501486/1‑93 (MQ=25)

                                          aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                                            >  2:4753793/1‑93 (MQ=25)

                                           ttGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                           >  1:2636221/1‑93 (MQ=25)

                                            tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTa                                          <  2:2750196/93‑1 (MQ=25)

                                            tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTa                                          >  2:813434/1‑93 (MQ=25)

                                             gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAc                                         <  2:2879011/93‑1 (MQ=25)

                                             gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAc                                         <  2:957621/93‑1 (MQ=25)

                                             gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAc                                         <  2:5609164/93‑1 (MQ=25)

                                               gggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcg                                       >  2:3775224/1‑93 (MQ=25)

                                               gggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcg                                       >  1:5114965/1‑93 (MQ=25)

                                                ggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgc                                      >  2:2432203/1‑93 (MQ=25)

                                                ggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgc                                      <  1:5285114/93‑1 (MQ=25)

                                                 gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcg                                     <  1:2732884/93‑1 (MQ=25)

                                                  ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                                    >  1:44243/1‑93 (MQ=25)

                                                   tGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTc                                   >  1:4700019/1‑93 (MQ=25)

                                                   tGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTc                                   >  2:3165379/1‑93 (MQ=25)

                                                    gCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                  >  2:3833523/1‑93 (MQ=25)

                                                     cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCtt                                 >  2:668663/1‑93 (MQ=25)

                                                      tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                >  2:3023482/1‑93 (MQ=21)

                                                      tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                <  2:1966195/93‑1 (MQ=25)

                                                       aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt                               <  2:4946849/93‑1 (MQ=21)

                                                       aTGCCTGATGCCGATTCGGAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt                               <  2:2489476/93‑1 (MQ=21)

                                                        tGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATc                              >  1:3307348/1‑93 (MQ=21)

                                                        tGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATc                              >  1:2078023/1‑93 (MQ=21)

                                                         gCCTGATGCGGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             <  2:5683565/93‑1 (MQ=21)

                                                         gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             >  1:5068824/1‑93 (MQ=21)

                                                         gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             >  1:1416444/1‑93 (MQ=21)

                                                         ccgtGATGCCGATTCGTAGCCCGAATAGGGCGCTCGCCC‑CGCATCGGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             <  2:1137364/90‑1 (MQ=21)

                                                         cccTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             <  1:3355112/92‑1 (MQ=21)

                                                           cTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                           <  1:5652262/93‑1 (MQ=25)

                                                            tGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc                          <  1:2830441/93‑1 (MQ=21)

                                                            tGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc                          >  1:1056853/1‑93 (MQ=21)

                                                            tGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc                          >  1:2587311/1‑93 (MQ=21)

                                                             gATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                         <  1:1639966/93‑1 (MQ=21)

                                                             gATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                         <  1:180224/93‑1 (MQ=21)

                                                             gATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                         >  1:5701110/1‑93 (MQ=21)

                                                              aTGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCt                        >  1:5654054/1‑93 (MQ=21)

                                                                gCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTAc                      <  1:648662/93‑1 (MQ=21)

                                                                 ccGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACt                     <  2:1691367/93‑1 (MQ=21)

                                                                 ccGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACt                     >  1:1766225/1‑93 (MQ=21)

                                                                     ttCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                 >  2:6020026/1‑93 (MQ=18)

                                                                     ttCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAACATGCCTACTGGGa                 >  2:3215399/1‑93 (MQ=18)

                                                                      tCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAg                >  1:456104/1‑93 (MQ=21)

                                                                        gTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGcc              >  1:3410956/1‑92 (MQ=21)

                                                                         tAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc             <  1:5760268/93‑1 (MQ=21)

                                                                             ccGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCtt         >  1:1226311/1‑93 (MQ=21)

                                                                               ggATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTa       >  2:4592307/1‑93 (MQ=21)

                                                                                gATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac      <  2:4269530/93‑1 (MQ=21)

                                                                                 aTAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTaca     >  1:673118/1‑93 (MQ=25)

                                                                                    aGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCg  <  2:4982862/93‑1 (MQ=25)

ATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCG‑‑G‑CGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCG  >  NC_000913/4448449‑4448623

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