BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,448,543	+C	intergenic (+48/‑90)	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,542	1	.	C	100.0%	455.7 / NA	161	intergenic (+47/‑91)	ytfQ/ytfR	putative LACI‑type transcriptional regulator/YtfR
Reads supporting (aligned to +/- strand): ref base . (0/0); new base C (108/53); total (108/53)

Only 100 of 176 total aligned reads displayed.

AAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCG‑‑G‑CGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCG  >  NC_000913/4448461‑4448623

aagaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCgg                                                                          >  2:305520/1‑93 (MQ=37)

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

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

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

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

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

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

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

   aaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGa                                                                       >  1:3357722/1‑93 (MQ=37)

    aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGag                                                                      >  2:92653/1‑92 (MQ=37)

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

    aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGAt                                                                      >  2:3668289/1‑93 (MQ=37)

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

      tAGAGAAGAAGAAAAAGATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATgg                                                                    <  2:3592331/93‑1 (MQ=25)

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

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

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

        gagaAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTg                                                                  <  1:1470951/93‑1 (MQ=37)

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

         agaagaaaaAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGc                                                                 <  2:4034139/93‑1 (MQ=25)

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

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

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

            agaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACt                                                              <  1:1845572/93‑1 (MQ=37)

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

               agaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGaa                                                           >  2:4561996/1‑93 (MQ=37)

               agaaAAAAGTGGGTTTTTGAGGGTTGCTTTGCAGTATGCCGGTTCGTGGCCGGGATAAGGCTCTCGCCC‑CCCATCCGGCGATGGTGCACTGaa                                                           <  2:2749707/93‑1 (MQ=255)

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

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

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

                    aaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTg                                                      >  2:3943339/1‑93 (MQ=37)

                     atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGTTGCACTGAAGCCTGa                                                     >  1:3141078/1‑93 (MQ=25)

                     atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                                     <  2:3828543/93‑1 (MQ=37)

                     atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                                     <  2:3722612/93‑1 (MQ=37)

                     atatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGa                                                     >  1:2585447/1‑93 (MQ=37)

                      tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    >  2:4382520/1‑93 (MQ=37)

                      tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    <  2:4129029/93‑1 (MQ=37)

                      tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    <  1:3673085/93‑1 (MQ=25)

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

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

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

                         ggggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCg                                                 >  2:1294656/1‑93 (MQ=25)

                          gggTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGa                                                <  1:3093006/93‑1 (MQ=25)

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

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

                               ttGAGGGTTGCTATGCCGGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                           <  1:3630552/93‑1 (MQ=25)

                                 gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAc                                         >  2:271434/1‑93 (MQ=25)

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

                                  aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                                        >  2:1187140/1‑93 (MQ=25)

                                   ggtttGCTATGCCTGATGCCGATTCGGAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcg                                       <  2:369737/90‑1 (MQ=21)

                                     gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcg                                     >  1:2059718/1‑93 (MQ=25)

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

                                       tGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGGACTGAAGCCTGATGCGTCGCTTACcgcggc                                   >  1:4006506/1‑91 (MQ=21)

                                       tGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTc                                   <  2:1508065/93‑1 (MQ=25)

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

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

                                          tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGCAGCCTGATGCGACGCTTACCGCGTCTTa                                >  1:2603567/1‑93 (MQ=21)

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

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

                                           aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCCGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt                               >  2:1452238/1‑93 (MQ=21)

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

                                             gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             <  1:4330273/93‑1 (MQ=21)

                                             gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             <  2:4662950/93‑1 (MQ=21)

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

                                             gCCTGAGGCCGATTCGTAGGCGGGATAGGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             <  1:4039304/93‑1 (MQ=21)

                                              ccTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt                            >  1:1410033/1‑93 (MQ=21)

                                              ccTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGACTGATGCGACGCTTACCGCGTCTTATCAt                            >  1:4390385/1‑93 (MQ=21)

                                              ccTGAGGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt                            <  1:3477592/93‑1 (MQ=21)

                                               cTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                           >  1:1070934/1‑93 (MQ=21)

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

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

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

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

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

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

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

                                                   tGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGTCTa                       >  1:2765624/1‑93 (MQ=21)

                                                   tGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa                       <  1:3865058/93‑1 (MQ=21)

                                                      cGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTg                    <  2:1363837/93‑1 (MQ=18)

                                                       gATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTgg                   >  2:1034431/1‑93 (MQ=21)

                                                       gATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTgg                   <  1:3451804/93‑1 (MQ=18)

                                                       gATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTgg                   >  1:2134313/1‑93 (MQ=21)

                                                       gATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTgg                   <  2:1948763/93‑1 (MQ=21)

                                                        aTTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTggg                  >  2:2299718/1‑93 (MQ=18)

                                                         ttCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                 >  1:2482458/1‑93 (MQ=21)

                                                         ttCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                 >  1:919694/1‑93 (MQ=21)

                                                         ttCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                 >  2:1938882/1‑93 (MQ=21)

                                                         ttCGTAGGCCGGATAAGGCCCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                 <  1:3202812/93‑1 (MQ=21)

                                                            gTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCa              <  1:3802300/93‑1 (MQ=21)

                                                            gTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCa              >  2:1800600/1‑93 (MQ=21)

                                                               ggCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGc           >  2:3249886/1‑93 (MQ=21)

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

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

                                                                    gATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGACTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTac      >  1:454447/1‑93 (MQ=21)

                                                                     aTAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTaca     <  2:2689588/93‑1 (MQ=25)

                                                                      tAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac    >  1:4344961/1‑93 (MQ=25)

                                                                        aGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCGTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGCGCACGCTTTACACCg  >  1:556692/1‑93 (MQ=25)

                                                                        aGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCg  >  2:4008804/1‑93 (MQ=25)

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

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

AAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCG‑‑G‑CGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCG  >  NC_000913/4448461‑4448623

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